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Preface 



This textbook has been compiled from a set of lecture notes developed while 
teaching courses in the theory of structures to civil engineering students at Addis 
Ababa University during the past seven years. The book is primarily intended for 
use as a text for instruction and contains sufficient material for a two-semester 
course in theory of structures. It may also be useful to the structural engineer 
who wishes to strengthen his background in structural mechanics. 

The purpose of this book is to present a balanced treatment of the funda- 
mental principles of structural mechanics, with their applications to the analysis 
of structural systems and their components. The coverage is selective, to allow a 
thorough treatment of the most common and useful analytical methods of 
structural analysis. 

An attempt is made to present the subject matter in a -unified, coherent and 
easy-to-understand manner which brings out the basic principles underlying the 
field of structural theory. The book is illustrated with ample example problems, 
to which solutions are presented to demonstrate the various methods, and also 
to widen the scope of the subject covered by the text. 

The author is indebted to the authors of the many books he has freely 
consulted in the preparation of this work. The author also wishes to acknowledge 
his debt to all his students who, over the years, checked out the examples and 
assignment problems. 

NEGUSSIE TEBEDGE 

Addis Ababa 
June, 1982 
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1. Introduction 



1.1 STRUCTURAL ANALYSIS 

Structural analysis is the process of determining the response of a structure due 
to specified loadings in order to satisfy essential requirements of function, safety, 
economy and sometimes aesthetics. This response is usually measured by calculat- 
ing the reactions, internal forces of members, and displacements of the structures. 

Structures may be classified into two general categories: statically determinate 
and statically indeterminate. A structure which can be completely analysed by 
means of statics alone is called statically determinate. It then follows that a 
statically indeterminate structure is one which cannot be analysed by means of 
statics alone. 

There are specific advantages and disadvantages in using one type of structure 
over the other. The primary advantage of a statically indeterminate structure is 
that it will generally have lower bending moment and shear force than a 
comparable determinate structure. Another advantage of a statically indeterminate 
structure is that it is generally stiffer for a given weight of material than a 
comparable determinate structure. Both of these advantages are a result of 
continuity of structural members acting to reduce stress intensities and displace- 
ments. A statically indeterminate structure can often furnish a compensation by 
redistribution within the structure in the case of overloads. On the other hand, 
however, indeterminate structures introduce computational difficulty in 
establishing the required equations. Another disadvantage is that indeterminate 
structures are, in normal cases, internally stressed due to differential settlement 
of supports, temperature changes and errors in the fabrication of members, 

1.2 STATICAL INDETERMINACY 

Consider a structure in space subjected to non-coplanar system forces. For the 
structure to be in equilibrium, the components of the resultants in the three 
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orthogonal directions must vanish. This condition constitutes the six equations of 
equilibrium in space which are written as 

ZF X = 0 2F y = 0 2F 2 = 0 

[1.1] 

XM X = 0 XM y = 0 XM Z = 0 

For a structure subjected to a coplanar force system, only three of the six 
equations of equilibrium are applicable. The three equations of equilibrium in 
the xy plane are 

XF X =Q 

2F r =0 [1.2] 
XM z =0 

When a structure is in equilibrium, each member, joint, or segment of the 
structure must also be in equilibrium and the equations of equilibrium must also 
be satisfied. As discussed earlier, a structure which can be analysed by means of 
the equations of equilibrium alone is statically determinate. This book deals with 
statically indeterminate structures in which the structures cannot be analysed by 
the equations of equilibrium alone. 

When a structure is statically indeterminate, there is some freedom of choice 
in selecting the member or reaction to be regarded as redundant. When the 
reaction is taken as the redundant, the structure is said to be externally 
indeterminate. On the other hand, when the member itself is regarded as the 
redundant, the structure is said to be internally indeterminate. It is also possible 
that the structure may have a combination of external and internal indeterminacy. 

The question of identifying external or internal indeterminacy is largely of 
academic interest. What is of primary importance in the analysis of indeterminate 
structures is to know the degree of total indeterminacy. Nevertheless, a separate 
discussion of external and internal indeterminacy is desirable as a method to 
evaluate the degree of total indeterminacy. 

(a) External Indeterminacy If the total number of reactions in a structure 
exceeds the number of the equations of equilibrium applicable to the structure, 
the structure is said to be externally indeterminate. The structures shown in 
Fig. 1.1 are examples of external indeterminacy. Each of the structures has five 
reaction components. Since there are only three equations of equilibrium, there 
are two extra reaction components that cannot be determined by statics. The 
number of unknown reactions in excess of the applicable equations of 
equilibrium defines the degree of indeterminacy . Thus the structures of Fig. 1.1 
are indeterminate to the second degree. An alternative approach to determine the 

degree of indeterminacy would be to remove selected redundant reactions until 
the structure is reduced to a statically determinate and stable base or primary 
structure. 
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(0 

Figure 1.1 



(b) Internal Indeterminacy A structure is internally indeterminate when it is 
not possible to determine all internal forces by using the three equations of 
static equilibrium. For the great majority of structures, the equation of whether 
or not they are indeterminate can be decided by inspection. For certain 
structures this is not so, and for these types rules have to be established. The 
internal indeterminacy of trusses will be first considered, and then that of 
continuous frames. 

It is evident that any truss developed by using three bars connected at three 
joints in the form of a hinged triangle, and then using two bars to connect each 
additional joint, forms a stable and determinate truss. This is because the shape 
of the triangle cannot be changed without changing the length of any of the 
members. For stable and determinate trusses, built up as an assemblage of 
triangles, there are two conditions of equilibrium for each joint, so that if there 
are/ joints, m members and r reaction components, a test for statical determinacy 
is: 

2/ = m + r [1.3] 

In this equation, the left-hand side represents the total possible number of 
equations of equilibrium, while the right-hand side represents the total number 
of unknown forces. 

The above equation is usually written in the form 

m = 2j-r [1.4] 
If there are more members than are indicated by the equation, then the 



3 



METHODS OF STRUCTURAL ANALYSIS 



structure is statically indeterminate; whereas if it has fewer members it is 
unstable. Caution must be exercised in applying the above equation because of 
the fact that the fulfilment of this equation is a necessary condition but not 
sufficient for internal stability of trusses. This may be summarised as 

m = 2/ -r (determinate if stable) 

m > 2j - r (indeterminate if stable) 

m < 2/ — r (unstable) 

The truss in Fig. 1.2(a) has m = 17 J = 10 and r = 3. Application of [1.4] gives 

(10 x 2) — 3 = 17 members, thus the structure is statically determinate. Referring 
to Fig. 1.2(b), there are 18 members, or one more member than is needed for a 
determinate structure; thus, the additional diagonal member is redundant and 
the truss is indeterminate to the first degree. Figure 1 .2(c) represents the omission 
of one diagonal member, keeping the same total number of bars, m = 17. Again 
the condition equation is satisfied. However, inspection of the truss indicates 
that the structure is unstable with one panel free to collapse, thus causing the 
entire truss to collapse. Hence, satisfaction of the above equation is not a 
sufficient condition for internal stability of trusses. Inspection of the structure 
and consideration of stress paths are more reliable approaches to settle the 
question of stability and internal indeterminateness of trusses. 

An alternative approach to determine the degree of indeterminacy of trusses 




Figure 1.2 
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is by removing the redundant quantities until a determinate and stable base 
structure remains. 

The number of rigidly jointed frames are subject to shearing forces, bending 
moment and axial force, so that there are three unknown internal forces for each 
member, or a total of 3m unknown components. Moreover, at each joint three 
equations of equilibrium can be written, giving 3/ equations in all. Therefore for 
a statical determinacy, it is necessary that 

3/ = 3w+r [1.5] 

or that the number of redundants n is given by 

n = 3m + r - 3/ 

When there is a roller or pin support, the degree of indeterminacy is reduced 
by one or two, respectively, for each support. 

An alternative approach, which in this case may be considered more 
instructive, is the method by inspection where the structure is cut until it 
becomes a determinate and stable base structure. Consequently, the total number 
of released internal force components corresponds to the degree of mdeterminacy. 



1.3 KINEMATIC INDETERMINACY 

When a structure is subjected to a system of forces, the overall behaviour of the 
members of the structure may be defined by the displacement of the joints. The 
joints undergo displacements in the form of translation and rotation. A system 
of joint displacements is known to be independent if each displacement can be 
varied arbitrarily and independently of the other displacements. The number of 
independent joint displacements that serve to describe all possible displacements 
of a structure is known as the number of degrees of freedom or degree of 
kinematic indeterminacy. 

In determining the degree of kinematic indeterminacy, attention is focused 
on the number of independent displacement degrees of freedom that the 
structure possesses. If a structure has n degrees of freedom, that is, n number of 
independent displacement quantities required to describe all possible displace- 
ments for any loading condition, the structure is said to be kinematically 
indeterminate to the «th degree. When these displacements are set to zero, the 
structure then becomes kinematically determinate. 

Consider, for example, the rigid-jointed plane frame shown in Fig. 1 .3, which 
is fixed at supports A and C and has a hinged support at D. Assuming that the 
axial deformations are negligible, there will be no axial displacements in the 
frame and the only unknown displacements are the joint rotations 0 B and 0 D at 
joints B and D, respectively. Since these displacements are independent of one 
another, the degree of kinematic indeterminacy of this structure is two. 
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It is observed that the degree of statical indeterminacy of the frame of Fig. 1.3 
is four since there are a total of seven possible unknown reactions and three 
equations of equilibrium. If, for instance, the fixed support at C is replaced by a 
hinge, the degree of statical indeterminacy is reduced to three since an additional 
equilibrium condition is introduced. However, the kinematic indeterminacy is 
increased by one since an independent rotation at C now becomes possible. In 
general, an introduction of a displacement release decreases the statical 
indeterminacy and increases the kinematic indeterminacy. 



The objective of structural analysis is to study the response of a structure to 
specified loadings after determining the external reactions and internal stress 
resultants. The forces determined must satisfy the conditions of equilibrium and 
the displacements produced by these forces must be compatible with the 
continuity of the structure and the support conditions. In determining the 
unknown forces in a statically indeterminate structure, the equations of 
equilibrium are not sufficient, and additional equations must be formulated 
based on compatibility of displacements. These supplementary equations that 
ensure the compatibility of the displacements with the geometry of the structure 
are known as the compatibility conditions. 

Two general methods of analysis are available for the solution of statically 
indeterminate structures. The first is the force or flexibility method. This 
method is simple and conceptually straightforward to understand and provides 
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Figure 1.3 



1.4 METHODS OF STRUCTURAL ANALYSIS 
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an effective method for certain types of structures. In this method the structure 
is made statically determinate by providing a sufficient number of releases by 
removing the redundant forces. Due to the given loading condition the primary 
structure undergoes inconsistency in geometry which must then be corrected by 
applying the redundant forces such that compatibility conditions throughout 
the structure are established. This method is sometimes referred to as the 
compatibility method. 

The second method of analysis of statically indeterminate structures is the 
displacement or stiffness method. This method is also simple and straightforward 
and provides an effective method for certain classes of structure. In this method, 
restraints are imposed to prevent displacement of joints until the structure 
becomes kinematically determinate and the forces required to produce the 
restraints are evaluated. Displacements are then permitted to take place at the 
restrained joints until the imposed restraining forces have been removed such 
that equilibrium conditions throughout the structure are established. This method 
is also known as the equilibrium method. 

Either the force or the displacement method can be used to analyse any 
structure. The choice of the method of analysis, either force or displacement, 
depends largely on the degree of statical or kinematic indeterminacy. In both 
methods, the analysis generally involves the solution of a system of simultaneous 
equations where the number of unknown variables must be equal to the degree 
of indeterminacy. If manual calculations are to be adopted, it would be logical 
to use the method that produces the smaller set of simultaneous equations. 

1.5 PROBLEMS 

Determine the degree of statical indeterminacy of the structures shown below. 



1.1 



1.2 
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INTRODUCTION 

Determine the degree of statical and kinematical indeterminacy of the structures 
shown below. 
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1.13 
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1.14 
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2. Methods of 

Consistent Displacements 



2.1 INTRODUCTION 

For statically indeterminate structures there will be an indefinite number of 
combinations of redundant forces which will satisfy equilibrium conditions. 
However, among them there will be only one set of values that will simultaneously 
satisfy the requirements of equilibrium and compatibility. Compatibility places 
constraints on the displacements of a structure to ensure continuity and that the 
structure conforms to the displacement boundary conditions prescribed by the 
supports. 

The methods of consistent displacement are based on the concept of 
equilibrium of forces and compatibility of displacements which may be stated as 
follows: Given a set of forces applied on a statically indeterminate structure, the 
reactions must assume such values that satisfy not only the conditions of static 
equilibrium with the applied loads but also the conditions of compatibility. The 
general method of consistent deformation is applicable for analysing all types of 
indeterminate structures. It is also applicable whether the structure is subjected 
to external loading, temperature changes, movements of supports, fabrication 
errors, or any other cause. Of course, there are other methods that are definitely 
superior for certain specific structures or loading conditions, but methods of 
consistent deformation are the most versatile and general. 

2.2 ANALYSIS OF BEAMS 

The principle of consistent displacement can best be illustrated by considering 
singly indeterminate structures. As a simple and classic example of this method 
consider a propped cantilever beam as shown in Fig. 2. 1 . The beam has three 
unknown reactions K A , M A and F B and is therefore statically indeterminate to 
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(a) Actual structure 




B 



V, 



B 



W 



(b) Base structure 



m 



(c) Displacement due to known 
loads 



A 



bo 




(e) 

Figure 2. 1 



the first degree. Any one of the unknown reactions may be taken as the 
redundant. A stable and determinate primary structure may be formed by 
determinate primary structure by selecting as the redundant the vertical reaction 
at the right support, as shown in Fig. 2.1(b). 

The displacement of the cantilever beam AB may be considered to consist of 
the superposition of two independent displacements: 

Afco ~ upward deflection at B of the base structure due to the known 
applied loads only 
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A bb = upward deflection at B of the base structure due to the redundant 

It may be noted that it is not possible to evaluate A bb prior to the evaluation 
of V B . However, by applying the principle of superposition such that A bb = 8 bb V B , 
where 

8 bb = upward deflection at B of the base structure due to a unit upward 
load at B 

then, the condition that the support at B is rigid requires its displacement, Ab , 
the algebraic sum of displacements due to the applied loads and the redundant, 
must be zero. This geometric condition, defined as the equation of consistent 
deformation, may be writtten as 

A B =5 6o +F B 5 b6 = 0 [2.1] 



or 



Mmdx 



8 bb [m 2 dx 



[2.2] 



FA 



where M is the moment in the base structure due to the applied loads and m is 
the moment due to a unit load acting at B. 

It is noted that if Kb acts in the same direction as A bb , a negative value is 
obtained which indicates that the assumed direction is wrong. Conversely, a 
positive value for V B indicates that the assumed direction is correct. In general, 
it must be noted that the magnitude of the true reaction K B is that required to 
restore the end B of the beam to its original position level with A. 

In a similar manner, if Af A which is the moment reaction at A is taken as the 
redundant, the applied loads will cause the tangent at A to rotate through an 
angle 6 ao . If the rotation due to a unit moment at A is taken as 9 aa the moment 
M A necessary to rotate the tangfent at A to the original horizontal position is 



Ma = 



Mm Q dx 

9 -<oJ f f23l 
O aa i mldx 12 31 

EI 



where m 0 is the moment due to a unit moment acting at A. 

The analysis of beams of higher degree of indeterminacy follows closely the 
procedure described above. For a beam with n degrees of indeterminacy, n 
redundants are selected which will be removed from the structure and replaced 
by n effectively equivalent redundant forces X X ,X 2 , • ■ ■ X n . All these redundant 



13 



METHODS OF STRUCTURAL ANALYSIS 



forces and the given external loads are applied on the base structure such that 
their magnitudes must cause the displacements at the points of application of 
the n redundants of the base structure to be equal to the displacement of the 
corresponding points on the actual structure. 

Consider the four-span continuous beam of Fig. 2.2. The beam has three 
redundant reactions which can be chosen in a variety of ways, one of which is 
shown in Fig. 2.2(a). 

At this stage it is convenient to follow a definite notation for the various 
redundant forces and displacements. The redundant forces X a , X b and X c are 
recognised and identified by single subscripts which denote their point of 
application. The displacements are identified by double subscripts: the first 
subscript denotes the point on the base structure at which the displacement 
occurs, and the second subscript is used to denote the force producing the 
displacement. If, for example, the points A, B, C, etc. are the points on the base 
structure where the redundants occur, then, 

X a = the redundant force at point A 

A ao = the displacement in the base structure at point A in the direction of 
X ai caused by the actual applied loads acting on the structure 

8 aa = displacement in the direction of X a in the base structure caused by 
X a = 1 and no other load acting 

&ab ~ displacement in the base structure at A in the direction of X a caused 



8 ac = displacement in the base structure at A in the direction of X a caused 
by X c = 1 acting alone 

Since the displacements at A, B, and C should be zero, the reactions X ai X b 
and X c must have values such that compatibility condition is satisfied. Thus, 
using the above notation in the superposition equations, which gives as many 
equations as there are redundants, the equations may be written as follows: 

&ao + X a b aa + X b 8 ab + X c 8 ac = 0 



Since 8 ab = 8 ba , 8 ac = 5 ca , etc. by Maxwell's principle of reciprocal deflections, 
[2.4] may be written as 



by X b = 1 acting alone 



&bo + Xa&ba + X b 8 bb + X c 8 bc 



0 



0 



[2.4] 




[2.5] 
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1 A|H||| IB 



n 




<-2 




c 



(a) Actual structure 



: 



TTTTT 



a jX b 

(b) Base structure 



LLLTJ 



&00 



1 



1 



1 



P 3 



^bo 



(c) Displacement due to applied load 



(d) Due lo X e = 1 



>ab 



'bb 



'cb 



(e) Due to X b = 1 



(f) Due to X c = 1 






Figure 2.2 
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The equations can be written in the following matrix form: 









+ 


A co . 





5 a6 
5 ac 



5 a6 
5 66 
5 6c 



"ac 
5 6c 
See 









0 




*6 




0 








0 



EXAMPLE 2 . 1 Determine the reactions and support momen t of the 
continuous beam shown in Fig. 2.3. 



A 



El = Cons font 



B 



-75 kN m 




5.0m 




(a) Actual structure 




(b) Base structure 



(c) M 0 - diagram 



DkN 



2.5m 



5.0m 



IOkN 




lOkN 



C 





+ 5kN m 



(d) m, - diagram 




+ /OkN m 



(e) m 2 - diagram 



Figure 2.3 
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The beam is indeterminate to the second degree, and the redundants chosen 
are the reactions at B and C. The moment diagrams due to the applied load 
Xi = 1 and * 2 = 1 are shown in Fig. 2.3(d) and (e) respectively. 

The elastic equations are 



The displacements are obtained by graphic multiplication method: 



Substituting the 5 terms into the elastic equations: 

-729.17+ 41.67*! + 104.17*2 =0 
-2109.38 + 104.17*, + 333.33* 2 =0 

The solution of the simultaneous equations is 

*, = 7.68 kN 
* 2 = 3.93 kN 

Thus, 

R B =7.68 kN (upward) 
R c = 3.93 kN (upward) 

From statics 

R A = 1.61 kN (downward) 
M A = 2.70 kN (clockwise) 



A 10 +*i5u +* 2 5 12 - 0 

^20 + J2 + *2^22 = 0 
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Alternative Solution 




The moment reaction and the vertical reaction at A are chosen as redundants. 
The moment diagram due to the applied loads *, = 1 and * 2 = 1 are shown in 
Fig. 2.4(d) and (e) respectively. 

The displacements are 

/2.5x2.5\/2 , \ /I2.5x2.5\ / 1 „ 

= 78.13 

/5.0x5.0\ (2 

= 83.33 

««-(-=2f=S)(U9-(=2f=5) fV*) 

= -20.83 



-_ ^ 



,._.5x2.5\/2 2.5\ /I2.5x2.5\/ , 2 
= - l ( J x - - ( ) I 2.5 + - 




= -15.63 



EJ8 22 = (1.0x5.0) (1.0)+ ^— 



x 5.0 \ / 2 




2 / \3 



= 6.67 

Substituting into the elastic equations: 
78.13 + 83.33*, - 20.83*2 = 0 

-15.63 - 20.83*, + 6.67* 2 = 0 
The solution of the simultaneous equations is 

*, =-1.61 kN 

* 2 = -2.70 kN m 
The reactions are 

/? A = 1.61 kN (upward) 
M A = 2.70 kN m (clockwise) 
R B = 7.68 kN (upward) 



R c = 3.93 kN (upward) 
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A 



EI = Constant 



B 




5.0m 



(a) Actual structure 



2.5m 




5.0m 




Xo = / 




El = Constant 




(b) Base structure 



IOkN 




(c) Af 0 - diagram 



f5.0kN 



(d) - diagram 



IQkN 




12.5 kN m 




-1.0 kN m 




(-) 



-/.0 




x 2 



(e) m 2 - diagram 

Figure 2.4 
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2.3 ANALYSIS OF TRUSSES 



A statically indeterminate truss with external redundant reaction or internal 
redundant member may be analysed by a procedure closely analogous to that 
followed in beams. The analysis of trusses with a redundant reaction consists 
of choosing a base structure by removing the redundant reactions. Acting on 
this base structure are the applied loading and the redundant reactions. Then 
the condition of compatibility is applied such that the displacements in the 
direction of the redundants become zero. In a similar manner, when the truss 
has redundant members, the base structure is obtained by cutting the redundant 
members and replacing it by a pair of forces and then applying the condition of 
compatibility. Take, for example, the truss shown in Fig. 2.5. The truss is 
internally indeterminate to the first degree. 




Figure 2.5 



In this truss, any member may be considered redundant. Choosing member 
AB as the redundant, the redundant member is removed by cutting it at any 
section. Due to the effect of the external load/* on the base structure, the two 
faces of the cut member AB will be displaced by A ao . Now, applying a pair of 
forces N AB as shown in Fig. 2.5(b), such that the relative displacement of the 
actual truss at the cut surface is zero, gives the following relationship: 

A«,+JV A bS« = 0 [2.7] 

where 5^ is the relative displacement of the cut faces due toN AB = 1. The 
internal force in the redundant member is 



N AB = ~ T— 

°aa 



[2.8] 
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But from virtual work principle, 

NnL 



EA 

n lk 

EA 



[2.9] 



where N = force in any member due to the external applied load acting on the 

base structure 

n = force in any member due to a unit pair of forces applied at the cut 
faces of the member. 



Thus 



7Va R = 



„NnL 



AB " nH 



[2.10] 



2 



EA 



Note that the summation in the denominator is taken over the whole truss, and 
the summation in the numerator applies only on the base structure. 

The analysis of trusses of higher degree of indeterminacy follows closely the 
procedure described above. Consider, for example, the truss shown in Fig. 2.6 
which is externally statically indeterminate to the second degree. If the supports 
at B and C are removed, a simple truss supported at A and D Will be the basic 
determinate truss. The deflected bottom chord due to the applied loading is 
shown in Fig. 2.6(b). The displacements at B and C are determined from the 
expressions 

Afto_S EA 

l2 - ,1] 

Figure 2.6(b) shows the displacements at B and C due to a unit load applied at 
B, and in Fig. 2.6(c) due to a unit load applied at C. The vertical displacements 
are determined from the expressions 

dbb -*EA 

8 bc =8 cb = z'^ [2.12] 
bcc = *EA 
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(a) 




(b) 




(0 




(d) 



Figure 2.6 



In the above expressions N stands for forces in the members due to the external 
applied loads on the base structure, and N b and N c are the forces in the members 
due to a unit, load applied at B and C, respectively. 

The conditions of compatibility required from which R% and R c can be 
determined are 



A co +fl b 5 fec +fl c 5 cc =0 



[2.13] 
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Figure 2.7 



Consider again the truss shown in Fig, 2.6. If the redundants are taken as the 
bar forces X\ and X 2 shown in Fig. 2.7(b), then the determinate truss is three 
independent simple span trusses. Due to the effect of the external loading on the 
base structure, the two faces of the cut members 1 and 2 will be displaced by 
A 10 and A 2 o» respectively. After applying a pair of forces F x = 1 and F 2 = 1 as 
shown in Fig. 2.7(c) and (d), the corresponding relative displacements of the cut 
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faces can be determined. The displacements are 

Nn*L 

A20 " 2 EA 

n\L 

«12 -«21-Z ^ 



[2.14] 



The conditions of compatibility from which X x and AT 2 can be determined are 

Aio + *1 5 11 + * 2 5l2 =0 

10 1 11 2 1 2 [2.151 

^20 + ^1^12 +^2^22 = 0 



EXAMPLE 2 .2 Find the reaction at B and the bar force in member BF, of the 
truss in Fig. 2.8. The cross-sectional area of the members in cm 2 are shown in 
parentheses. E is constant. 

The given truss is indeterminate to the second degree; it has one redundant 
member (internal indeterminacy) and one redundant reaction (external 
indeterminacy). 

A base structure is obtained by removing the reaction at B and cutting the 
diagonal member BF. The two conditions of compatibility are: 

A B +fl B 5 bf) +F BF 5 v =0 
A F +R B 8 bf +F BF 8 ff =0 

The displacements are computed in tabular form as shown in Table 2.1. Substitut- 
ing the displacements: 

-1680.8+ 134.9 \R B + 54.52F BF =0 
- 202.3 + 54.52K B + 118.1F BF =0 
The solution of the simultaneous equations is 

R B = 14.44 tonnes (upward) 



F BF = —4.95 tonnes (compression) 
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(b) Base structure 




(c) Due to external redundant R B = 1 



-0.80 




Figure 2.8 
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2,4 ANALYSIS OF FRAMES 

A framed structure is composed of an interconnected assemblage of beams and 
columns. A frame is said to be rigid if the members are rigidly connected. The 
basic analysis of statically indeterminate frames by the method of consistent 
deformation is essentially an extension of the same principle encountered in 
dealing with beams. 

The members in frames are usually subjected to both axial and bending 
stresses; however, the axial stresses in the members of rigid frames are in most 
cases small compared with that of bending stresses. Thus, in computing the 
displacements in rigid frames for the conditions of consistent deformation, the 
effects of the axial stresses are usually neglected and the effects of bending 
stresses only are considered. This, however, does not mean that there are no 
axial forces in the members even if the change in the length of the members of 
rigid frames has insignificant effect on the values of the redundants. 

To formulate the equations for the general case of multiply redundant 
structures, consider the frame shown in Fig. 2.9, which is triply statically 
indeterminate. Let the three support reactions at A be chosen as the redundants. 
When these redundants are removed, A will be displaced vertically and 
horizontally and will also rotate. 

It will be seen that it will be convenient to adopt a slightly different notation 
with numerical subscripts for the redundants and displacements, which are 

defined as 

A 10 , A 20 , A 30 = displacements at A in the directions of X\ 9 X2 and A3 

respectively, due to the applied loads on the base structure 

5 1 1» ^21* S31 = displacements at A on the base structure in the directions 

of X x , X 2 and X 3 respectively, due to X x = 1 acting alone 

5 12 , 622* 5 32 = above displacements on the base structure due to 

X 2 = 1 acting alone 

6 13 , 5 2 3, 633 = the above displacements on the base structure due to 

X 3 = 1 acting alone 

If it is known that there are no support displacements, the equations of 
consistent deformation are 

A,o+^i5ii+^25i2+^ 3 5i3 = 0 

A 20 +X 2 5 12 +* 2 6 22 + X 3 6 23 = 0 [2.16] 
A20 + X 3 8 l3 + A r 2 5 23 +X 3 5 33 = 0 

The general equation for a structure with n redundants may then be written in 
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Figure 2.9 
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matrix form as 







r«n+«i2+-- 


. + «!„ 




*i 




0 


- • | 


+ 


• ■ • • • 


* 




X 2 

# * a 




0 

# ■ ■ 


. A *0. 




Sin + &2n + ■ ■ 


• + 5„„_ 




x„. 




0 



[2.17] 



These equations, sometimes referred to as the elastic equations, form the basis 
for several different methods of analysing statically indeterminate structures. 
The coefficients 6 { x , 6 12 - - - of the redundants on the base structures, which 
are the displacements due to unit loads are known as flexibility coefficients or 
influence coefficients. 

Whenever support displacements occur, the right-hand side of the equations 
may be suitably adjusted before solving the simultaneous equations. 

In the general case where deflections occur as a consequence of flexural and 
axial deformation of members of the structure, displacements in the base 
structure due to the applied loads may be written in the form 



I Mmdx 
A>o=j— + Z 



NnL 
EA 



[2.18] 



and the flexibility coefficients are 
5 *-J El + 1 ^EA 



[2.19] 



EXAMPLE 2.3 Determine the reaction components at support D of the frame 
of Fig. 110. 

Since the horizontal and vertical displacements and the rotations at support D 
must be zero, the compatibility equations are 



[Aio 






+ 


A30. 





oil 

621 
831 



5 12 5 13 
5 32 6 33 









0 




x 2 




0 




*3. 







The displacement coefficients are evaluated using the graphic multiplication 
method: 



EIA 



10 



.(_M 



X 1 M (2.5) + (-5.0 x 6.0)(3) + I (3 0) 



= -136.125 
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EIA 20 = I- 30 * 15 1 (-3) + (0.5 x 6X^-0.5) + ^^~ (1.0) 



= +8.25 



EIA 30 = 



- ( M ^- L£ ) (-1.0) + (-5x6)(-1.0) + 



3x9 



(-1.0) 



= 45.75 

EI&u = I ^T^" 1 (2.0) + (3.0)(5.0)(3.0) 



= 49.5 

Elb l2 =( l - A ^ 1 ) (-3.0) + (1.5)(-0.5) 



= -14.25 = £/5 21 
^/6 13 = I^T^-I (-1.0) + (1.5x5)(-1.0) 



= -17.25 =£Y6 31 



EI8 22 = ( - 



3 x3 \ , ^ / 15 x3 \ , ™ I 3x3 . , 
— J (-2.0)+ I-— J— J (-3.0)+ (^- — J (-2.0) 



2x2\ 4 



= 34.17 




2x2 



)(-».0)+(-^)(-1.0)+(-^)(-l. 

) (-i.o) 



0) 



= 11.5=£76 32 

EI8 33 = (- 1 .0 x 3)(- 1 .0) + (- 1 .0 x 3)(- 1 .0) + (- 1 .0 x 5)(- 1 .0) 
= 9.5 

Substituting these values into the elastic equations: 



[-136.125] 




1" 49.5 


-14.25 


-17.25 


8.25 


+ 




34.17 


11.5 






-17.25 


11.5 


9.5 





\ Xl ] 




0] 




X 2 




0 




*3. 




.0 



30 



METHODS OF CONSISTENT DISPLACEMENTS 




S 




4kN 



CVJ 



27 



E 

o 



B 



E 
O 

to 



3kN 



1 



4kN 



9. 



(a) Actual structure 



(b) Base structure 




(c) M Q - diagram (d) M 0 /EI - diagram 



3.0 




{+) 



1.5/EI 




X, = l 



3.0 



3 /EI 



(e) - diagram 



(f) m,/£7 - diagram 



Figure 2.10 
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(i) m 3 - diagram (i) m 3 /EI - diagram 

Figure 2.10 contd. 

Solving the simultaneous equation: 

V D =X, =2.34kN (upward) 
Hn = X-i - 1 .56 kN (to left) 

Md ~%3 = 2.45 kN m (counter-clockwise) 



2.5 THE ELASTIC CENTRE METHOD 

The elastic centre method is a special method of solving statically indeterminate 
structures of the one-loop form. Rigid-jointed portal frames, single-bay gabled 
bents, single-span arch, closed or ring structures are examples of the type of 
problems easily solved by this method. 

Consider a fixed arch as in Fig. 2.1 1(a) under an arbitrary loading which 
produces bending moments M 0 in the primary structure. 
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X 2 =H A 



(a) Real structure 



B 




B 



x,= v A 



(b) Base structure 



Figure 2.11 



The primary structure is taken to be a cantilever as in Fig. 2.1 1(b), where the 
left support is removed and the redundant reactions X x = V A ,X 2 =H A and 
X 3 = M A are applied at the support point. The compatibility equations are 



A.ol 




[5ii 


5 21 




+ 


S21 


5 22 


A30 




.«3i 


6 23 



531 

5 23 
633 

The displacements are evaluated as 

















0 








0 



10 



20 



_ { M 0 m x ds _ f M 0 xds 
"J EI ~ J EI 



M 0 m 2 ds _ f M Q yds 
EI "J EI 
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[M 0 m 3 ds 
EI 



_ C M Q ds 
J EI 



[2.20] 



[2.21] 



since m l =x, m 2 -y and m 3 = 1. 
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Also 



_ f m\ds f 
11 J EI "J 



x 2 ds 
EI 




m\ds f y ds 




EI J £/ 



if 



5=5 - [ ^1^2^ _ [ sygfr 

12 21 J £Y J £7 



bl2=831= )~EI~ = )ET 



5 23 " 




m 2 m 3 ds _ ?yds 
EI ~)Jl 

If ds/EI is considered as an elemental area of length ds and a width normal to 
the arch axis of l/EI, then the following interpretations may be made: 

5 [ j and 5 22 - moment of inertia of I /EI area about the j> axis and jc axis, 

respectively 

5 33 = total l/EI area of the arch 

5 i2 = product of inertia of \/EI about the given axis 

6 23 = statical moment of \/EI area about thej> axis and* axis, 
respectively 

If the origin of the axes can be transferred to the centroid or elastic centre 
of the elastic area, ds/EI, the computation may be simplified by the virtue of the 
fact that 6 13 and 5 2 3, being the statical moments of elastic areas, disappear. Also 
if the axes through the elastic centre are the principal axes, b n > being the 
product of inertia, also vanishes. 

It is statically possible to transfer the forces X x - F A and X 2 - #a to any 

point, provided X 3 = M A is properly modified, since any force may be replaced 

by an equal parailel force acting through any arbitrarily chosen point and a 

couple. Accordingly, the redundants may be applied at point 0 (x 0 ,y 0 ) which is 

attached to A by a perfectly rigid arm (Fig. 2.12). It is clear that this arm does 

not fundamentally change the structure, since, being rigid, it makes no direct 

contribution to the deflection of the arch. 
Taking 0 as the origin of coordinates, 



M A * X 3 + X&o + X&o 



[2.23] 
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B 



Xo 



Rigid arm 



Figure 2.12 

If the structure is symmetrical, the centroidal-principal axis will coincide with 
the axis of symmetry, and hence the product of inertia will be zero. Thus 



21 



13 



31 



23 



32 



The corresponding compatibility equations are 
Aio+*i5n =0 

A 2 o + ^2^22 ~ ^ 

A30 +^3633 = 0 



[2.24] 



where the redundant reactions are taken at the elastic centre. 

If the structure is not symmetrical where the* and y axes are not the principal 
axes, 6 12 will not vanish, while 5 13 and 5 23 disappear. The compatibility 
equations for this case are 



A 10 *X\b\\ + -Y 2 5 2I =0 
A 2 o + X\ 5 12 + -^2^22 = 0 

a 30 +^3633 = 0 

Solving the simultaneous equations 



[2.25] 



X-7 = - 



&10 -(§12/522)^20 

5 11 — ( 5 12/ 5 22> 5 12 

^20 - ( S 12/5n)A 10 
5 22 - Qllfillfiu 



[2.26] 



X Z - ~ A 3o/ 5 33 



EXAMPLE 2.4 Determine the reaction components at D of the rigid frame of 
Fig. 2.13, using the elastic centre method 

The elastic centre is located by taking moments about AB for the x coordinate 
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(a) Real structure (b) Base structure 

-3 




(e) m 2 - diagram (f) m 3 _ diagram 

Figure 2.13 
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and about BC for they coordinate: 

(3)(3) | (3)(1.5) 

1 2 
*o = 3 = 1.18 m 

5+- + 3 

(5) (2.5) | (3) (3) 

1 2 
y Q = — = 1-79 m 

The displacement coefficients are computed using graphic multiplication method: 

£7A 10 = ^^y^j(-0.68) + (-6x5)(-1.18)+ ("^) (-1.18) 
= 52.87 

*/A,. = (-^) (.,9) + (-6x5) (=23^12) ♦ (-^)(-2.2 
= 47.11 

EIA 30 = (-^y^) (-1.0) + (-6x5)(-l)+ ^-^ij(-i.O) 
= 47.75 

£»„ =(-1.18x5)(-1.18) + l (-fx Lis) 

4(-^)(-fxl.l8) + (3xl.82)(1.82) 
= 18.18 

/ 3.21 x 3.21 \ / 2 , \ /l.79 x 1.79\ /2 

j— ) (-3x3.21)^—^—) 

4(1.79.3X1.79)* (-1^2) (-fxl 
/ 1.21 x 1.2 \ /2 \ 

= 20.25 

EI8 33 = (5 x 1)(1) + 1 (3 x 1)(1) + (3 x 1)(1) 
= 9.5 
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Since the axes through 0 are not principal axes, 5 j 2 will not vanish. Thus, 

^-.fi311ia K , ( ^ + Iflfi=L!fi ll , aT9| 

t M9- L 2 lx3(1 _ 82) 



= 6.63 



The elastic equations are 



r A, 0 j 






+ 


A 30 





6 n 5,2 

5 12 5 22 

0 0 



0 
0 



33 









0 








0 








0 



substituting 



f 52.87] 




47.11 


+ 


47.75 





18.18 6.63 



0 



6.63 20.25 0 



0 



0 



9.5 









°1 




x 2 




0 


. 






0 



The solution of the simultaneous equations is 

X x = -2.34 kN 
X 2 =-1.56 kN 
X 3 = -4.82 kN in 
The values of the redundants at support D are 
V D = 2.34 kN (upward) 
H D = 1.56 kN(left) 
M D =4.82 - 2.34 x 1.82+ 1.56 x 1.21 



= 2.45 kN m (counter-clockwise) 



2.6 THE THREE-MOMENT EQUATIONS 

Consider a continuous beam with n spans as shown in Fig. 2.14. This beam is 
indeterminate to the (w - l)th degree when the support reactions are taken as 
the redundants, each of which contain all the unknowns. However, when 
support moments are used as the redundants, although the same number of 
equations must eventually be solved, each equation contains only three of the 
unknowns. The latter choice of redundants localises the loading conditions on 
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vfy///. ' V^SS* '/fy///. '.'/$///. >/fy/M "J}}///,. 

2 I J K n n+l 



h-H 11 -t- " 



Figure 2.14 

the base structure and the resulting relationship between the redundants permits 
the equations to be written in a simple and systematic manner. These equations 
express the three-moment equations first presented by the French engineer 
Clapeyron. 

The three-moment equation expresses the relation between the bending 
moments at the three successive supports of a continuous beam. The relation is 
derived from the continuity of the elastic curve when the compatibility equations 
are obtained in terms of the support slopes of adjacent spans. 

Consider two adjacent spans IJ and JK of a continuous beam shown in 
Fig. 2.15. The moment of inertia is considered constant between I and J and 
equal to /, and likewise constant between J and K and equal to/,-. The beam is 
assumed to be initially straight, and the support settlements amounting to A,-, Ay 
and A k take place at support I, J and K respectively, as indicated by the heavy 
line in Fig. 2.15(e). 

Compatibility equations are written at each support expressing the equality 
of end slopes at adjacent spans. The condition of continuity of the slope gives 

A/ - A,- + Sf _ A k -Aj-8 k 
Li " Lj 
Rearrange the equation as 

5; 6 fc A; - A; A k — A; 

-± + -A = -L-— >+-*-—> [2.27] 

Lj Lj L t Li 

But from the Second Theorem of the Area-moment Method, 

5 ' = ^-[^ + (|^) (V3) + (|jW£/) (2Z.//3)] 



5„ = 



=~ ^A t Xf + (±M k Lj j (Lj/3) + {^MjL^ (2L,/3)j [2.28] 
Combining [2.27] and [2.28] gives Clapeyron's Equation of Three Moments: 
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(c) Base structure bending moment 




\- * 4. f -I 

(e) Beam after displacement 

Figure 2.15 
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In the particular case when /, =/,•=/ and there are no support settlements, the 
equation becomes 

MjLj + 2Mj{Li + Lj) + M k Lj = -6 [2.30] 

Li 

In the special case when the continuous beam has a fixed support, the bending 
moment at the fixed support is a redundant for which a new equation must be 
established. This case may be treated by considering the fixed support equivalent 
to an outer imaginary beam of finite span with infinite stiffness. This is illustrated 
in Fig. 2.16. 



P 2 [P* W 





1 1 1 I I I I 1 I I 1 IXUJ (a) 




I — -h — H — I 



X * 




a, P I"' t 





w 

I I I II 1 II I I 1TTTT] (b) 




Li I L 2 , L 3 





1 

Figure 2.16 

The three-moment equation, when applied to span A 0 A and AB is 



0+2Ma ^ +Mb ^ = _ 2 



6Ax 
LI 



This reduces to 



6Ax 

2M A L+M B L = -I, 



EXAMPLE 2.5 Find the support moments of the continuous beam shown in 
Fig. 2. 17 using the three-moment equations. 

The three-moment equations are written for: 
Span A 0 Aand AB 

„ (*<o\ (Lo 4 \ », / 4 \ 6x12x7/3 
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4 m 
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8kN 
2kN/m 

nmn 
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8m 
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4kN 



Figure 2.17 




Span AB and BC 
Span BC and CD 



+ M, 




[(6) (12) x 5] 3 _ (6) (149.33)4 
(2) (4) (8) (8) 



+ M 




(6) (149.33)4 (6) (18)3 



(8) (8) 



(3) (6) 



Simplifying: 



2A/ A + M B = -10.5 
2M A + 6M B + M c = -7 1 



M B +6M c = -58 



Solving: 



M A =+0.01 kN m 
M B = -10.52 kN m 
M c = -7.91 kN m 



2.7 THE METHOD OF ELASTIC WORK 

Consider the indeterminate beam shown in Fig. 2.18 which has two redundant 
reactions. 

Taking the base structure as a simply supported beam, as shown in Fig. 2.18(b), 
subjected to the loads and the redundant reactions M A and V B , it is obvious that 
the base structure is statically equivalent to the actual structure. Compatibility 
condition furnishes two additional equations, namely, the rotation at A and 
zero. These displacements by Castigliano's First Theorem are 

0a= — = 0 



A B = 



dU 



3F, 



= 0 



[2.31] 



B 
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(a) 





M t I I I I I I 



w 




(b) 

Figure 2.18 



Therefore, the condition for determining the redundantsiW A and K B satisfying 
the displacement conditions is such that the total internal energy will be a 
minimum, or strictly speaking a maximum. It is evident that it must be a 
minimum value, because there is no maximum value when the stored energy 
increases, as the redundant forces increase indefinitely. The above equations may 
also be interpreted as follows: Among all possible sets of values that redundant 
forces in the system may assume, the correct set of values is that which makes 
the strain energy U a minimum. Since one equation is obtained for each 
redundant force, a set of equations corresponding to the compatibility will 
result. 

Therefore, Castigliano's Theorem of Least Work, also commonly known as 
Castigliano's Second Theorem, may be stated as follows: 

In any loaded statically indeterminate structure with no temperature changes 
or support displacement, among all possible equilibrium states the correct 
values of the redundants are those for which the total elastic strain energy, 
resulting from the application of a given system of loads, has a minimum 
value. 

Consider the beam shown in Fig. 2.18. The bending moment at any point on 
the beam is given by 



M = M 0 + M A w A + V B m B 



[2.32] 



where M a 



m A 



moment in the base structure due to the applied loads 
moment resulting from a unit couple at A 
moment resulting from a unit vertical load at B 
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Using Casligliano's Theorem of Least Work 

dU .Am dx 



e,- — -M\ 



M A J dM A EI 

f . , dx f 0 dx f dx ^ 

= ) M a m A — + M A \ mi — + R B | m A m B — = 0 

dx . . f dx __ f 2 t/x _ Q 



f etc f d* „ f - 

= \ M 0 m B — +M A \ m B m A —+V B \m B 



EI 



From the above two equations, the integrals are evaluated for the base 
structure and the redundants M A and/? B are tnen solved. 
Writing the above equations in terms of displacements: 

6 A =A a +M A 8 aa + V B 8 ab = 0 

A B =A b +M A 8 ba + V B 8 bb = 0 [2.34] 

Note that these equations are identical in form with those obtained from the 
virtual work method. 

EXAMPLE 2.6 Find the reactions of the propped cantilever beam shown in 
Fig. 2.19 using Castigliano's Theorem of Least Work. 



A 



L/2 1 L/2 





(a) Actual structure 



PL -V D L 



M vI____I 




V A = P-V B 

(b) Base structure (U B as redundant) 





V B 2 L 



(c) Base structure (M A as redundant) 

Figure 2.19 
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(a) Taking the reaction at B, K B , as the redundant (Fig. 2.19(b)), the redundant 
F B is determined from 



3F B 



Table 2.2 



Origin 

Section * = 0 at Limits M dM/bV B 




AC C x = 0 to V B | x + - 1 - Px | x + 

x = L/2 

CB B x = 0 to T/ 

x = 1/2 Fb * 



3K B J 9 Kg £/ 




o 



1 |~1 / L\ 3 Px 3 PLx 2 ] 



L/2 



EI 




1 ["K B I 3 _5PL 
L 3 



£7 3 48 



= 0 



Thus 



The other reactions are 

R A =P-V B =^P 

PL 3 
Af A = — - V B L = — PL (negative moment) 

2 16 
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(b) Taking M A as the redundant (Fig. 2.19(c)): 



Table 2.3 



Origin 

Section x = 0 at Limits M dM/bM A 

m a\ , * 

AC A x = 0to -M A + - + — * -1+7 




x = L/2 \ 2 1 I 1 



M A \ x 
CB B x = 0to |-- — \x 



x = L/2 




L 



9A/a EIU\\ Ma + \2 + T U 



■4TJ(*-¥-f-*»(^)?)> 



+ — 

EI 




£7J 0 \ 2 L J L 
3 ~ 16 




= 0 



3/>Z, 

M A = — — (negative moment) 
16 

The other reactions are 

2 L 16 
B 2 Z, 16 

Note that a moment that causes tension on the top side of the beam is 
considered negative moment. The positive value obtained for M A indicates that 
the assumed negative moment is correct. 

EXAMPLE 2.7 Find the maximum moment and the decrease in the vertical 
diameter of the ring shown in Fig. 2. 20. 

Due to symmetry, one-quarter of the ring is considered. 
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(a) Actual structure 




(a) Maximum Moment 

W n 

Oa = — = 0 
A *M 0 

bU _ M P dM 
dM 0 ~El)M 0 

PR 

M= y (1 -cosd)-M 0 
dM = l 



9M 0 
dx = Rdd 
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Substituting, 



3M Q 



1 f*/ 2 (PR , , \ 

= -Ti) 0 ( T 0-«»*>-J<,)jw«o 




PR 2 . PR 2 . . . |»/2 

o 



2 0 + — sin 0 4 #M O 0 




PR I 2 

The moment at any point is 

J?R 

M = — (1 -cosfl)-^ 
= — | - - cos 6 

The maximum moment is at the point of load application, i.e. when 6 = n/2 

M ~ PR 

7T 

(b) Deflection 

B dP J EI \ BP 

„ p * I 2 
M = — --cos0 

2 Vtt 

m r 12 

— = — I — - cos 0 

ap 2 \7r 

etc = .fa/0 
Substituting: 

4 ["I 2 PR 2 (2 V 

4 4 11 
0 - - sin 0 + - 0 + - sin 2d 
it 2 4 





Castigliano's Theorem of Least Work provides a more suitable method for 
the analysis of partially or completely articulated structures. To illustrate the 
method, consider the three-span continuous truss (Fig. 2.21). 
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1 L ' I 



L 2 



(a) 




t-3 




(b) 

Figure 2.21 



By removing the interior supports and replacing them by the forces R B and 
R c which are equal in magnitude and direction to the respective reactions, the 
structure of Fig. 2.21(b) becomes statically equivalent to the actual structure. 
The bar force in any member of the structure is 



N=N 0 + R B n B +R c n c 



[2.35] 



where 



N Q = bar force in any member of the base structures due to applied loads 
n B = bar force in any member due to a unit vertical force applied at B 
tiQ = bar force in any member due to a unit vertical force applied at C 

The expression for the total internal (strain) energy is 



N 2 L 



[2.36] 



The equations for determining the redundants are written from Castigliano's 
Theorem of Least Work: 



bU NL dN 



dR B EA dR B 



= 0 



[2.37] 



3(7 NL dN 



dR c EA dR c 



= 0 
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Note that 



W _ 
3*b = " B 



= n c 



Substituting [2.37] into [2.35] gives 



EA EA EA 

or, in terms of displacements, 

A b +R u 8 bb +R c 8 bc = 0 
A c +R B 8 cb +R C 8 CC = 0 



[2.38] 



[2.39] 



These equations are identical in form with the equations obtained from the 
virtual work principle. 



EXAMPLE 2.8 Find the reaction at B and the bar force in member BF using 
Castigliano 's Theorem of Least Work (Fig. 2. 22). The cross-sectional area of the 
members, in cm 1 , are shown in parentheses. 



The truss is indeterminate to the second degree with one redundant member 
(internal indeterminacy) and one redundant reaction (external indeterminacy). 
Taking member BF and the reaction at B as redundant, the two condition 
equations of least work are 



dR B dR B A 

dU bN L n 

= 27V— — - = 0 



3F BF dF g p A 

All computations necessary for the evaluation of the terms of the above equations 
are shown in Table 2.4. 

The required equations of least work are: 

135.09/? B + 54.53F BF - 1680.8 = 0 
54.52K B + 1 18.2F BF - 202.3 = 0 
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(c) Due to external redundant R B 



0.0 



-0.8F BF 




(d) Due to internal redundant F BE 

Figure 2.22 
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The solution of the simultaneous equations is 
R B = 14.44 kN (upward) 

F BF = -4.95 kN (compression) 



2.8 PROBLEMS 



2.1 Evaluate the moments at supports A and B of the beam shown in Fig. P2.1 



(Ans: M A = -0.0626L 2 , M B = -0.0036 \L 2 ) 

IkN/m 

A '&U11 1 1 1 1 1 I I 1 1 I 1 I I 



B 



21 
L/2 




I 

U3_ 



Figure P2.1 

2.2 Determine the reactions and bar forces for the truss shown in Fig. P2.2. 



5 KN 




Figure P2.2 

2.3 Determine the bar forces in the truss double diagonal system shown in 
Fig. P2.3. The area of all top chord members is 2A and the area of the 
remaining members is A . 

(Ans: N DF = -3.05 kN N BE = -5.56 kN A^ BD = -3.65 kN 
N EF = 3.00 kN JV CF = -2.80kN yV BF =-3.76kN 
W A d = -2-76 kN 7V AF = 4.60 kN N CE = 4.66 kN) 



lOkN 




10 kN 



Figure P2.3 
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2.4 Construct the shear and bending moment for the continuous beam shown in 
Fig. P2.4. 



2 kN /m 

MI1M11IIIL1II MINIUM 




B 



5m 



4m 




ZkN 



Figure P2.4 



2.5 Draw the bending moment of the frame shown in Fig. P2.5. 

(Arts: M A = +5.40 kN m M D = -2.65 kN m 

M B = -6.75 kN m M E = +2.71 and -7.29 kN m 



M c = +7.09 kN m) 



5kN 



B 



21 



IOkN mO 



31 



'///////. 



E 



T 



E 



wwi F 



4m 



6 m 



Figure P2.5 

2.6 Find the support moments of the continuous beam in Fig. P2.4 using the 
method of three moments. 

2.7 Find the joint moments of the frame shown in Fig. P2.6 using the elastic 
centre method. 



04ns: M B =-M D =0.42kNm) 



B 



21 



10 kN m 



7Z> 77a E 



5" ^ 1 



Figure P2.6 
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3. The Slope Deflection 
Method 



3.1 INTRODUCTION 

The slope deflection method was presented by G. A. Maney in 1915 as a general 
method of analysis for continuous beams and rigid-jointed frames. One of the 
main values of the slope deflection method lies in that it forms the basis of other 
methods such as the Cross and Kani moment distribution methods. These 
methods are numerical iteration solution of the slope deflection simultaneous 
equations. 

In the slope deflection method the rotations and relative joint displacements 
are the unknowns and the moments at the joints are found in terms of rotations, 
joint displacements, stiffness and length of members which are determined from 
the solution of simultaneous algebraic equations. The method assumes all joints 
to be rigid; that is, the angles between the members at the joints do not change 
under applied loads. Also, displacements due to axial and shearing forces, being 
very small, are neglected. Consequently, the joint is considered to rotate as a 
whole and a single angle of rotation is sufficient to define the rotation of all 
members at a joint. By properly determining the joint rotation and displacements 
at the ends of the individual members, the joint moments are formulated from 
the equations of equilibrium. 

3.2 DEVELOPMENT OF SLOPE DEFLECTION EQUATIONS 
3.2.1 Sign Convention 

There are different sign conventions used in different textbooks. In all 
subsequent discussion a statical sign convention suitable for beam and frame 
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analysis is adopted, as follows. 

(a) Moments An end moment M is considered positive if it tends to rotate the 
member counter-clockwise or the joint clockwise . That is, if the member tends to 
rotate the support clockwise, the moment is considered clockwise. 

(b) Joint Rotation The rotation 0 of a joint is positive if the tangent turns in 
a clockwise direction. 

(c) Member Rotation The rotation of the chord connecting the ends (A/Z-) is 
positive if the member turns in clockwise direction. 

3.2.2 The Slope Deflection Equations 

To develop the slope deflection equations, consider a base structure which is a 
straight prismatic restrained beam AB of span L and moment of inertia/ 
subjected to a system of loads. Due to the effect of loads and the adjacent 
members, the beam is displaced as shown in Fig. 3.1(b). With reference to the 
original position, the ends rotate through 0 A and 0 B at A and B, respectively, and 
a relative vertical displacement of A of the member ends. The end moments 
produced at A and B are, respectively M AB and Af BA . The member AB is 
equivalent, both statically and kinematically, to the sum of those shown in 
Fig. 3.1(d)— (f)- The effect of each displacement component to determine the 
end moments at A and B is studied separately. Note that the end and member 
rotations are assumed clockwise and therefore positive. 
These components are : 

(a) Fixed-end moments 

For the applied load system the fixed-end moments M F AB and M$ A are first 
determined. 



(b) End moments due to rotation 0 A 

The member AB is propped cantilever with clockwise rotation 0 A (Fig. 3.1(d)). 

Using the conjugate beam method, the rotation and displacement at support A 
are 



(M AB +M BA )L 



2EI 
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H L - -| 

(a) Original position 




(c) Fixed-end moments 




(f) Relative displacement 

Figure 3.1 
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From which 



Therefore 



_ M AB L 



Hence 



Mab - 4E ™ A and M BA - - 



(c) End moments due to rotation 0 B ( Fi g- 3.1(e)) 

In a similar manner as above, 

_ « B +A/ BA )Z. 

0B "~ Wl 



and 



_ _M* K L ( L\ M'^L IlL 
B 2£7 \ 3 ) 2EI \ 3 




From which 



^ab - -^"^ba 



Therefore 



a _ ^baZ- 



Hence 



_ 4£Y0 B _2£Y0 B 
M BA — — - — and /W AB - — 



f<y moments due to relative displacement A (Fig. 3.1(0) 

The slope at support A is 

_ (M A ' B +M BA )L 
A 2EI 
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Thus 



Af AB = -M B A 



Using the conjugate beam method 



Mr A / L \ Ma'r I 2L 



. _ mba ( _ \ + ab 
2£Y \3/ 2£Y 

_ M BA L 2 
6£Y 



Hence 



<b=--^- and < A =+-^- [3.3] 

The separate moments of Fig. 3.1(c)— (f) may be superimposed to represent 
the true end moments for beam AB according to the bending moment sign 
convention: 

M BA =< A +M BA +M BA +M B ' A 
or in terms of the end rotations and displacements: 

IF I 

M AB =< B + — (20 A + 0 B - 3A/L) 

1EI 

M BA = Mg A - — (20 B + d A - 3A/L) 

Adopting the statical sign convention and introducing K = I/L the equations are 
revised as 

M AB = M£ B - 2EK(26 A + 0 B - 3 A/Z,) 

[3.4] 

^ba = ^ba - 2£ff(20B + 0a - 3 A/Z,) 

The above equations are known as the slope deflection equations. They 
express the end moments of a member in terms of its end displacement. Note 
that the sign of the fixed-end moment is shown as plus. However, its correct sign 
will be determined by the direction it tends to rotate the joint following the 
adopted sign convention. 

These equations can further be generalised as a single equation in the form 

M jm = Mfrn - 2^(20,- + 6 m - 3A Jm /L fm ) [3.5] 

where the subscripts / and m represent the near and far end of member JM 
respectively. 

The slope deflection equations given by [3.5] represent two equations with 
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six unknowns. The unknowns are solved by applying boundary conditions for 
the displacements and continuity conditions for the moments. 



3.2.3 Fixed-End Moments 

When a beam whose supports are completely fixed against rotation or translation 
is subjected to transverse loads the beam is called a fixed-end beam. The bending 
moments at the supports of such a beam are called fixed-end moments. The 
values of the fixed-end moments for common types of loadings are given in 
Table 3.1. 



3.3 APPLICATION OF SLOPE DEFLECTION EQUATIONS TO BEAM 
PROBLEMS 



3.3.1 Beams With No Support Settlements 

The slope deflection equations for a member with no relative lateral displacement 
between the ends are: 

M AB = M£ B - 2EK{2Q A + 0 B ) [3 6] 

M BA =M B 7 A -2EK(2d B +6 A ) 

Consider the two-span beam shown in Fig. 3.2. 



Fixed-end Moments 

The fixed-end moments depend on the applied loading. The fixed-end moments 
are M AB and M BA for span AB and M BC and Mqb for span BC. 



Slope Deflection Equations 

The slope deflection equations are written for each member as follows: 
M AB =M A T B +2EK AB (2d A +d B ) 
M BA =Mg A +2EK AB (d A +2d B ) 
Mbc =Mg c + 2EK BC (2d B + d c ) 
M C b =M£ B + 2EK BC (d B + 29 c ) 
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Table 3. 1 



A 



w 



w 



M 



B 




B 





B 







/to 2 


ft 2 * 


A 2 


L 2 




wi 2 


12 


12 


u* 2 


wZ. 2 


30 


20 


h-L 2 / O o fl 2 \ 




60i V /J 


-£(-»f«£) 


5wL 2 


5wL 2 


96 


96 


-*<2a-ft)p 


- fl (2fc - a) ^ 



Equilibrium Condition of Joints 

The equilibrium equations are written for each joint by taking the free-body 
diagram of the joints. These equations are 

2M A =M AB =0 

2M B = M BA + M BC = 0 [3.8] 
2A/ C =M CB =0 
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Figure 3.2 

The equilibrium equations [3.8] may be written in terms of the rotations 
given in [3.7]. Hence 

M F B + 4EK AB 6 A + 2EK AB 6 B = 0 

Mba +m£ c + 2EK AB 6 A + (4EK AB + 4EK BC )d B + 2EK BC 6 C = 0 [3.9] 

Mq B + 2EK BC 6 B + 4EK BC d c = 0 
The above equations may be written in matrix notation as follows: 









2EK AB 0 1 








[0 


^BA + ^BC 


+ 


2EK AB 


4£(a: ab +a:bc) 2EK BC 




0B 




0 


Mq B 




0 


2EK&Q 2EK$c 




0r 




0 



[3.10] 



The above equation may be written in compact form as 

\M F \ + [K){d\= \0\ [3.11] 

where I M F \ is the vector of the fixed-end moments, [A'] is the member stiffness 
matrix, and \ d] is the vector of end rotations. In [3.10] the unknown quantities 
are the end rotations which can be determined from the solution of the 
simultaneous equations. These rotation quantities can then be substituted into 
[3.7] to find the support moments. 

The procedure for analysing continuous beams is as follows: 

(a) Determine the fixed-end moments in each span. 

(b) Use the slope deflection equations to express the end-moments and end- 
rotations. 

(c) Establish the equilibrium equations of moments where the support 
rotations are the unknowns, at each joint capable of rotating: the sum of the 
end moments of all members at the joint is zero. 
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(d) Evaluate the rotations by solving the simultaneous equations. 

(e) Substitute the rotations back into the slope deflection equations to compute 
the end moments. 

EXAMPLE 3 . 1 Determine the support moments for the beam shown in 
Fig. 3.3. 



100 kN 

r30 k N/m 
MINI II 1111 Mill 

6 " ' 




Figure 3.3 



Relative Stiffness 



(9 ■ (5) ■ 



(*ab):(* BC )= TT : 77 =5:8 






6m | 4m | 


| m 8m m j 





Fixed-end Moments 



30 x 8^ 

m£b = -M£ A =—^—= 160 kN m 

__p 100x6 x4 2 _ ivt 
A*bc= J7j2 = 96 kNm 

,~ 100x4x6 2 

^CB = - ^2 =-144 kNm 



Sfope Deflection Equations 

M AB = 16O-5(20 A +0 B ) 

= 160 - 50 B (since 0 A = 0) 

M BA =-16O-5(20 B ) 
M BC = 96 - 8(20 B + 0 C ) 



A^ B = _144-8(0 B + 20 C ) 
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Equilibrium Equations 

AtjointB: M BA +M BC =0 

At joint C: Mq B = 0 
Hence 

260 B + 80 c = -64 
80 B + 160 c = -144 
In matrix form 



End Moments 

Substitution of 0 values into the slope deflection equations yields the support 
moments according to the statical sign convention: 

M AB = 160 - 5(0.364) = 158.18 kN m 
M BA = -160 - 5(0.726) = -163.64 kN m 
M BC = 96 - 8(0.728 - 9.182) = +163.64 kN m 
M CB = -144 - 8(0.364 - 18.354) = 0.0 kN m 

3.3.2 Beams With Support Settlements 

Support yielding in continuous beams induces bending moments at all supports 
of the beam. The general slope deflection equations are used to analyse 
continuous beams. The equations are 



L 8 16J L0 c J L-144J 
Solution of the above equation yields 






[3.4] 
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EXAMPLE 3 .2 Determine the support momen ts of the con tinuous beam 
shown in Fig. 3.4. The support at A rotates through 0. 15 radian in a clockwise 
direction and the support at C settles down 10 mm; E = 2.1 \ 10 s N/mm 2 and 
1 = 4x10* mm 4 

El = 2.1 x 10 s x 4000 x 10" 6 = 8400 kN/m 2 



^ Iq a = 0.15 rad Settlement 




A = / Omm 



" - ~ Ac 1 



e 

6m , 4m 




Figure 3.4 



Slope Deflection Equations 



2 x 8400 x 2 

M AB = - (20 A +6 B ), where 0 A = + 0.15 radian 



= -1680 - 56OO0 B 

2 x8400x2^ 
M BA = (20 B + 0.15) 



= -112O0 B -840 
2 x 8400 x 3 

M BC = (20 B + 0 C - 3 A/6) (where A = 10 mm 



= -168000b - 84OO0 C + 42 



= 0.01 m) 



= -168O0 C -84OO0 B +42 
Mcd = - 2 X ! 4Q ° 20 c +0 D - |- 3X ' " 1 I (where0 D =0) 



= -84000c -31.5 



_ 2 x 8400 

m dc — : 



(«c + — ) 



= -42000c - 3 1 .5 
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Joint Equilibrium Conditions 

(i) Joint B : M BA + Mqq = 0 

(ii) Joint C :Mc B +M CD =0 

Substituting the moment expressions into the equilibrium equations gives 



-28000 -840C 
8400 -25200 




7981 
1050j 



Solving the simultaneous equations, 
0 B = -0.03181 radian 
0 C = +0.01101 radian 

End Moments 

M AB = -1680 - 5600 (-0.03181) = -1501.9 kN m 

M BA = -1 1200 (-0.03181) - 840 = -483.8 kN m 

M BC = -16800 (-0.03 181) - 8400 (0.01 101) + 42 = 483.9 kN m 

Mq B = -16800 (0.01 101) - 8400 (-0.03181) + 42 = +124.2 kN m 

M CD = -8400 (0.0 1 1 0 1 ) -3 1 .5 = - 1 24.0 kN m 

M DC = -4200 (0.01 101) -31.5 = -77.7 kN m 

Converting the statical sign convention into bending moment convention 
yields 

M A = +1501.9 kN m (tension at bottom fibre) 
M B = -483.9 kN m 
Mb =+124.1 kNm 

M n = 77.7 kN m 



3.4 APPLICATION OF SLOPE DEFLECTION EQUATIONS TO FRAMES 



3.4.1 Frames Without SMesway 

The slope deflection equations for frames without sidesway are: 
M AB =M A 7 B -2EK(2d A +d B ) 
M BA =Mg A -2EK(2d B +d A ) 
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Under this case all joints of the frame remain fixed in location during loading. 
Such frames are either physically held against sidesway or are symmetric and 
subject to symmetric loading system. 

In the slope deflection equations the joint rotations are considered as the 

unknown while the joint moments are determined from joint conditions of 
equilibrium. In effect, the slope deflection solution of a frame without 
sidesway is essentially the same as that of a continuous beam. For frames 
without sidesway, there are always as many conditions of equilibrium of joints 
as unknown rotations. After the rotations are determined, the joint moments 
can be found from the slope deflection equations. 

3.4.2 Frames With Sidesway 

When loads are applied to frames, there are cases in which lateral movement 
occur through unknown distances, although usually in known direction. The 
effect of such translation is to cause the joints to deflect relative to initial 
unloaded positions. In such cases, the joint rotations and relative displacements 
are the unknowns in the slope deflection equations. The governing slope 
deflection equations for members subjected to a sidesway are 



Assuming that all axial deformations are so small as to cause no change in 
member lengths, the relative sidesway of the joints may be evaluated depending 

on the geometry of the frame. 

Figure 3.5 shows a two-column frame subjected to a lateral force P applied at 
B. The frame will deflect to the right while joints B and C rotate clockwise. Thus, 
the frame has three unknown displacements namely 6%, 9 C and A. 



M AB =M£ B - 2EK{ld A + 9 B - 3A/L) 
M BA = M BA - 2EKQB* + 0 A - 3A/L) 



[3.4] 




(a) 



(b) 



Figure 3.5 
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With three unknown displacements, three condition equations based on static 
equilibrium are required. Two equilibrium equations can be written for the two 
joints B and C. The third independent equation can be formulated by considering 
the horizontal static equilibrium of the frame. Taking the columns as free bodies, 
as in Fig. 3.5(b), 



#ab ~ 



^AB + M BA 



[3.12] 



_ M CD +M DC 

«bc 1 

Applying the equilibrium condition ZH = 0 for the whole frame, 

P-H AB -H DC =0 [3.13] 

This additional equation is generally known as the sway equation and sometimes 
known as the shear equation or the bent equation. 

EXAMPLE 3 .3 Find the joint moments using the slope deflection equations 
of the frame shown in Fig. 3. 6 

Relative Stiffness 



Member 


AB 


BC 


CD 


Moment of inertia 


/ 


21 


/ 


Span (m) 


5 


3 


3 


I/L 


ns 


27/3 


//3 


Relative K 


3 


10 


5 



4kN 



o 




3kN 



A '////////, 



3.0 



Figure 3.6 
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Fixed-end Moments 



<B= 2 g^ = +1.44kNm 
_ 3(2) (3) 2 _ 

M BA " Y$ " 

«Sc = -A^ F B=-^r-=1.5 kNm 

S/ope Deflection Equations 

M AB = 1.44-6(0,, -0.6A) 
Af BA = -2.16-6(20 B -0.6A) 
M BC = 1.5 -2O(20 B + 0c) 
M CB =-1.5-2O(0 B +20 C ) 

A4bD = -lO(20 c - A) 
M DC = -1O(0 C - A) 

/owir Equilibrium Equations 

(i) At joint B : M BA + M BC = 0 

(ii) At joint C : M CB + M CD = 0 



SAiar Condition 



Substituting 

„ _ M AB +M BA +6 
//a- 5 



and 



Mqd +M dc 



into the shear condition 3.0 - H A - H D = 0: 





M 



DC 




45 - 18 - 3(M AB +M BA ) - 5(M CD + M DC ) = 0 



540 B + 15O0 C - 121.6A = -29.16 
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or 

-3.60 B - 10.00c + 8.1 1A = 1.944 

Substituting the moment expressions into the joint equilibrium equations gives 
the following simultaneous equations: 



56 


20 


-3.6 1 








-0.6671 


20 


60 


-10 




Be 




-1.500 


-3.6 


-10 


8.11 

i 




A . 




[ 1.944 



Solving the three simultaneous equations: 

0 B = -0.00195 
0 C =0.0195 
A = 0.263 

End moments 

M AB = 1-44 - 6 [(-0.00195 - 0.6(0.263)] = 2.40 kN m 
M BA = -2.16 - 6[2(-0.00195) - 0.6(0.263)] = -1.19 kN m 
M BC = 1.5 - 20[2(-0.00195) + 0.0195] = +1.19 kN m 

M CB = -1.5 - 20[-0.00195 + 2(0.0195)] = -2.24 kN m 
Mc D = -10[2(0.0195) - 0.263] = 2.24 kN m 
M DC = -10(0.0195 - 0.263) = 2.44 kN m 

3.5 SWAY EQUATIONS 

In simple frames the sway equation was obtained by considering the horizontal 
static equilibrium of the frame. In more complicated frames, such as two or more 
storey frames or gabled frames, the equilibrium of the appropriate parts of the 
frame must be examined. To achieve this the number of independent sway 
modes must be identified. This may be done by replacing the rigid joints by 
hinged joints; the number of independent sway modes is then equal to the 
number of independent kinematic mechanisms. Since an independent mechanism 
corresponds to the action of different loading systems, each mechanism then 
corresponds to an independent equation of equilibrium. The sway equations are 

therefore obtained by considering the equilibrium of the parts of the frame 
corresponding to the independent mechanism in turn. 

A two-storey frame as shown in Fig. 3.7 is used to illustrate the derivation of 
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I I I I □ 




* H 

(a) Actual structure 




»? A 




M 



B 



'AC 



CA 



i 




DB 




*CA 



t 




(c) Sway 




Figure 3.7 
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sway equations. The frame has a total of six unknown displacements, namely 
6 A , 0 B > 0 C , 0 D , A B and A D . The four joints furnish four equations of equilibrium 
and the additional two equations of equilibrium may be obtained from 
independent sway modes. 

The frame has two independent modes of sway identified by the displacements 
A x and A 2 as shown in Fig. 3.7(b) and (c). The two independent sway equations 
are obtained by considering the equilibrium of the portion of the frame corres- 
ponding to each mode of sway. 



Sway I 

In column AC, from SM A = 0 

_ M AC +Mc A 
HcA ~ h, 

In column BD, from 2A/ B = 0 

_ M BD +M PB 

From 2// = 0 of frame ABCD, 
#ca +#db =p i 

or -}-(M AC +M CA +M BD +M DB )=P l [3.14] 



Sway II 

In a similar manner, the relationship between the lateral forces on the whole 
frame and the end moments is 

^(M CE +M EC ) + ^(M DF +M FD )=/^ +P 2 [3.15] 

For the case of a gabled frame shown in Fig. 3.8(a), there are two independent 
mechanisms corresponding to two sway modes as indicated in Fig. 3.8(b) and (c). 

"/'DC = '/'CD =+^i 

2a 

^AB = ^BA =+—^l 
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Since beams BC and CD are equal in length, 

'/'bc = ^cb = -'/'cd ■ -'/'dc = -^1 
Similarly for Sway II 

'/'BC = '/'CB = "'/'CD = "'/'DC = 

2a 

0DE =^ED =+ y^2 

Some of the slope deflection equations are: 



Af AB = -2EK AB (« B -3y h 



M BC = -2£-A: BC [20 B + do ~ 3(-0, + i// 2 )] 
VWcd = -^cd [(2»c + 0d - 3(<//, - <// 2 )] 



A/ DE = -2FA: DE ( 20 D - 3 y * 2 



The sway equations are determined from the free-body diagrams shown in 
Fig. 3.8(d) and (e). 

The equilibrium equations corresponding to Sway I are: 
From 2M B = 0 for member AB, 



_ M AB + M BA 
Ha~- - 

From SM C = 0 for member CD, 



[3.161 



V D L\ 



H D =-^ (mc D + M DC -^-j [3.17] 

From 2M B = 0 for member ABCD 

\ (VL 

V ^ = -l\Y +M dc-M ba ) [3.18] 

Substituting [3.17] and using "LH = P + H A + H D = 0, the sway equation corres- 
ponding to Sway I is 

//_ VL = Mab +A/ba i Mba +Mdc i McD [3 191 

4a b 2a a 

In a similar manner for Sway II, 

VL _ M BC +M DE | M DE +M ED | Mq B 
Aa 2a b a 
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\ L/2 ■ | ■ L/2 ■ | 

(a) 





i.e. 



4 S Wi 




(b) Sway I 



(c) Sway 











| 1 


L/2 


| L/2 


(d) 


L/2 | 


(e) 





Figure 3.8 



74 



THE SLOPE DEFLECTION METHOD 

EXAMPLE 3.4 Determine the joint moments of the gable frame of Fig. 3. 9 
using slope deflection equations. 

Relative Stiffness 



Member 


AB 


BC 


CD 


DE 


Moment of Inertia 


/ 


27 


27 


7 


Length (m) 


6.0 


10.2 


10.2 


6.0 


I/L 


m 


7/5.1 


7/5.1 


7/6 


Relative K 


1.70 


2.0 


2.0 


1.70 



c 




j m IQ.Om | m 10.0 m j 

(a) 



-4.50 




(b) 

Figure 3.9 
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Member Rotations 



Sway I: 



2x2 

0AB = ^ba = + ~T~ 01 = +0-6670, 



0BC = 0C B = -01 



0CD - 0DC ~ 01 



Sway II: 



0BC - ^CB " 02 
0CD = 0DC = -02 

2x2 

0DE = 0F.D = ~~7 02 = 0.66702 



Slope Deflection Equations 

Since supports A and E are fixed, 0 A = 0 E = 0 

A/ AB = -1.7(0 B -3x0.6670,) 
^ba = -1-7(20 B - 3 x 0.6670,) 
M BC = -2.0 [20 B + 0 C - 3(-0, + 0 2 )] 
Mqb = -2.O[20 C + 0 B - 3(-0, + 0 2 )] 
M CD = -2.O[20 C + 0 D - 3(0, - 0 2 )] 
M DC = -2.O[20 D + 0 C - 3(0, - 0 2 )] 
M DE = -1.7(20 D - 3 x O.6670 2 ) 
M ED =-1-7(0 d -3 x O.6670 2 ) 



Sdfor Equilibrium Conditions 

(i) Joint B: A/ BA +M BC =0 

(ii) Joint C: Aic B +M CD =0 

(iii) J oint D : M DC + M DE = 0 
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Shear Conditions 



Using [3.19] and [3.20]: 

1Q = M AB + ^BA , M BA + M D C | M CD 

6 2x2 2 

, v _ ^BC + ^DE | M)E + A*ED ■ M CB 

K) 2x2 6 2 

Substituting the moment expressions into the equilibrium equations gives the 
following simultaneous equations: 

-7.40 B - 2.O0 C - 2.60! + 60 2 = 0 
-2.O0 B - 8.O0 C - 2.O0 D = 0 

[3.21] 

- 2.O0 C - 7.40 D + 60, - 2.60 2 = 0 

1O.20 B + 150 c + 120 D - 38.90! + 270 2 = -60 

120 B + 150 c + 1O-20D + 270i - 38.902 = 0 
Solving the five simultaneous equations, 

0 B =3.350 

0 C = - 1 .948 
0 D = 4.443 
0! = 7.804 
02 =6.864 

End Moments 

Substituting the computed values of rotation into [3.21] gives 
M AB = -1.7(3.350 - 3 x 0.667 x 7.804) = +20.84 kN m 
Mba = —1 -7(2 x 3.350 — 3 x 0.667 x 7.804) = +15.14 kN m 
M BC = -2.0[2x 3.350- 1.960 - 3(-7.804 + 6.864)] = -15.14 kN m 
Mc B = -2.0[-2 x 1.948 + 3.355 - 3(-7.804 + 6.864)] = -4.55 kN m 
M CD = -2.0[-2 x 1.948 + 4.465 - 3(7.804 - 6.864)] = -4.55 kN m 
M DC = -2.0[2 x 4.443 - 1.960 - 3(7.804 - 6.864)] = -8.23 kN m 
M DE = -1-7(2 x 4.443 - 3 x 0.667 x 6.864) = +8.23 kN m 
M ED = -1.7(4.443 - 3 x 0.667 x 6.864) = +15.78 kN m 

The final bending moment diagram is shown in Fig. 3.9(b). 
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3.6 PROBLEMS 



3.1 Find the joint moments of the box structure shown in Fig. P3.1 
(Am: M A =M D = -0.24 kN m 
M B =M c = -0.88 kN m) 




4 m 



Figure P3.1 



3.2 Determine the joint moments of the frame shown in Fig. P3.2 

n IkN/m n 

1 1 I I 1 1 1 I I 1 1 1 1 I 1 1 1 




21 



5kN 



E 



21 



V V//. 



4m 



6m 



Figure P3.2 



^ 



3.3 Draw the bending moment diagram of the frame shown in Fig. P3.3 
(Arts: M A = -6.60 kN m 
M B = +6.08 kN m 

Mq = -10.09 kN m and -0.09 kN m 
M D =+7.19kN m) 
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5kN 



T 



B 



21 



IOHN 



21 



6 m 



Figure P3.3 



3.4 Find the joint moments of the frame shown in Fig. P3.4. 




lOkN 



41 



B 

21 



h — * 



Figure P3.4 



3.5 Find the joint moments of the gable frame shown in Fig. P3.5 
(4ns: M B =+4.66kNm 

M B = -6.15 kN m and -8.15 kN m 
M c = +4.45 kN m) 

I kN/m 

II II H 1 1 I 1 I HI 1 1 1 1 1 M 1 1 1 [I I 1 1 I 1 1 I 





v///////. 



o 

* 



2m 




6 m 




6 m 



2 m 




Figure P3.5 
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3.6 Draw the bending moment diagram of the gable frame shown in Fig. P3.6 




lOkN 



El 



Constant 

E V777?. 




4 m 




Figure P3.6 
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4. The Cross Method of 
Moment Distribution 



4.1 INTRODUCTION 

The moment distribution method was developed by Hardy Cross and was 
formally presented in 1930. The method became dominant for continuous 
beams and rigid frame analysis for the next thirty years. The method remains 
very important for manual solutions of frames. The method of solution places 
emphasis on physical behaviour, which is performed in tabular form, and is quick 
and easily remembered. Because of the nature of the concept and procedure, 
which readily lends itself to the development of design judgement, the method 
has transformed engineering thinking and design office practices. 

The moment distribution method starts from the same basic assumption made 
in the slope deflection method. In the analysis of continuous beams and frames 
all joints are assumed fixed and the moments are then corrected by means of the 
slope deflection equations. While in the slope deflection method the correction 
is effected by solving a set of linear simultaneous equations, in the moment 
distribution method the correction is dealt with by successive approximations 
through the use of arithmetic only where the results can be obtained to any 
desired accuracy. 

4.2 ITERATIVE SOLUTION OF SLOPE DEFLECTION EQUATIONS 

The slope deflection equation involves the formulation of equations of joint 
equilibrium in terms of rotations, displacements, stiffness and length of members. 
These equations of equilibrium form a set of linear simultaneous equations which 
may be solved by methods such as elimination methods, iteration methods, or 
relaxation techniques. 

One of the iteration methods for obtaining the solution of a set of simultane- 
ous equations is by the Gauss— Seidel method. To illustrate the method, consider 
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the problem of Example 3.1. The governing equations of equilibrium are 



One means of obtaining an approximate solution to the equations is to make 
a guess. Substituting the assumed values in the right-hand sides of [a] and [b] 
furnishes a new set of values on the left. If the first assumed values do not work, 
the new set of values are selected as a better guess and the process is repeated 
until a consistent set of values are obtained that satisfy the equations. To 
furnish a more rapid convergence of the iterative process, the Gauss— Seidel 
method presupposes the most recent cycle of values at each stage of calculations. 
This method is explained below by solving the above equations. 

Cycle 1 Assume 0 C = 0 in [a] and solve; 0 B = -2.462. Using 0 B = -2.462 in 
[b], solve; 0 C =-7.769. 

Cycle 2 Using the new value of 0 C = —7.768 in [a] obtain a new value for 
0 B (-0.069). Use 0 B = -0.069 in [b] to solve; 0 C = -8.965. 

This procedure continues, as shown in Table 4.1, for a sufficient number of 
cycles until two successive cycles agree within acceptable error-limits. The final 
results are 

0 B = +0.366 

0 C =-9.183 

which agree closely with those obtained in Example 3.1 
Table 4. 1 



Cycle 


0B 


0 C 


Initial 


0 


0 


1 


-2.462 


-7.769 


2 


-0.069 


-8.965 


3 


+0.299 


-9.150 


4 


+0.356 


-9.178 


5 


+0.365 


-9.182 


6 


+0.366 


-9.183 


7 


+0.366 


-9.183 



260 B + 80 c = -64 
80 B + 160 c = -144 



The above equations may be written as 
0 B = -2.462 - O.3O80 C 



[a] 
[b] 



0 C = -9.0 - O.50 B 
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4.3 INTERPRETATION OF THE ITERATIVE SOLUTION 



si 



Each iterative cycle shown in Section 4.2 has a physical interpretation on the 
beam of Example 3.1. It is recalled that the slope deflection equations are 

M AB =16O-50 B 

M BA =-16O-1O0 B 

A/ BC =96 - 160 B -0 C 
M C b = -144 - 80 B - 160 c 
and the equilibrium equations are 

^BA + ^BC ~ 0 

Mcb=0 



[c] 



[d] 

[e] 



The values of 0 B and 0 C for each step of the iteration scheme as evaluated in 
Table 4.1 are given in Table 4.2. Also shown are the moments of the slope 
deflection equations corresponding to the values of 0 B and 0c evaluated in 
Table 4.1 



Table 4.2 



Cycle 


013 


6c 




Mba 




Mcb 


Initial 


0 


0 


160 


-160 


96 


-144 


1 


-2.462 


0 

-7.769 


172.31 


-135.38 


135.39 
197.54 


-124.30 
0.0 


2 


-0.069 


-7.769 
-8.965 


160.35 


-159.31 


159.26 
168.82 


-19.14 
-0.01 


3 


+0.299 


-8.965 
-9.150 


158.51 


-162.09 


162.94 
164.42 


-2.95 
+0.01 


4 


+0.356 


-9.150 
-9.178 


158.22 


-163.56 


163.50 
163.73 


-0.45 
0.00 


5 


+0.365 


-9.178 
-9.183 


158.18 


-163.65 


163.58 
163.62 


-0.07 
-0.01 



A close review of the joint rotations and the corresponding support moments in 
Table 4.2 reveals the following: 

(a) At the initial stage none jof the equilibrium equations, [d] or [e], is satisfied. 

83 



METHODS OF STRUCTURAL ANALYSIS 



(b) In the first phase of Cycle 1, joint B is permitted to rotate while joint C 
remains clamped. Consequently, [d] is satisfied while [e] is not. 

(c) In the second phase of Cycle 1, joint B is clamped in its rotated position 
while joint C is allowed to rotate until Mqb becomes zero, thus satisfying 
[ej. However [d] is not satisfied. 

(d) In Cycle 2, as in Cycle 1 , the same pattern of permitting rotation at one 
support while clamping at the other is repeated. 

(e) The process of permitting joint rotations at one support while clamping the 
other is repeated as many times until the equilibrium equations are satisfied. 
At this stage the correct support moments are obtained. 

This process of balancing operation by alternately permitting and preventing 
joint rotations at each joint until the equilibrium condition is attained is the 
physical process which corresponds to the iterative solution. This physical 
process is known as the moment distribution method, and it is developed 
hereafter. 



4.4 FUNDAMENTAL FACTORS USED IN MOMENT DISTRIBUTION 

Prior to developing the method of moment distribution, certain definitions and 
fundamental relationships need to be considered. 



Rotational Stiffness 

Consider a beam with constant moment of inertia, /, and no interior load 

(Fig. 4.1). If end B is fixed and there is no relative displacement of the ends, the 

slope deflection equation is written as: 




A 



h L - H 

Figure 4.1 
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M AB = 2EK{26 A ) 



or 



^-4EK [4.1] 

Since the value of M AB /d A is a measure of the resistance of the beam to 
rotation, it is defined as the stiffness of the member. 

When a beam has its far end hinged instead of fixed, less moment is required 
to rotate the end through a given angle. Figure 4.2 is simply supported at the far 
end and a moment M AB is applied at the near end. By the slope deflection 
method, 

M AB = -2EK(29 A - 6 B ) 

M BA =-2EK(2d B +9 A ) = 0 



or 



e B = -8 A /2 



Thus 



M AB = 3EKd A [4.2] 

It is noted that the stiffness of the beam at the near end when the far end is 
hinged is three-fourths of the stiffness when the far end is fixed. The significance 
of this ratio is that it permits the establishment of distribution factors consistent 
with true end hinged taken as 3EI/L; the far end can be released in the first 
balancing operation and then left free during subsequent calculations, or, no 
moment is carried over to the hinged end. The adjusted relative stiffness makes 
this unnecessary and the technique substantially reduces the calculation. 





B 



Figure 4.2 



Carry-Over Factor 

The slope deflection equation gives the corresponding moment at the far end, 
Fig. 4. 1 (support B), as 

M BA = 2EK6 A 
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The ratio of the induced moment M BA to the applied moment M AB 

m ba _ 1 
M AB 2 

is defined as the carry-over factor. 
Distribution Factor 

If there are several members framed into a joint, and if a moment is applied to 
joint A as shown in Fig. 4.3, then 

M A b + M AC + M\d -M Q = 0 

or, 

M 0 = 4Ed A (K AB + K AC + K AD ) 
= 4E6 A ZA" 

The moments of the members at joint A are 

M AB = 4EK AB d A 

M AC =4EK AC 6 A 

M AD =4EK AD 9 A 
Substituting 4E6 A = M 0 fZK 



M ab =y^M 0 = (DF)abM 0 



M AC =^M 0 =(DF) AC M 0 [4.4] 



M ad =^M 0 = (DF)adM 0 




B 

Figure 4.3 
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Thus the applied moment M Q is distributed to the members in proportion to 
the relative stiffness K. The term KfZK is known as the distribution factor, 
which is denoted asDF. It is self-evident that the sum of all distribution factors 
at a given joint must equal unity. 



4.5 MOMENT DISTRIBUTION METHOD FOR BEAM ANALYSIS 

The following example is used to illustrate the how and why of the moment 
distribution method. The statical moment sign convention is used; that is, 
positive moment tends to rotate the joint clockwise or the member counter- 
clockwise. 



EXAMPLE 4. 1 Determine the support moment of the continuous beam in 
Fig. 4.4. The beam has constant moment of inertia I. 

A step-by-step solution is given consisting of (i) locking all joints to prevent 
rotation, (ii) computing fixed-end moments, (iii) distributing unbalanced 
moments in proportion to the relative stiffness, (iv) carrying over one-half of the 
distributed moments to the opposite ends of the members. 



C 

4 m i 




Relative Stiffness Values and Distribution Factors 
The relative values of K = I/L are: 

SpanAB: K AB =7 x 12 = 2 

6 

SpanBC: K BC = t x 12 = 3 

4 

The distribution factors are: 



^AB 


2 


ZK 


2 + oo 


^BA 


2 




2 + 3 



Wba=^ — = 0.4 
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^c=^ = ^ = 0.6 

Note that (DF) AB = 0 since a fixed joint is considered infinitely stiff and 
(D/Ocb = 1 0 since CB is the only member at hinge C. The DF values are 
entered at the top of the calculation scheme in appropriate places. 

Fixed-End Moment 

Assume that all joints are locked against rotation before the loads are applied on 
the beam. The fixed-end moments produced by the loads are 

M B F A = 12X . 3 2 X 32 =+9.0kNm 

o 

ML = -9.0 kN m 



Balancing of Joint Moments and Carry-Over Moments 

Distribute the unbalanced moment +9.0 at support B in proportion to the DF 
values by applying +3.6 on member BA and +5.4 on member BC. Consequently, 
half the moment of 3.6 kN m is carried to A while half the moment of 5.4 kN m 
is carried to C. Since support A is fixed, any moment carried over to A remains 
as carried over. However, joint C is released by applying -2.7 kN m to cancel 
the moment of +2.7; this results in zero moment at C which must be the case for 
a hinge. As C rotates under the applied moment, B remains fixed and the 
unbalanced moment at B is —1.35 kN m. The procedure of distributing and 
carry-over is repeated until all joints are balanced to the desired accuracy. After 
the final balancing operation the end-moments are algebraically added to get the 
final end moments shown in Table 4.3. 
The support moments are: 

M AB = 11.11 kNm 

M BA = -4.77 kN m 

M BC = +4.77 kN m 

Mc B =0 

(see the numbers in italic in Table 4.14) 

The balancing computations may also be performed in tabular form as given 
in Table 4.4. 
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Table 4.3 



Relative K 



Distribution Factor 



0 



0.4 



0.6 



1.0 



Fixed-end moment 

Distribute 

Carry-over 

Distribute 

Carry-over 

Distribute 

Carry-over 

Distribute 

Carry-over 

Distribute 

Carry-over 

Distribute 

Carry-over 

Distribute 



+9.0 



+ 1.8 



+0.27 



+0.04 



-9.0 
+3.6 



+0.54 



+0.08 



+0.01 




-0.20 
+0.12 



-0.03 
+0.02 



+0.40 
-0.40 



+0.06 
-0.06 



Final Moment 



+ 11.11 



-4.77 



-4.77 



0.0 



Table 4.4 



Joint 



B 



Member 



AB 



BA 



BC 



CB 



K 



DF 



0.4 



0.6 



1.0 



Fixed-end 
moment 



+9.00 
+ 1.80 



+0.27 



+0.04 



-9.00 
+3.60 



+0.54 



+0.08 



+0.01 



+5.40 



-1.35 



+0.81 



-0.20 



+0.12 



-0.03 



+0.02 



+2.70 



-2.70 



+0.40 



-0.40 



+0.06 



-0.06 



Total 



+ 11.11 



-4.77 



+4.77 



0.00 
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EXAMPLE 4.2 Determine the support moments of the beam given in Example 
4. 1 using the simplified treatment of a hinge. 



Relative Stiffness Values and Distribution Factors 

_ 12 _ 
* AB = 6 = 2 

*bc=^=2.25 



(WOba = 



2 + 2.5 



= 0.47 



The moment distribution is performed in tabular form as shown in Table 4.5 
Table 4. 5 



Joint 


A 


B 


C 


Member 


AB 


BA | BC 


CB 


K 


2 


2 


2.5 


2.5 


DF 


0 


0.47 


0.53 


1.0 


Fixed-end 
moment 


+9.0 
+2.11 


+4.23 


+4.77 




Total 


+ 11.11 


-4.77 


+4.77 


0.0 



For beams with a loaded overhang at the hinged end the stiffness is 3EI/L 
since the stiffness of the cantilever is zero. When there is no load on the cantilever 
the moment at the hinged end will be zero, but if a load is applied on the 
cantilever the bending moment will be statically determinate and is independent 
of loads on adjacent spans. A joint with an overhang, therefore, may be treated 
as a beam with a hinged end as far as the moment distribution is concerned. 



Beams with Relative Displacement of Supports 

The moment distribution may be used to analyse continuous beams with relative 
displacement of supports or rotation of supports. The effects of relative 
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displacement of supports can be considered by referring to Fig. 4.5. The beam 
AB is fixed against rotation at A and B. The end B is deflected downward an 
amount A while both ends are restrained against rotation. 

The slope deflection equation is used to determine the value of Af AB . 

M AB »Jfe A =-2£K(-3A/Z.) 

6EKA 



The corresponding /zxed-erttf moments are 
M F - m f _ 6g/A 

M AB "MBA - ^ 2 

Similarly, the fixed-end moments due to rotation 0 A at A are 

AEId K 



[4.5] 



<B = 2Mf A =-^ [4.6] 
A beam with one end hinged which deflects A is shown in Fig. 4.6. 




Figure 4.6 



The slope deflection equation gives 
M BA = -2EK(2d B - 3A/L) = 0 
0 B = 3A/2L 
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Thus 



( 

= -2EK\d B -- 



3EKA 



The corresponding fixed-end moment is 

m f -^L 

AB " /2 

[4.7] 

M BA =0 

The appropriate fixed-end moment corresponding to the relative displace- 
ments will then be introduced into the moment distribution scheme. These 
moments are treated in exactly the same manner as those due to the applied 
loadings. 



4.6 MOMENT DISTRIBUTION METHOD FOR FRAME ANALYSIS 



In applying the moment distribution method to frame problems, there is no 
inherent limitation which prevents the use of the basic concepts developed so 
far. The method can be applied to frame problems with some changes in the 
book keeping scheme. For frames with fewer number of joints, the tabular 
scheme of calculations may be more convenient to adopt. However, for frames 
with numerous joints the sketch of the structure itself may be more practical 
to use. One system of recording the moments is to place the computed value 
under the beam on the left and over the beam on the right. Similarly, the column 
moments are recorded to the left on the top of the column and to the right on 
the bottom. This arrangement of showing the calculations on the structure 
itself is shown in Fig. 4.7. 

Frame problems are classified in two categories: 

(a) Frames without sidesway 

(b) Frames with sidesway 



4.6.1 Frames Without Sidesway 

Frames without sidesway are analysed in the same manner as continuous beams. 
In the case of frames there are frequently more than two members meeting at a 
joint so that the joint moments distribute among all members according to the 
appropriate distribution factors. Consequently, for such frames there is a need of 
a more systematic arrangement of computations. 
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T 








T 


■j] — 








JL 





Figure 4.7 



EXAMPLE 4.3 Determine the joint moments of the frame shown in Fig. 4.8. 



Relative Stiffnesses and Distribution Factors 



K A .=f-4 

2x 12 
K BC =—- - = 4 



3 1-5 x 12 „„„ 
At E =-x— =3.375 

Wba=^ = 0.5 



93 



METHODS OF STRUCTURAL ANALYSIS 



(PF)BC = 



4 + 4 



= 0.5 



»B = 



4 + 4 + 3.375 



= 0.352 



11.375 



= 0.352 



w - = i&h - 0296 



T 




8 



3 kN/m 




8kN 


lllllllllllllllllllllllllllll 






27 


C '51 





6 m 




4m 



Figure 4.8 



Fixed-End Moments 

* 3(6) 2 
Mf c = -^B = -^- = 9kNm 

F 8(3)(1) 2 
A/ce = — ^2 = 1.5 kN m 

„ fc . sage . kN m 

16 

The moment distribution is perfonned in tabular form since the number of 
joints is few. 
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Table 4. 6 



Joint 


A 


B 


c 


D 


E 


Member 


AB 


BA 


BC 


CB 


CE 


CD 


DC 


EC 


K 


4 


4 


4 


4 


3.375 


4 


4 


3.375 


DF 


0 


0.5 


0.5 


0.352 


0.296 


0.352 


0 




Fixed-end 
moment 


— 2.4o 
-0.11 


-4.96 


+9.00 

+0.92 
-4.96 


-9.00 
+ 1.85 


+ 1.50 
+2.25 

+ 1.55 


+ 1.85 


+0.93 
+0.44 
+0.02 


-4.50 
+4.50 


-2.48 
+0.87 


+0.74 


+0.87 


-0.22 


+0.44 
-0.11 


+0.04 


+0.03 


+0.04 




+0.02 


Total 


—2.59 

1 


-5.19 


+5 19 


—8.83 


+6.07 


+2.76 


+ 1.39 


0.0 



4.6.2 Frames With Sidesway 
Rectangular Frames 

When a frame is loaded laterally, or in the case when the loading or the frame 
itself is unsymmetrical, sidesway or joint translations occur. As the joint 
displacements are unknown the fixed-end moments due to the displacements 
cannot be calculated- The solution for frames with a single mode of sway is 
indicated in Fig. 4.9. In applying the moment distribution method the joints are 
first assumed to be held against sidesway by introducing artificial restraint R, at 
the appropriate joints. The fixed-end moments caused by the applied loads are 
then distributed to obtain the non-sway balanced end moments. Next, the 
magnitude of the reaction./?, at the artificial restraint is determined from the 
considerations of equilibrium of the members. The effect of the restraint R is 
determined by permitting an arbitrary sway to take place, unaccompanied 
initially by rotation at the joints. These fixed-end moments are then distributed 
and subsequently, the force F necessary to maintain the frame in its swayed 
position is determined from the equilibrium condition of the column shears. The 
actual magnitude of lateral force, F, consistent with the condition necessary to 
eliminate the artificial restraint R 9 and therefore the moments caused by the 
sidesway, are determined from the condition: 

R = kF [4.8] 
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B 



R, a 




1 



1 



1 



(a) 






kF A. 



(b) 

Figure 4.9 



+ Ar x |J 
II 



(0 




The final set of end moments are obtained by adding the original set of 

moments to the sway moments multiplied by the correction factor k. This 

process is illustrated by the frame shown in Fig. 4.9. 

Note that the frame shown in Fig. 4.9(a) is equivalent to the sum of Fig. 

4.9(b) and Fig. 4.9(c). In Fig. 4.9(b) the frame is prevented from sidesway by the 

artificial restraint/? applied at joint C. The arbitrary displacement A occurs at C 

with the joints held against rotation and when the fixed-end moments are 

distributed the horizontal force that caused the sidesway F is determined from 

equilibrium conditions. The shear condition of the given frame is 



H' K +H{>-P 1 =R 



But 



//; = 



^ab+^ba , P 2 a 



„, _ M C T> + M DC 
D " h 



Thus 



h\ hi h 2 



Similarly 



or 



^AB + ^BA , McD +Mq C 

hi h 2 



= F 



[4.9] 



But 



R=kF 



[4.10] 
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Therefore by superposition condition 
^ab =Af AB + Wab 
M$a =^ba +*Mba. etc. 



[4.11] 



EXAMPLE 4.4 Determine the joint moments of the frame shown in Fig. 4. 1 0. 



1.5m 



T 



B 



o 

CM 



E 

Q 

ro 



3kN 



4KN 



E 

O 



3.0m 



(a) 



3 .0 



2/2 




154 



154 





YA'1.07 



146 




12 



(b) 



Figure 4.10 



Relative Stiffnesses and Distribution Factors 



^CD ~ "J" ~ 5 



0.48 




0.24 




V B =0.24 



024 
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10 + 3 



= 0.231 



10 



(WObc-jj- 0-769 



CD/0 



CB 



10 



10 + 5 



= 0.667 



Wcd = ^ = 0.333 



Fixed-End Moments 



M 



^CB " ~^BC _ 



25 



F _ 



kN m 



4(1.5)(1.5) : 



= 1 .50 kN m 



M Moment Distribution without Sides way 



Table 4. 7 Distribution without Sidesway 



Joint 


A 


B 


C 


D 


Member 


AB 


BA 


BC 


CB 


CD 


DC 


K 


3 


3 


10 


10 


5 


5 


DF 


0.0 


0.231 


0.769 


0.667 


0.333 


0.0 


Fixed-end 


+ 1.44 


-2.16 


+ 1.5 


-1.5 






moment 






+0.5 


+ 1.0 


+0.50 


+0.25 




+0.02 


+0.04 


+0.12 


+0.06 














-0.04 


-0.02 


-0.01 


Total 


+ 1.46 


-2.12 


+2.12 


-0.48 


+0.48 


+0.24 
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By using the end moments on the free-body diagrams as shown in Fig. 4. 10(a) 
the end shears are determined to be 

3x2+1.46-2.12 , n ^ n , V1 
V A = = 1 .068 kN (left) 

Kb . . 0 . 240 kN 0eft) 

The artificial joint restraint is 

R =3.0- 1.068-0.24= 1.692 kN (left) 

(b) Moment Distribution with Sides way 

Assume the frame to sway an arbitrary amount A such that the fixed-end 
moments in the columns are 



6£YA 



6£YA 



= 0.24£Y 



6£YA 



Mcd = +M£ C = — = 0.661EI 



Taking EI = 10, distribution of the fixed-end moments due to sidesway is 
shown in Table 4.8. 



Table 4.8 Distribution with Sidesway 



Joint 


A 


B 


C 


D 


Member 


AB 


BA 


BC 


CB 


CD 


DC 


DF 


0.0 


0.231 


0.769 


0.667 


0.333 


0.0 


Fixed-end 
moment 


+2.4 
-0.27 


+2.4 
-0.55 


-1.85 


-0.93 


+6.67 


+6.67 








-1.92 


-3.83 


-1.91 


-0.95 






+0.44 


+ 1.48 


+0.74 












-0.24 


-0.49 


-0.25 


-0.12 




+0.03 


+0.06 


+0.18 














+0.03 


-0.06 


-0.03 




Total 


+2.38 


+2.36 


-2.36 


-4.48 


+4.48 


+5.59 
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The column shears are 

K AB = ^P = 0.948 ,N 

Koc= 4 -^ = 3.356 ,N 

The net shear is 

F=V AB + V DC 
= 0.948 + 3.356 =4.305 kN 

The correction factor is 

, R 1.692 A „ rt „ 

k -?-^rs- +om 

The final end moments are obtained by multiplying the moments from the 
sway solution and then adding the moments from the non-sway solution to 
obtain the moments in the original structure. 



Thus, 



M AB = 1.46 + 2.38(0.393) = +2.40 kN m 
M BA = -2.12 + 2.36(0.393) = -1.19 kN m 
Mcb = -0.48 - 4.48(0.393) = -2.24 kN m 
M DC = 0.24 + 5.59(0.393) = 2.44 kN m 



Frames with Inclined Members 

The superposition method described for the case of rectangular frames with 
sidesway, can also be used in the analysis of frames with inclined legs. For 
example, for the single storey portal frame with inclined legs shown in Fig. 4.1 1, 
the frame is held against sidesway by artificial joint restraint R (Fig. 4.11(b)) 
and, after the moment distribution calculations are carried out, the value of R is 
determined. In the second analysis, the consistent joint force F is computed due 
to an arbitrary lateral displacement A. 

After determining the consistent joint force F, the correction factor k is then 
obtained to calculate the final moments. 
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METHODS OF STRUCTURAL ANALYSIS 

EXAMPLE 4.5 Find the joint moments of the frame shown in Fig. 4.12 by the 
moment distribution method. 



Relative Stiffness and Distribution Factors 



30 

*ab= y-10 



^BC ~ ~, = 1 5 



^cd = T" = 6 
Wba =^ = 0.4 
(PF)bc =^ = 0.6 
(Wcb-JJ- 0-714 
Wcd = |f = 0.286 



Fixed-End Moments 



< c =-M&=^P=5.0kN 



faj Moment Distribution without Sidesway 



The member shears are 

1.51+3.04 
*ab = ; =1.52 kN 



Kdc . Mlll^ . 0 . 62 kN 



Figure 4.12(b) shows the forces and reactions on the frame. The artificial 
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lOkN 




Figure 4.12 



joint restraint/? may be found by taking moments about the point of inter- 
section of members AB and DC. Thus, 

5.33/? + 10(2) + 0.62(1 1.67)+ 1.51 - 5(5.33)- 1.52(8.33)- 1.02 = 0 
K = 2.17kN 
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Table 4. 9 Distribution without Sidesway 



Joint 


A 


B 


C 


D 


Member 


AB 


BA 


BC 


CB 


CD 


DC 


K 


10 


10 


15 


15 


6 


6 


DF 


0.0 


0.40 


0.60 


0.714 


0.286 


0.0 


Fixed-end 
moment 


- 1 .0 


-2.0 


+5.0 
-3.0 


-5.0 
-1.5 












+2.32 


+4.64 


+ 1.86 


+0 93 




-0.46 


-0.93 


-1.39 


-0.70 












+0.25 


+0.50 


+0.20 


+0.10 




-0.05 


-0.10 


-0.15 


-0.08 












+0.03 


+0.06 


+0.02 


+0.01 






-0.01 


-0.02 








Total 


-1.51 


-3.04 


+3.04 


-2.08 


+2.08 


+ 1.04 



(b) Moment Distribution with Sidesway 

The loads are removed and the frame is permitted to sway to the right through a 
horizontal distance A. As shown in Fig. 4.12(c) the relative displacements are 

Member AB : 1 .OA 
Member BC : 0.75A 
Member CD : 1.25A 



Fixed-End Moments 



F _ 6EIA _ 6£7A 
~~ " = (3) 2 



p _ 6£(27)(0.75A) 

(4) 2 



M BC 



p _ 6EI(\. 2 5 A) 
(5) : 
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Assuming EI A = 1 0 

Mg B = +6.67 kN m 
Mg c = -5.63 kN m 
Mqq = +3.00 kNm 



Table 4.10 Distribution with Sidesway 



Joint 


A 


B 1 


C 


D 


Mp mhpr 

1 VI W ill UWl 


AR 


BA 


BC 


CB 


CD 


nr 


DF 


0.0 


0.40 


0.60 


0.714 


0.286 


0.0 


Fixed-end 


+6.67 


+6.67 


-5.63 


-5.63 


+3.00 


+3.00 


moment 


-0.21 


-0.42 


-0.62 


-0.31 












+ 1.05 


+2.10 


+0.84 


+0.42 




-0.21 


-0.42 


-0.63 


-0.31 












+0.11 


+0.22 


+0.08 


+0.04 




-0.02 


-0.05 


-0.07 


-0.03 














+0.02 


+0.01 




Total 


+6.23 


+5.78 


-5.79 


-3.94 


+3.93 


+3.46 



The member shears are 

r 6* + SM mAMW ( , eft) 



The consistent joint force F required to produce the set of moments given in 
Table 4.10 is found by taking moments about the point of intersection of AB and 
CD. Thus, 

5.33F+ 6.23 + 3.46 - 4(8.33) - 1.48(1 1.67) = 0 
F=7.67kN 
The correction factor is 

*-§-!£- 0383 
F 7.67 
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Final Moments 



The final moments are determined by algebraically adding the results of the 
no-sway solution to the products of k and the corresponding results of the sway 
solution. Thus 

M AB = -1.51 + 0.283(+6.23) = +0.25 kN m 
^ba =-^bc = -3.04 + 0.283(+5.78) = -1.40 kN m 
Mqb = -^cd = -2.08 + 0.283(-3.94) = -3.20 kN m 
M DC = +1.04 + 0.283(+3.46) = +2.02 kN m 



4.6.3 Frames With Multiple Degrees of Freedom With Respect to Sidesway 

A rigid frame which has n independent joint translations is said to have n degrees 
of freedom with respect to sidesway. For example, in the case of the two-storey 
frame shown in Fig. 4.13(a), with joint B deflecting A B and joint C an independ- 
ent deflection of A c , the frame is said to have two degrees of freedom with 
respect to sidesway. 

By applying the principle of superposition, the two-storey frame (Fig. 4.13) 
may be analysed in three separate steps. First, the frame is completely prevented 
from sidesway by introducing artificial supports as shown in Fig. 4.13(b). In this 
case all the given loads are applied and a regular no-sway moment distribution is 
carried out to obtain the artificial joint restraint (assumed positive to the left) 
R\o and R 20 ln the second step, a translation of the frame of an arbitrary 
displacement A, is introduced at joint C (Fig. 4.13(c)). While translation is 
introduced the joints are locked against rotation and initial moments are 
developed in members AB and EF as shown by the solid lines. To permit the 
joints to rotate, shown by the dotted elastic curve, a moment distribution is 
performed and the consistent joint forces R { { and R x2 are calculated. Finally, 
a similar solution is conducted for an arbitrary displacement A 2 at joint C 
(Fig. 4.13(d)). 

With the separate trial solutions, it now remains to determine how much of 
the two sidesway solutions should be superimposed to the first case to obtain 
the final results for the actual given problem. It must be possible to find the 
final values by obtaining multiplying factors k x and k 2 from a linear combination 
of the obtained results. These factors are obtained by solving two simultaneous 
equations formulated from the superposition equations for the reactions at D 
and E: 

[4.12] 

R 20 +k x R 2{ + k 2 R 22 =0 
After finding the values of the proportionality factor k x and k 2 , the moments in 
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Figure 4.13 

the frame are determined from the equation 

M = M 0 +k 1 M l +k 2 M 2 [4.13] 

where M 0 represents the moments from the non-sway moment distribution, and 
M x and M 2 are the moments due to sidesway, shown in Fig. 4. 1 3(c) and (d) 
respectively. 
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The procedure described above can also be used in the analysis of a frame 
having n degrees of freedom with respect to sidesway. For such a frame there 
will be one non-sway moment distribution and n separate sidesway cases, in each 
case only one independent sidesway is permitted. After completing the required 
(n + 1) moment distribution analyses, the n superposition equations for the 
artificial joint restraints and consistent joint are 

^10 
^20 

§ * ■ 

.^30 

or the simultaneous equations may be written in index notation as 

n 

*/o+L = 0 /=l,2,...n [4.15] 
/=l 

The solution of the above simultaneous equations gives the values of the multi- 
plying factors k x , k 2 , . . . k n . These values are used to find the final moments. 

M = M 0 + k x M x + k 2 M 2 + . . . + k n M n 

n 



Split-Level Frames 

The frame of Fig. 4.14 has two degrees of freedom with respect to sidesway. 
The horizontal displacements of the two roof levels are designated as A B and A D 
as shown in Fig. 4.14(a). Since A B and A D are unknown, two separate solutions 
must be obtained where independent sidesways are permitted as shown in 
Fig. 4.14(c) and (d). With the consistent joint forces R 10 , R 1 2 , R 2 1 and R 22 
determined, the proportionality factors k x and k 2 are obtained from equilibrium 
equations as described above. 

Gabled Frames 

Gabled frames of a single span have two degrees of freedom with respect to joint 
translation; accordingly, two artificial joint restraints are required to prevent the 
joints from moving. The steps required for the analysis of gabled frames by 
moment distribution are shown in Fig. 4.15. 
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R n R 



R 



21 



R 



12 



22 



Km 



« | * 



• • • 



R 



nl 



R 



nn . 



k 2 



0 
0 



0 . 



[4.14] 
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10 (d) 

Figure 4.14 

4.7 CANTILEVER MOMENT DISTRIBUTION 

Unlike the conventional moment distribution, the cantilever moment distribution 
permits sidesway to occur during the moment balancing process. The method, 
therefore, makes it possible to evaluate the moments in one distribution without 
requiring artificial joint restraints. Within its range of applicability the cantilever 
moment distribution provides a simple but powerful method of analysing 
symmetric frames subjected to lateral loads at storey heights or to antisymmetric 
loads. The derivation of the basic equations necessary for the development of 
the cantilever moment distribution is given below. 

(a) Stiffness Factor of a Cantilever 

Consider a uniform cantilever beam AB subjected to an end moment A/ab ^ 
shown in Fig. 4.16. 
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A 




Figure 4.15 




Figure 4.16 



Mba 




It is seen that the angle of rotation is given by the expression 



0 AB - 



_M AB L 
EI 



or 



M AB _EI 



9 



AB 



[4.17] 



which is defined as the rotation stiffness of the cantilever beam. Notice that the 
stiffness of beam AB is one-fourth of the same beam with support A hinged. 
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(b) Carry-over Factor 

The bending moment at end B of the beam (Fig. 4.16), in accordance with the 
statical sign convention, is equal in magnitude but opposite in direction to the 
applied moment at support A. The translational carry-over factor is therefore 

-1.0. 



(c) Stifpiess Factor of Beams under Antisymmetric Bending 

Consider a uniform beam (Fig. 4.17) under antisymmetric bending moments at 
the ends. 




Figure 4.17 



The slope-deflection equation for the beam is 

2EI 

M AB =—-(2d A +d B ) 



or 



M AB _ 6EI 



[4.18] 



which is the rotation stiffness of a simply supported beam under antisymmetric 
end moments. 



The steps in applying the cantilever moment distribution method may be 
summarised as follows: 

(a) Evaluate the stiffness value for a member parallel to the axis of symmetry 
from EI/L. The carry-over factor is —1 .0 when the far end of the member is 
fixed. 

(b) The stiffness factor for a member perpendicular to the axis of symmetry is 

6EI/L. 

(c) Compute the fixed-end moments from the condition that the joints are 
locked against rotation but free to translate. For more than one-storey 



111 



METHODS OF STRUCTURAL ANALYSIS 



frames, the fixed-end moments are computed to be directly proportional to 
the product of the storey-shear and storey-height. 

(d) Balance the moments for one-half of the structure. 

(e) Determine the correction factor to satisfy horizontal equilibrium condition. 

(f) Compute final joint moments by multiplying with the correction factor the 
moments obtained in step (d). 



EXAMPLE 4.6 Determine the joint moments of the frame in Fig. 4. 18 using 
(a) the conventional moment distribution method (b) the cantilever moment 
distribution method 



lOkN fc B 



I 



21 



6.0m 



Figure 4.18 



(a) Conventional Moment Distribution Method 
Relative Stiffness Values and Distribution Factors 



I 



1 



E 

o 



^AB ~~ Km - T - 0.25 



KftQ — — - 0.33 



Wba = (DF) CD = 



0.25 



0.25 + 0.33 



= 0.431 



(DF) BC = (DF)cb = ^ = 0.569 



Fixed-End Moments (Relative Values) 

M£ B =Mg A =M c f d =Mg c = ^ = 10.0 kN m 
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Table 4.11 



Joint 


A 


B 


C 


D 


Member 


AB 


BA 


BC 


CB 


CD 


DC 

^^^^ 


a: I 


0.25 


0.25 


0.33 


0.25 


0.25 


0.25 


DF 1 


0.0 


U.tJ 1 


U.JO" 


0 


0 431 




Fixed-end 1 
moment I 


+ 10.0 


+ 10.0 


-2.85 


-5.69 


+ 10.0 
-4.31 


+ 10.0 
-2.15 




-1.54 


-3.08 


-4.07 


-2.04 












+0.58 


+ 1.16 


+0.88 


+0.44 




-0.13 


-0.25 


-0.33 


-0.16 












+0.04 


+0.09 


+0.07 


+0.03 




-0.01 


-0.02 


-0.02 








Total 


+8.32 


+6.65 


-6.65 


-6.64 


+6.64 


8.32 



The base shear is 



V=-^(M AB +M BA +Mcd +M DC ) 



4.0 



(8.32 + 6.65) = 7.485 kN 



The correction factor is 

The final end moments are 

^ab =M DC = 1.336 x 8.32 = 1 1.12 kNm 
^ba =^cd = 1-336 x 6.65 = 8.88 kN m 

(b) Cantilever Moment Distribution Method 
Relative Stiffness Values and Distribution Factors 



K AB ~ K CD ~ ~j~ ~ 4 q _ 0 25 
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6.0 



2.0 



2.0 + 0.25 



= 0.11 



(PF)bc = ^ = 0.89 



Fixed-End Moments (Relative Values) 

A/ab=^ba=^=10.0 kNm 

The distribution is carried out in tabular form as shown in Table 4.12. Notice 
that the fixity at joint B does not exist, it is therefore released by applying at 
joint B a balancing moment of -10.0 kN m. The moment at joint B is distributed 
as -10.0 x 0.89 = -8.9 kN m to member BC and -10.0 x 0.1 1 = -1.1 kN m to 
member BA. The carry-over factor for column BA being -1.0, the moment 
carried over to joint A from B is +1 .1 kN m. 

Table 4. 1 2 



Joint 


A 


B 


Member 


AB 


BA 


BC 


DF 


0.0 


0.1 1 


0.89 


Fixed-end 


+ 10.0 


+ 10.0 




moment 


+1.1 


-1.1 


-8.9 


Total 


-11.1 


+8.9 


-8.9 



The base shear is 



K = — (1 1 .1 +8.9)= 10.0 kNm 
The correction factor is 

The final joint moments are 

Mab = ^dc = 1.0(11.1)= 11.1 kNm 
Mba = M C v = 1 0(8.9) = 8.9 kN m 
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It is noted that the computations are carried out on one-half of the frame 
only since the moments on either sides of the axis of symmetry of the frame are 
identical. 

EXAMPLE 4.7 Find the join t moments of the two-storey frame shown in 
Fig. 4.19 using the cantilever moment distribution method. 



Relative Stiffness Values and Distribution Factors 

/Car =Xff = -= 0.50/ 



^BC ~ ^DE ~ T ~ 0.33/ 



^BE _ ^CD " 



3 

6(20 



= 3.0/ 



(DF)be = 



0.50 



0.50 + 0.33^3.0 
3.0 

= 0.783 



= 0.130 



3.83 



(DF)ec = ^§ = 0.087 

(DF) CB =-^— = 0A0 

0.33+3.0 



ma, = W3 - 0.90 



4KN 



lOkNB 



I 



21 



21 



21 



4.0m 



21 



Figure 4.19 



o 



- 



E 

O 
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Fixed-End Moments (Relative Values) 

The fixed-end moments are taken as the product of the storey-shear and the 
storey-height since the columns behave as a cantilever. Therefore, 

MSc = M F 

Subsequently, in the upper storey the fixed-end moments are found by 
multiplying the storey-shear (4.0 kN) by the column height (3.0 m). In the lower 
storey the total horizontal force is 14.0 kN. 



Thus 



Mbc = Mg B =M£ D =Mg E = 4x3 = 12kNm 
Af£ B = Afg A =M£ F =Mf E = 14 x 4 = 56 kN m 



Table 4. 13 



Joint 



B 



Member 



A B 



BA 



BH 



0.50 



0.50 



BC 



3.0 



0.33 



CB 



0.33 



CD 



3.0 



DF 



0.0 



0.130 



0.783 



0.087 



0.10 



0.90 



Fixed-end 
moment 



Total 



+56.00 
+8.84 



+56.00 
-8.84 



-53.24 



+ 12.00 
-5.92 



+0.23 



-0.23 



-1.40 



+ 1.79 
-0.16 



+ 12.00 
+5.92 

-1.79 



+65.07 



-0.00 



+46.93 



-0.02 



-54.66 



+0.02 
-0.00 

+7.73 



+0.16 
-0.02 



+ 16.27 



-16.13 



-16.27 



Correction Factor 



Upper storey shear 



1 



Vy =j(M CB +M B c+M DE +M ED ) 



= - (7.73 + 16.27)= 16 kN 
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Lower storey shear 

F 2 = | (65.07 + 46.93)= 56 kN 

A correction factor must be applied at each level to satisfy horizontal 
equilibrium. 

Since the relative fixed-end moments are taken as the product of storey-shear 
and storey-height, the correction factor becomes common to both storeys and is 
given as 

4 14 

correction factor = — = — = 0.25 

16 56 

The final end moments are obtained by multiplying those obtained in Table 4.13 
with the correction factor. 



Final End Moments 





= 0.25(+65.07) ' 


= +16.27 kN m 




= 0.25(+46.93) ' 


= +11.73 kNm 


Mbe 


= 0.25(-54.66) 


= -13.67 kNm 


Mbc 


= 0.25(+7.73) = 


+ 1.93 kN m 


Mcb 


= 0.25(+16.27) 


= +4.07 kN m 




= 0.25(-16.27) 


= -4.07 kN m 



4.8 ARBITRARY LOADING ON SYMMETRIC FRAMES 

Any arbitrary loading on a symmetric frame can be resolved into symmetric and 
antisymmetric loading system components. The frame may then be analysed by 
using the conventional moment distribution method to the symmetric loading 
system, and the cantilever moment distribution method for the frame subjected 
to the antisymmetric loading. The final end moments are then obtained by 
adding algebraically the results of the two solutions. Figure 4.20 shows a 
symmetric frame subjected to any arbitrary loading. The same frame is also 
shown loaded by symmetric and antisymmetric (Fig. 4.20(b) and (c)) loading 
systems whose algebraic sum furnishes an equivalent system to the original 
loading. 
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P/2 JPA 



= B 



W/2 



I I I I I I I 



D C 



1 P/2 



E + B 



F A 



P/2 



W/2 



i_±l 



-W/2 



4% 



(a) Arbitrary loading 



(b) Symmetric loading 



(c) Antisymmetric loading 



Figure 4.20 




Figure 4.21 
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EXAMPLE 4.8 Determine the joint moments of the two-storey frame in 
Fig. 4.21. 

(a) Symmetric Loading 

The frame is analysed with the conventional moment distribution method on 
half of the frame only, since the frame and the loading are both symmetrical. 



Relative Stiffness and Distribution Factor 



KAB=K E F=f = 0.50/ 
Kbc = Kde={ = 0.33/ 

Kbe =K cd =^=0.50/ 

(DF) R a = (DF) nv = = 0.375 

K ;ba v ;be 0.50 + 0.33 + 0.50 

0.333 

0.333 n „ 

- o^TKTo - a40 
WcD = S = °' 60 



Fixed-End Moments 

mF 5(2.5)(1.5) 2 5(1.5)(2.5) 2 

WD = -^DC = (4 Q) 2 + (4 0) 2 " 4 ' 69 kN m 

Mg E = -Mg D =!^£= 1.33 kN m 

Table 4.14 shows the distribution of the fixed-end moments carried over the 
left half of the frame. The moments carried from the right half are indicated in 

italics. 
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Table 4. 14 



Joint 


A 


B 


C 


Member 


AB 


BA 


BE 


BC 


CB 


CD 


K 


0.50 


0.50 


0.50 


0.333 


0.333 


0.50 


DF 


0.0 


0.375 


0.375 


0.250 


0.40 


0.60 


Fixed-end 
moment 


-0.25 


-0 50 


+ 1.33 
-0.50 


-0.33 


-0.17 


+4.69 








+0.25 


-0.90 


-1.81 


-2.71 




+0.12 


+0.24 


+0.24 


+0.17 


+0.08 


+1.35 








-0.12 


-0.29 


-0.57 


-0.86 




+0.07 


+0.14 


+0.14 


+0.09 


+0.04 


+0.43 








-0.03 


-0.09 


-0.19 


-0.28 




+0.03 


+0.06 


+0.06 


+0.05 


+0.02 


+0.14 








-0.03 


-0.03 


-0.06 


-0.08 




+0.01 


+0.02 


+0.02 


+0.02 






Total 


-0.02 


-0.04 


+ 1.32 


-1.27 


-2.66 


-2.66 



(Z>y Antisymmetric Loading 

The cantilever moment distribution method is used to analyse the frame, 



Relative Stiffnesses and Distribution Factors 



21 

Kab = J = 0.50/ 



K BC = - = 0.33/ 

^BE = ^CD = = 3-0/ 



0.50 



0.50 + 0.33 + 3.0 



= 0.130 



(DF) BC = m| = 0.087 
BL 0.83 
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Wbe= |^ = 0.783 

0.33 

(DF)cr = = 0.10 

K ;cB 0.33 + 3.0 

3.0 

£DF)bc = — = 0.90 
Fixed-End Moments 

The fixed-end moment for a beam with a uniformly distributed load extended 
kL from the left end (support A) is 



Aff n = wL 2 



AB - . 



j (6 - 8* + 3* 2 ) 



Mg A =-wL 2 l^)(4-3*) 
For the problem at hand, k = 0.5. Thus 

Jtft = 1.0(4.0) 2 [6 _ 8( o.5) + 3 (0.5) 2 ] 

-1.0(4.0) 2 ^(4-3(0.5)) 

= +0.500 kN m 

F _ 5(1.5)(2.5) 2 5(2.5)(1.5) 2 
^ D -- (4.0) 2 + (4.0) 2 

= -1.172 kNm 

These are distributed using the cantilever moment distribution method as shown 
in Table 4.15. 



Final End Moments 

M AB = -0.02 + 0.050 = +0.030 kN m M DC = +2.66 - 0. 1 56 = -2.504 kN m 

M BA = -0.04 - 0.50 = -0.090 kN m M DE = +2.66 + 0. 146 = +2.806 kN m 

M BC = -1.27 - 0.146 = -1.416 kN m M EB = -1.32 + 0.196 = -1.124 kN m 

M BE =+1.32 + 0.196 = +1.516 kNm M ED = +1.27 - 0.146 = +1.124 kN m 
Mqb = -2.66 + 0.146 = -2.514 kN m = +0.04 - 0.050 = -0.10 kN m 

Mc D = -2.66 - 0.156 = -2.816 kN m M FE = +0.02 + 0.050 = +0.070 kN m 
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Table 4.15 



Joint 


A 


B 


C 


Member 


AB 


BA 


BE 


BC 


CB 


CD 


K 


0.50 


0.50 


3.0 


0.33 


0.33 


3.0 


DF 


0.0 


0.130 


0.783 


0.087 


0.10 


0.90 




+0.065 
-0.015 


-0.065 


+0.500 
-0.392 


-0.043 


+0.043 
+0.1 13 


-1.172 
+ 1.016 


+0.015 


+0.088 


-0.113 
+0.010 


-0.10 




Total 


+0.050 


-0.050 


+0.196 


-0.146 


+0.146 


-0.156 



4.9 PROBLEMS 



4.1 Find the support moments of the continuous beam shown in Fig. P4.1. 



{Am: M A = -12.6 kNm 
M B = -14.4 kNm 

M c = -2.2 kN m) 




2kN/m 

1 1 1 1 1 1 1 1 n m 




8kN 
4kN/m 



TT 



21 



B 




1 



lOkN 



6m 



-HH- 



3m _j 


2m | 


6m 









Figure P4.1 



4.2 Determine the joint moments of the frame shown in Fig. P3.2. 

4.3 Draw the bending moment of the frame shown in Fig. P4.3. 
(4ns: M A = -14.72 kNm 

M B =+10.19 kNm 
Mq = -6.73 kN m) 
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lOkN 



E 
10 



B 



A '///////. 



21 



I 



31 




6m 



Figure P4.3 



4.4 Calculate the support moments of the frame of Fig. P4.4. 



'///////. £ 



lOkN 



B 



21 



A 7/)////,. 




E 



6m 



I m 4m 




Figure P4.4 

4.5 Find the joint moments of the split-frame shown in Fig. P4.5 

(Ans: M A = -7.93kNm 
M B = +6.08 kN m 

Mc = -6.65; -1.69; -4.96 kN m 
M D = +7.37 kN m 
M E = +3.49 kN m 



Afp = -2.48 kN m) 



3kN E 



E 

CM 



4kN 




21 




8m 



Figure P4.5 
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4.6 Find the joint moments of the two-storey frame shown in Fig. P4.6 using 
the cantilever moment distribution method. 



!5kN 




3.5 m 



£ 

CM 




E 



Figure P4.6 

4.7 Find the joint moments of the two-storey frame shown in Fig. P4.7 using 
the cantilever moment distribution method. 



04 aw: M ab = +0.63 kN m 



M D = -5.35 kN m 



M BA = +0.19 kN m M ED = +1.69 kN m 
M BE = -3.25 kN m M EB = +2.46 kN m 
M BC = +3.44 kN m A*e F = -0.78 kN m 



M c = -3.6 kN m 



M FE = -0.34 kN m) 




Figure P4.7 
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5. Kani Method of Moment 
Distribution 



5.1 INTRODUCTION 

The moment distribution by the Cross and Kani methods are both iterative 
procedures to solve the slope deflection equations. However, the Cross method 
obtains the unknown end moments by iterating the moment increments, while 
Kani's method iterates the end moments themselves as the unknowns. Kani's 
method consists of carrying out a single operation applied repeatedly at the 
joints of a structure in an arbitrary sequence. Using Kani's method the results 
may also be obtained to any desired accuracy by continuing the calculations a 
sufficient number of times. In addition to its simplicity, the method has the 
advantage of acquiring a built-in error elimination scheme. Moreover, the method 
is more suitable for frames with high degree of redundancy, including frames, 
with sidesway, since only one set of computations is necessary. For such frames 
the required computational effort by Kani method is minimal when compared to 
other methods. 



5.2 FRAMES WITHOUT SIDESWAY 

Consider member f-m as integral part of a frame (Fig. 5.1) and that there are 
many such members meeting at joint / so that m is the general designation, for 
the far ends. Let M jm and M mj - be the end moments due to the applied loads at 
joints / and m , respectively. 

The slope deflection equation for member j-m at joint / which is considered 
part of a frame without sidesway is written as 



M im =Mf m -2EK im (26 i + d m ) 



[5.1] 



The same equation may be written in the form 
Mjm = Mf m + 2M' jm + M' mj 



[5.2] 
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(a) Loading system 



M 



in 




■ 

J 



m 



(b) Fixed-end moments 



M 




m 




1 




(c) Moment at / 





(d) Moment atm 

Figure 5.1 



where 



M' jm = -2EK jm dj 
M' mj = -2EK jm d m 



ML ■»£ 



m 



Since M ; >„ constitutes the contribution by 0, to the total moment M jm it is 
referred to as the rotation contribution of the end m to M ]m . 

For any joint / where the number of m members are joined to it, the joint 
being in equilibrium of end-moments gives 



F/m = 0 



or 



2^ + 2(2^+^ = 0 



[5.3] 



Defining M ; the algebraic sum of the fixed-end moments at joint / as the 
restraint moment 



Mj = 2 Mf m 
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Then [5.3] may be written as 

2M; =-\(Mf+XM' mf ) [5.4] 



m z m 



Since M jm for any member must be proportional to the relative stiffness of the 
member, the moment in any member j-m is 

M'jm = ^jr- 3>n 
LK jm m 

Denoting the term -\{K jm I^Kj m ) as the rotation factor R^, [5.5] may be 
written as 

M' jm =R jm (M j + ZM' mj ) [5.6] 



Equation [5.6] forms the basis for Kani's method for frames without sidesway 
where the rotation contributions are evaluated for member /-m. These 
contributions and the fixed-end moments are then added algebraically to 
determineMy m (see [5.1]). ButM^ must initially be evaluated. However, their 
final values may be determined by successive approximations using the Gauss— 
Seidel iteration scheme where the rotation contributions proceed from estimated 
values (starting with initial zero-approximation), the subsequent values being 
obtained with better approximation. The iteration is terminated when the latest 
approximation furnishes a value acceptably close to the preceding result. 
Kani's method for frames without sidesway may be summarised as follows: 

(a) Determincthe fixed-end moments of all members. At each joint evaluate 
the resultant fixed-end moment, restraint moment, XMj. 

(b) Calculate the relative stiffnesses (K values) of all members using the Gauss- 
Seidel iteration scheme. 

(d) Evaluate the final end-moments using [5.2] . 



Hinged-End Members 

In a structure which contains a hinge located at one end, the stiffness of the 
member becomes three-fourths of the stiffness of a corresponding beam fixed at 
both ends. Such members, after the fixed-end moments are determined based on 
the actual member length and introduced into the calculation scheme, are 
replaced by fictitious members fixed at the hinged ends with the K values which 
are three-fourths of the actuals values of the members in the original structures. 
This substitution is justified since the end moment required to produce a unit 
rotation in the original member is the same as the substitute with three-fourths 
stiffness. The hinged-end is then left free and its final end-moment is set to zero. 
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EXAMPLE 5. 1 Find the joint end moments of the symmetric frame shown in 
Fig. 5.2. 

Fixed-End Moments 

Mg E =Mef = --Wed = -M& = ^ff = 9.0 kN m 

Mq H = A/hi = -Mhg = —Mm = 7T = 6.0 kN m 

The fixed-end moments are recorded at the corresponding member ends in the 
scheme of calculation (Fig. 5.3). 



2.0kN/m 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


I H 

1.51 

3.0kN/m 


I 

1.51 1.51 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 


iii iii ii 1 1 1 : i 


21 E 

21 

B 


21 

21 21 

C 



E 
O 

rn 




E 

o 




6.0m 



6.0m 



Figure 5.2 



Restraint Moments (M, = TMfa) 

m 

Mq = +6.0 kN m 
M H = (-6.0 + 6.0) = 0.0 
Mi = -6.0 kN m 

M D = +9.0 kN m 

M E = 0.0 

M F = -9.0 kN m 
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m 




r i 

H 


o 




O 






f6.0 






o 


0.0 


o 


+fi no 




d 


-0.75 


+ 0.08 


d 


d 


+0.08 




0.375 


-0.60 


-0.02 




0.300 


-0.02 




0.00 




-0.58 


-0.00 




0.00 




0.00 




2.25 




+ 0.57 


0.00 


+ 


0.23 




0.00 




1.80 








+ 


0.07 






1 .73 










0.0/ 








1 .74 










0.00 







/r 



-600 ^ 



+ 0.74 
+ 0.60 
+ 0.58 
+0.57 




0.00 
+ 2.22 
+ / .78 
+ / .74 
+ /.74 



-/.36 
-/.37 
- /.36 
-1.27 
0.00 





0.00 
0.00 
0.02 
0.09 
0.00 



0.00 
0.00 



0.00 
+-0.09 

+ 0.02 
0.00 
0.00 



0.00 
0.00 



+ / .36 
+ / .36 
+ /.35 
+ f.26 
0.00 



-0/875 




c 




0./50 




r 


•0.1675 




+ 9.0 




+ 9.00 


-9.00 




0.0 


O 


+ 9.00 


-0.90 


m 
d 


-9.0 






d 


- 0.84 


+ 0.06 




d 


+ 0.06 


+ 0.84 




-0./875 


-0.9/ 


+ 0.0/ 


-0./50 


i 


+ 0.0/ 


+ 0.90 


-0./875 



+ 0.9/ 
+ 0.9/ 



0.00 
+ / .26 

+ i .35 
+ / .36 
+ / .36 



6 

Figure 5.3 



Relative Stiffness Values (K^ = I^/L^) 



&AD ~ ^BE ~ ^CF ~ T ~ 0-50 



^DG ~ ^EH " ^FI " ~ " 0-50 



^DE ~ ^EF ~ T ~ 0-333 

o 



*gh=*hi=T = 0.167 
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Rotation Factors (Rp„ = -Q.5K jm l'LK jm ) 

0.50(0.5) 

^ DA = " 0.5 + 0.333 + 0.50 = "°- 1875 

0.5(0.333) 

*«— U33 =-° 125 

Check. ItRj m = -0.1875 - 0.1875 - 0.125 = -0.50 
Similarly 



^GD 




-0.375 


^GH = 


-0.125 


^HG 




-0.10 


R HE = 


-0.30 


^IH 




-0.125 


^IH = 


-0.375 


^ED 




-0.10 


^EH = 


-0.15 


^EB 




-0.15 


^EF = 


-0.10 


^FE 




-0.125 




-0.1875 



R m = -0.10 



The rotation factors are recorded at the corresponding ends of the members in 
the computational scheme (Fig. 5.3). 



Rotation Contributions 

The contribution to the end moment M jm by the rotation 0 y at joint / is given by 

M; m =R Jm (M j + I,M' mi ) 

The calculation of the rotation contribution may be started at any joint and 
continued at other joints in any chosen sequence. The sequence adopted here is 
GHI-DEF. 

First Cycle 

(a) Joint G Since the joints H and D are initially locked, the contributions to 
the end moments at G are zero. Thus, the initial values of these joint moments 
are set to zero or M^ G = M& G = 0 and the rotation contributions at G are 

M GH = -0.125 (+6.0 + 0.0 + 0.0) = -0.75 kN m 

M GD = -0.375 (+6.0 + 0.0 + 0.0) = -2.25 kN m 
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These rotation contributions are entered at the joint G below the respective 
fixed-end moment as shown in Fig. 5.3. 

(b) Joint H At this joint again, M E h = M{ H = 0.0. But, Mq H = —0.75, as 
computed above at joint G. Therefore 

Mh G = -0.10(0.0 - 0.75 + 0.0 + 0.0) = +0.08 kN m 
Similarly 

MAe = -0.3(0.0 - 0.75 + 0.0 + 0.0) = +0.23 kN m 

Mhi = -0.1 (-0.75) = +0.08 kN m 

In the same manner the calculations are performed at the other joints until the 
first cycle is completed. 

Second Cycle 

The results obtained from the first cycle are used to obtain better approximations 
to the rotation contributions. For example, consider the rotation contribution at 
joint D: 

Mdg = -0. 1 87(+9.0 - 1 .80 + 0.06) = - 1 .36 kN m 

A/de =-0.125(9.0- 1.80 + 0.06) = -0.91 kN m 

Mn A =-0.188(9.0- 1.80 + 0.06) = -1.36 kN m 

A similar calculation is performed at all joints until the second cycle is 
completed. This procedure is performed until the two successive cycles furnish 
values differing by only a small acceptable amount. In this particular example, 
four cycle operations were considered adequate and the complete calculation 
scheme is shown in Fig. 5.3. 

Final End Moments 

The final end moments are determined from the relation 

Substituting the fixed-end moment values and the rotation contributions 
obtained from the scheme of calculation: 

M AD = 0.0 + 2(0.0) - 1 .36 = - 1 .36 kN m 

Af DA = 0.0 + 2(- 1 .36) + 0.0 = -2.72 kN m 

M DE =+9.0 + 2(-0.91) + 0.0 = +7.18 kN m 

M DG = 0.0 + 2(- 1 .36) - 1 .74 = -A .46 kN m 



131 



METHODS OF STRUCTURAL ANALYSIS 

M ED = -9.0 + 2(0.0) - 0.91 = -9.91 kN m 
M GD = 0.0 + 2(-1.74) - 1.36 = -4.84 kN m 

M GH = +6.0 + 2(-0.58) + 0.0 = +4.84 kN m 
M HG =-6 + 2(0.0) - 0.58 = -6.58 kN m 

^EH = ^HE = ^EB = ^BE = 0 

Notice that the final end moments to the right of the structural axis are equal in 
magnitude but opposite in direction due to the symmetry of both the loading 
system and the structure itself. 



5.3 FRAMES WITH SIDESWAY 



When a frame is either structurally unsymmetric or is symmetric with 
unsymmetrical loading, joint translation or sidesway occurs. Figure 5.4 shows 
member j—m of a frame with lateral displacement. The rotations at joints / and 
m are 6^ and 0 m/ - respectively, and A /m is the relative lateral displacement 
between; andm. 

The slope deflection equation for the member j-m is 

= Mfr - lEKj^lOj + 9 m - ILpnlLfr) [5.7] 

or 

M jm = Mf m + 2M; m + M' mi + Ml' m [5.8] 
where 

Mfm = -lEK^dj 

M'mj = -ZEKjmQm 

The symbols Mj m andM' mj ; respectively, define the rotation contributions of 

the joints/ and m to the total moment M jm . In a similar manner, M'j m constitutes 

the contribution of M jm by the displacement A /m and is therefore defined as the 

displacement contribution. 

The algebraic sum of the end moments of all members meeting at joint / is 
zero. 



XM fm = 0 



or 



XMf m +-L(2Mj m +M' mj +M^) = 0 [5.9] 



m m 
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M 




jm 



} 



f rrftTI 



■yVn 



(a) Original system 



(b) Fixed end moments 




(c) Moment at/ 




(d) Moment atm 



m 





M 



M 




(e) Relative displacement 
Figure 5.4 




Defining the algebraic sum of the fixed-end moments at joint / as the restraint 
moment, 

then [5.9] may be written as 



[5.10] 



m 



Equation [5.10] forms the basis for Kani's method for frames with sidesway 
where the rotation and displacement contributions must be evaluated for 
member f-m. These contributions and the fixed-end moments are added 
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algebraically to determine the final joint moment M jm (see [5.8]). However, 
while the evaluation of the rotation contributions are easily computed as 
described in Section 5.2, the determination of the displacement contribution 
Mfa could be involved depending on the type of loading on the structure and 
whether there are columns of different length. The cases arising in the determina 
tion of displacement contributions are discussed below. 

5.3.1 Vertical Loading 

Consider the frame shown in Fig. 5.5(a) subjected to vertical loading only. The 
frame undergoes sidesway due to the unsymmetrically placed vertical loads. The 
frame analysis may be carried out in two steps as was done in the Cross moment 
distribution (Section 4.6). 

(a) No-sway Solution 

As shown in Fig. 5.5(b), artificial joint restraints are applied at storey heights to 
prevent sidesway. The analysis of the frame without sidesway follows the same 
procedure presented in Section 5.3. After the joint moments are evaluated the 
artificial joint restraints may be determined if required, from equilibrium 
considerations of the shear forces at every storey. 

fb) Sway Solution 

Since the artificial joint restraints introduced in the no-sway solution do not 
actually exist in the given frame, their presence may be nullified by applying a 
consistent force system (Fig. 5.5(c)) whose forces are equal in magnitude but 
opposite in direction to the respective artificial joint restraints. 

By cutting horizontally through all columns at the rth storey and from the 
consideration of equilibrium conditions, the algebraic sum of the column shear 
must be zero, 



where V jm = shear in column j-m of the rth storey. Let h r represent the column 
height of the rth storey, then 




0 




[5.11] 



But 




M mj = + 2M' mj + M; m + Mj' m 
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Since vertical loading only is considered, there are no intermediate loads on the 
column j-m and the fixed-end moments are zero. Thus 



M^=M'^ = 6EK ' mAim [5.12] 
Substituting into [5.11], 

Vjm = I Wjm + 3M' mi + 2l£) [5.13] 
Taking the sum of the column shears at the rth storey, 

ZV /m = 2 f [W'jm +M' m j) + = 0 [5.14] 

Solving for the sum of the displacement contribution, 

pfc--f septic*) 

Hence, the algebraic sum of the displacement contribution of all columns at the 
rth storey is determined to be —1.5 times the sum of the rotation contributions 
of the column ends of the same storey. Since all the columns at this storey 
undergo the same lateral displacement, the displacement contribution Mj' m for 
any individual column is proportional to its relative stiffness value (see [5.12]). 
The displacement contribution of any column is, therefore, obtained by 
distributing the sum I.Mj' m among the columns in the rth storey in proportion 

to their stiffness values. Thus, for any column j-m 

[^)^(Mpn+M' mJ ) [5.15] 

where ZK jm = the sum of the K values of all columns of the rth storey. 

In order to make the computation more convenient, and to make an analogy 
to rotation factor, a displacement factor D jm is defined for column j-m as 




[5.16] 



Equation [5.15] may be written as 

M;' m =D jm X(M; m +M' mj ) [5.17] 

Since the displacement contribution Mj^ and the rotation distribution M' jm and 
M' mj are interrelated, [5.17] may be solved by using the same convenient Gauss- 
Seidel iteration scheme. Once the displacement contributions are known the 
rotation contributions Mj m are obtained from [5.6] and [5.9] as 

M fm =R im \M / + 2<M' mi +M i H m )] [5.18] 
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and from these, again, the following approximation of the displacement and 
rotation contributions are calculated until the results of a desired accuracy are 
obtained. 

The final end moments are then obtained by using [5.8]. 



EXAMPLE 5.2 Find the joint moments of the symmetric frame subjected to 
unsymmetric vertical loading shown in Fig. 5. 6. 
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Figure 5.6 



Fixed-End Moments 



M F - M F -WO.- 
m EF — JVI FE j 2 



= 9.0 kNm 



o 



e 
o 



Restraint Moments (Mj = *LMf m ) 

Tfl 

M G =M D =0 

M H = -My = 6.0 kN m 

Me = -M F = 9.0 kN m 

The relative stiffness values and the rotation factors have the same values as in 
Example 5.1. 
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Displacement Factors (D jm = -1.5K Jm /2K jm ) 

Aud=Aje=-Dcf = -1-5 (y|) = - 0 - 5 

Z>DG=Am=/>Fi = -1.5 (t| I = -°- 5 

Notice that the restraint moments Mj and the rotation factors R jm are 
recorded in the usual manner (Example 5.1) in the computational scheme shown 
in Fig. 5.7, The displacement factors are also entered at the centre of the 
relevant columns in the same computational scheme. 

First Cycle 

The rotation contributions are first computed which are then used to evaluate 
the displacement contributions. 

(a) Rotation Contributions, Mj m = R im [Mj + Z(M' mJ + Mj' m )] 

m 

Since M' mj and M-' m are initially known, they are set to zero at the first cycle. 



(i) Joint G: 



SetA^ G =MA G =MGD = 0.0 



M G = 0.0 

Thus 

M GD = A/ GH = -0.375(M G + + A& G + Mfc) = 0.0 
(ii) Joint H: 

Set M^ H = ^fh = Meh = 0.0 

M GH = 0.0 (as found in (i)) 
M H = +6.00 

Thus 

Mhg = -0. 100(6.0 + 0.0 + 0.0 + 0.0) = -0.60 kN m 
Similarly 

Mhe = -0.300(6.0 + 0.0 + 0.0 + 0.0) = - 1 .80 kN m 

Mhi = -0.100(6.0 + 0.0 + 0.0 + 0.0) = -0.60 kN m 

In the same manner all rotation contributions are entered below the relevant 
fixed-end moments. 
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Figure 5.7 



ffij Displacement Contributions, Mj' m = Dj m 2(M} m +M' m j) 

The displacement contributions are computed by multiplying the algebraic term 
of the rotation contributions of all columns of a storey with the displacement 
factors of the individual columns. 
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Thus 

Mug =Me H +M& = -0.50(0.0 + 0.0 - 1.80 - 1.08 + 2.48 + 1.36) 
= -0.48 kN m 

Similarly 

Ma D = Mbe = M CF = -0.50(0.0 + 0.0 - 1 .08 + 0.0 + 1 .36 + 0.0) 
= -0.14kNm 

The displacement contributions are recorded at the middle of the relevant 
columns as shown in the scheme of calculation (Fig. 5.7). 

Second Cycle 

The results obtained from the first cycle are used to obtain better approximations 
to the rotation and displacement contributions. 

(a) Rotation Contributions 

(i) Joint G: 

M GD = -0.375(M G +Md G + M HG +^dg) 

= -0.375(0.00 + 0.00 - 0.60 - 0.48) = 0.41 kN m 

Similarly 

M GH = -0.125(0.00 + 0.00 - 0.60 - 0.48) = 0.14 kN m 
(ii) Joint H: 

AfcH = +0. 1 4 kN m (as in (i) above) 
Meh =-1.08 kNm 

Mm = +0.83 kN m 
Meh = -0-48 kN m 

M H = +6.0 kN m 

Therefore 

M HG =MAi = -0.100(6.0 + 0.14- 1.08 + 0.90-0.81) 
= -0.54 kN m 

M HE = -0.300(6.0 + 0.18 - 1.08 + 0.83 - 0.48) 
= -1.62 kNm 
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Similarly, all rotation contributions of the second cycle are entered below the 
relevant fixed-end moments. 

(b) Displacement Contributions 

M^g =M"eh =Mh = -0.50(0.41 +0.17 - 1.62- 1.17 + 2.12+ 1.56) 
= -0.74 kN m 

Kd =A*be =M C¥ = -0.50(0.17 + 0.00 - 1.17 + 0.00 + 1.56 + 0.00) 
= -0.28 kN m 

This completes the second cycle operations. The procedure is then repeated until 
two successive cycles furnish sets of values differing by a very small acceptable 
amount. In this particular example, the scheme of computation shows four cycles 
to be sufficient (Fig. 5.7). 

Final End-Moments, M jm = M f m + 2M? m + M' mj + Mj' m 

M AD = 0.0 + 2(0.0) + 0.33 - 0.45 = -0.12 kN m 
M DA = 0.0 + 2(0.33) + 0.0 - 0.45 = 0.21 kN m 
M DE = 0.0 + 2(0.22) + (-0.75) = 0.3 1 kN m 
M FE = -9.0 + 2(1 .12) + (-0.75) = -7.5 1 kN m 
and so on. 



5.3.2 Horizontal Loading 

The principle in determining the displacement contributions in frames subjected 
to horizontal loading remain the same as in vertical loading. However, the 
presence of horizontal loads on the frame requires an additional effort in the 
computation. 

Consider a frame subjected to horizontal loads applied as shown in Fig. 5.8(a). 
Again, the analysis may be carried out in two steps: 



(a) No-Sway Solution 

Artificial joint restraints are applied at storey heights as shown in Fig. 5.8(b) to 
prevent sidesway. These joint restraints may be determined from no-sway 
solution (section 5.3). 
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(b) Sway Solution 



Since the artificial joint restraints do not actually exist, they may be eliminated 
by applying a consistent force system (Fig. 5.8(c)) whose storey shear at any 
section is equal in magnitude but opposite in direction to the algebraic sum of 
the applied horizontal force above that section. Designating the sum of the 
restraint forces above the rth storey as the storey shear V r , the horizontal 
equilibrium condition above the rth storey requires that 

V r =H x +H 2 +...H r = ZH 

= X-(M Jm +M mf ) [5.19] 

Assuming the columns are not subjected to intermediate horizontal loads 
(such as by applying equivalent loads at the storey heights that will give the same 
global effect on the frame), all fixed-end moments become zero. Using [5.8] to 
expand the end-moments M jm andM m; -, [5.19] may be written as 

Rearranging: 

™i'm = - f [ - ^ + + Mm/)] [5-20] 

The quantity V r h r /3, which is one-third of the product of the storey shear and 
storey height, is defined as the storey moment M' r . 

3 

Equation [5.20] gives the sum of the displacement contributions of all columns 
in the rth storey. As explained in the case of vertical loads, the moment in any 

column /-m is obtained by distributing this sum in proportion to their K values. 
Thus 

MU =D im [M r +X(Mj m +M' m ,)] [5.21] 

Notice that the displacement factors are the same as for vertical loading, so 
that the displacement contribution in the case of horizontal loading (see [5.21 ]) 
differs in the extra term M r from the case of vertical loading (see [5.17]). 
Therefore, the analysis of frames subjected to horizontal loading differs from the 
analysis of frames with vertical loading only by the extra term M r which must be 
calculated for each storey and be added algebraically to the sum of the rotation 
contributions of the two ends of the columns of the storey considered. 
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EXAMPLE 5 .3 Find the joint moments of the frame subjected to horizontal 
loads as shown in Fig. 5. 9. 

The relative stiffness values, rotation factors, and displacement factors are the 
same as in Example 5.2 and are recorded in the computational schemes as usual 

(Fig. 5.9). 
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Storey Shears and Storey Moments 
(i) Second storey 



Storey shear, V r = 3.0 kN 

Storey moment, M r = - ^ = " 3, °^ 3 '°^ = -3.0 kN m 



(ii) First storey 

Storey shear, V r = 3.0 + 8.0 = 1 1 .0 kN 

11.0(4.0) 1JlMiM 
Storey moment, M r = = -14.67 kN m 

The storey moments are recorded in the computational scheme (Fig. 5.10) at 
the centre of the relevant storey. 
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Figure 5.10 
Displacement and Rotation Contributions 

In frames subjected to horizontal loads, the displacement contributions are 
usually significantly larger than the rotation contributions. Hence, the displace- 
ment contributions are calculated first as illustrated in the following calculations. 
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First Cycle 

(a) Displacement Contribution, Mj' m = D jm [M r + 2 (Mj m + M' mj )] 

Since the rotation contributions are initially zero, the displacement contributions 
are 

AfrjG =^eh -Jfi = -0.5(-3.0 + 0.0 + 0.0) 

= +1.50 kN m 
M A ' D =Mb' e +M:f = -0.5(- 14.67 + 0.0 + 0.0) 

= +7.33kNm 

(b) Rotation Contributions, Mj m = R jm [M' mj + S (M' mj + Mj' m )] 

TTX 

(i) Joint G: 

At this joint M G = M^ G =Mq G = 0 

andA/cD = -1.50kNm 
Mod = -0.375(0.0 + 0.0+1 .50) = -0.56 kN m 
M GH =-0.375(0.0 + 0.0+ 1.50) = -0.19 kN m 

(ii) Joint H: 

At this joint, M H = M^ H = M{ H = 0 

M GH =-0.19kNm 
and A/eh =+1.50 kN m 

Thus 

M HG =-0.100(0.0-0.19+0.0 + 0.0+ 1.50) = -0.1 3 kN m 
Similarly 

A/he = -0.300(0.0 -0.19 + 0.0 + 0.0+ 1 .50) = -0.39 kN m 
A/A, = -0.100(0.0 - 0.19 + 0.0 + 0.0+1 .50) = -1.3 kN m 

In the same manner all the rotation contributions are calculated until the 
first cycle is completed. 
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Second Cycle 

(a) Displacement Contributions 

The displacement contributions are obtained by using [5.21] where the results 
from the first cycle are used to obtain approximations. 

Mvg 'Mm =Mn = -0.5(-3.0 - 0.50 - 1.55 -0.39 - 1.11 -0.51 - 1.43) 
= +4.28 kN m 

Mad =Mbe =Mcf = -0.50(-14.67 - 1.55 - 1.11 - 1.43) 
= +9.38 kNm 

(b) Rotation Contributions 

There are no particular points to be noted here and similar calculations are 
performed until the second cycle is completed. The procedure is repeated until 
two successive cycles furnish sets of values differing by a very small acceptable 
amount. In this particular example, the scheme of computation shows seven 
cycles to be sufficient (Fig. 5.10). 

Final End Moments, M jm = Mf m + 2M' jm + M' mj + Mj m 

M AD = 0.0 + 2(0.0) - 3.03 - 1 1.35 = +8.32 kN m 
M DA = 0.0 + 2(-3.04) + 0.0 - 1 1 .35 = +5.29 kN m 
M DE = 0.0 + 2(-2.02) - 1.37 = -5.41 kN m 
M DG = 0.0 + 2(-3.04) - 1.62 + 7.78 = +0.10 kN m 
Mod =0.0 + 2(- 1.62) -3.03 + 7.78 = +1.51 kNm 
M GH = 0.0 + 2(-0.54) - 0.47 = - 1 .55 kN m 
and so on. 

5.3 .3 Frames With Columns of Unequal Heights 

For a frame with a storey containing columns of unequal heights, the calculations 
of the rotation contributions in all storeys remain the same as described earlier. 
Also, the computation of the displacement contributions for those storeys with 
equal heights are not altered. However, in establishing the governing equations 
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for the calculation of the displacement contribution for the storey with unequal 
heights, supplementary consideration is needed. 

Consider the frame shown in Fig. 5.1 1 . In the storey which has columns of 
unequal heights, an arbitrary column that appears most frequently is taken as 
the storey height. Let 

h r = storey height in the rth storey which has columns of unequal 
heights 

hj m = height of any other column j-m in the rth storey 





T 



Figure 5.11 




jm 



Writing the equilibrium condition of the horizontal forces at the rth storey 



V r + 2V Jm =Q 



[5.22] 



or 



1 



V r + X — (M jm + M mi ) 



V r h r + 2- I -(M im +M mj ) = Q 



r h 
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Introducing a factor defined as the height reduction factor for the f-m 
column in the rth storey, 

c 

hjm 

Equation [5.22] is written as 



[5.23] 



Qrh r + XC jm (M jm +M mj ) = 0 [5.24] 
Substituting the values of andM m/ - given by 

M jm = Mf m + + M' mj + M- m 
M m , = Af£, + 2M' mj + M) m + Mj' m 

into the shear equations and taking into consideration that the fixed-end 
moments are zero, 

VJi r + ZC jm (3M' jm + 3M' mj + 2M'/ m ) = 0 

Therefore 

XC im M;' m = -1.5(M, + ZC jm {M' jm +M' mj ) [5.25] 

where 



VJi 

M r = - J — L = storey moment 

3 



Since 



6 EK jm &jm 
M im— h 



jm 

them Mj' m is proportional to Kj m /h jm and also to C jm K jm . Also, since A jm is the 
same for all columns of the storey under consideration, 

-LC jm M jm = XC? m K im t 5 - 26 ! 

From [5.25] and [5.26], the basic equation for determining the displacement 
contribution M'^ n may be written as 

M'; m =D' jm [M r + SC) m (M^ +M' m/ )] [5.27] 

where 

v ^ 2 i = displacement factor 

fC jm K jm I 

For the storey with unequal column heights, the following changes must be 
noted : 
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(a) a reduction factor C jm must be introduced; 

(b) a modified displacement factor Dj m must be used. 

Having introduced these factors, [5.26] and [5.27] along with [5.8] and [5.18] 
are used to determine the end moments. 

5.4 PROBLEMS 

5.1 Solve problem 3.1 using Kani's method of moment distribution. 

5.2 Solve problem 3.2 using Kani's method of moment distribution. 

5.3 Solve problem 3.3 using Kani's method of moment distribution. 

5.4 Solve the problem of Example 4.7 using Kani's method of moment 
distribution. 

5.5 Solve the problem of Example 4.8 using Kani's method of moment 
distribution. 
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6. Influence Lines for 
Indeterminate Structures 



6.1 INTRODUCTION 

The determination of the maximum and sometimes the minimum structural 
effects of the appropriate load system is an important preliminary step in the 
analysis and design of structures. Structures subjected to moving or movable 
loads invariably involve the calculation of the maximum or minimum values of 
the bending moment, shear force, deflection and so forth, by preparing the 
influence lines for the various structural effects. 

The influence lines for statically determinate structures may be drawn by 
connecting a few key ordinates with straight lines. However, influence lines for 
indeterminate structures are curved and therefore cannot be drawn so easily. 
The first step in preparing the influence lines for the various functional values is 
to determine the influence lines for the redundants. The next step then is to find 
the influence lines for any other reaction, moment or shear, etc. that can be 
computed by statics. The influence lines for different functions in statically 
indeterminate structures may be obtained using the Miiller-Breslau principle 
backed by computational techniques such as the conjugate beam principle, Cross 
moment distribution, and energy methods. 



6.2 STRUCTURES WITH A SINGLE REDUNDANT REACTION 

The influence lines for indeterminate structures may be constructed by using 
either statical or kinematic methods. When using the statical method considera- 
tion of equilibrium alone is utilised- This may be demonstrated by considering a 
propped cantilevered beam of uniform cross-section. It is desired to prepare an 
influence line for the vertical reaction at support B shown in Fig. 6.1 . 
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h- 




1 




B 



Actual 

structure 



Primary 
structure 






M diagram 



4 



t 




m diagram 



Figure 6.1 

The reaction R B is determined from the compatibility condition 

^B^bb + &bo ~ 0 

from which 



[6.1] 



Rb = ~ 



>bb 



[6.2] 



The deflections b bo and 5 66 are determined from the^Vf andm diagrams (Fig. 

6.1). 

(Mmdx 

1 (x 2 \ (3L-x\ _ x 2 (3L-x) 



EI \ 2 



6EI 
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and 

hb ') EI 

1 (L 2 \ (2L\ L 3 



EI \ 2 J \ 3 / 3EI 
The reaction at B is therefore 
_ 5 6o _ x 2 (3Z,-x) 

This equation gives the value of the reaction i? B when the unit loadP = 1 is 
applied at any position along the beam which is therefore the influence line for 
the reaction at B. Using statics, the influence line for any other reaction, shear or 
moment may be determined. 

Suppose it is required to find the influence line for the shear at the midspan 
of the beam- From statics, the following may be determined: 

(Kc)left = -R* = - X2{ u* X) f ° r ° <X <L ' 2 



(^c)right = 1-*B = 1 X) for L/2<x <L 



2 — forO<x<L/2 



The influence line for the bending moment at the midspan of the beam may 
also be evaluated in a similar manner. 
The bending moment at midspan is 

Mc = Rb~ 

x 2 (3L-x) 
4L 

L I L 

x 2 (3L-x)-4L 2 x + 2L 3 
= for L/2 <x<L 

The ordinates of the influence line at quarter points for the reaction at B, R B , 
the shear at midspan, V c , and the bending moment at midspan, Mq , are given in 
Table 6.1 and also shown in Fig. 6.2. 

To demonstrate the kinematic method of constructing influence lines consider 
the same propped cantilever beam of Fig. 6.1 where the reaction R B is taken as 
redundant. Compatibility condition gives 

* B 5w>+ 5*o=0 [6.3] 
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Table 6. 1 Ordinates of influence lines 



X 


Kb 


V C 


M c 


0 


0 


0 


0 


0.25Z, 


0.0859 


-0.0859 


0.0430L 






[ -0.3125 




0.50L 


0.3125 


or 


0.1563L 






I +0.6875 




0.75L 


0.6328 


+0.2975 


0.0664/. 


L 


1.0 


0 


0 



EI = Consfonf 




I L for R t 



IL for V, 



IL for M, 



Figure 6.2 
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In this case the actual loading on the beam which produces the displacement 
8 bo is the unit load P = 1 at the movable point n, so that the second term can be 
written as 8 bn . Hence 



Using Maxwell's Reciprocal Theorem by replacing the displacement 8 bn by 8 nb , 

*B--f* [6-5] 

The above equation forms the basis for determining the ordinates of influence 
lines since the displacements in both the numerator and denominator are due to 
the same unit vertical load acting at B. The value of 8 nb represents the displace- 
ment in the direction of R B due to a unit load P = 1 travelling along the beam. 
The variation of 8 nb represents the elastic curve of the beam due to a unit load 
at B in the direction of R& Any of the methods of displacement computation 
may be used to find the shape of the elastic curve of the beam under a unit load 
at B. By dividing the ordinate of this curve by a constant factor (— 8 bb ) gives the 
ordinates of the graph representing R B when a unit load P = 1 traverses the 
beam. This curve by definition describes the influence line of R B . 

The elastic curve 8 nb may be obtained using the conjugate beam method as 
shown in Fig. 6.3. 




[6.4] 
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B 



x 



t 
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1.0 



(c) IL for /?, 



Figure 6.3 
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The displacement at any point n due to a unit load at B is 



&nb 



Lx \ (2k \ + (L - x)x ( x 




EI 
x 2 

(3L - x) 



2 \3 2 \3 



6EI 

In a similar manner the displacement at B due to Rq = 1 is 

hb ~ 3EI 

Hence, the equation of the influence line for the reaction at B becomes 

_ 8 nb _ x 2 (3L-x) 
RB " ~8 bb ~ 21? 

The kinematic method may also be used to construct the influence line for 
internal stress resultants such as moments and shears. For the propped cantilever 
beam, the bending moment at midspan may be taken as the redundant reaction 
moment by introducing a hinge as shown in Fig. 6.4 if it is required to draw the 
influence line of the moment at midspan. 

The redundant reaction moment is determined from compatibility conditions 
at the hinge such that the change of rotation of the two continuous sections to 
the right and to the left of the hinge must be zero. Thus 



or 



°CC 



[6.6] 



Therefore, the influence line for the bending moment at midspan will have the 
same shape as the deflection of the beam with a hinge at midspan. The conjugate 
beam method will be sufficient to determine the moment at C, which will be 
equal to 5 CC and the deflections at selected points along the span. For example, 
at point n (0<x<L/2) 

8 n = moment at n of the conjugate beam 




2(L - x)x 2 2x 



2L 3L 

x 2 (3L-x 
EIL 
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Figure 6.4 



8 CC = Reaction at C of the conjugate beam 



■a 



(2)(L) 



I 1 / 1 



4L 
3EI 



The bending moment at point n is 



(d) IL for M c 



M c = 



nc 



'cc 



= x 2 (3L-x) 
4L 2 



forO<JC<L/2 
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Similarly, when x > L/2 

x 2 ( 3L-x \ 41 

The corresponding bending moment is 

w x 2 (3L-x)-4L 2 x + 2L 2 c . ^ ^ r 
M c =— J— 2 (otL/2<x<L 

Notice that the above expressions are identical with those obtained by statical 
method. 

When using the kinematic method to determine the influence line for the 
shear at the midspan of the propped cantilever beam, the beam may be made 
determinate by first cutting at C and inserting a two-bar linkage as shown in 
Fig. 6.5. This linkage system cannot carry shear force and distorts into a 
parallelogram shape thus establishing an equal rotation of the left and right 
tangents at C. 

The displacement of the beam at point n (0 < x <L/2) determined from the 
conjugate beam is 



2 

— ~{3L-x) for0<x<L/2 
obi 



i r (l\ in 





EI N ' \ 2 / V 3 / 3EI 



The shear force when the load is at point n is 

Kc= ir 

_ x 2 (3L -x) 
2L 3 

Similarly, when x > L/2 

The corresponding shear force is 

°cc 

, x 2 (3L-x) , , 
= 1- , — - ioxL!2<x<L 



2L 
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Figure 6.5 

It may be concluded, therefore, that the kinematic method permits a 
simplified approach for the determination of the shape of the influence Une for 
any action. This shape is identical to that of the elastic curve of the corresponding 
primary structure loaded by a unit force or moment at the point the redundant 
has been removed. This analogy, first recognised by MUller-Breslau, provides the 
most widely used and convenient method of computing the influence lines of 
indeterminate structures. The principle of Miiller-Breslau may be stated as 
follows: 



The ordinates of the influence lines of any function (such as reaction, moment, 
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shear) of any structure is represented, to some scale, to those of the displace- 
ment curve which is obtained by removing the restraint corresponding to the 
relevant function imposing in its place a unit distortion such as rotation or a 
linear displacement. 



6.3 INFLUENCE LINES FOR MULTIPLE REDUNDANT STRUCTURES 



6.3.1 Influence Line for Bending Moment 

Applying the Miiller-Breslau principle to construct the influence line for moment 
at any point E between the support B and C of Fig. 6.6(a), a hinge is inserted at 
E so that the moment capacity of the beam is removed without impairing its 
shear capacity. 




/ 



Hinge 





B 




D 



(a) 





/ 





(b) 



6.6 



f 'nE/ S EE 




(0 



The beam is subjected to two force systems as shown in Fig. 6.6(b). A unit 
load is applied at n with the beam deflecting as shown in Fig. 6.6(a). The unit 
load is removed, and equal and opposite moments M are applied on either side 
of the hinge. The deflected shape of the beam is shown in Fig. 6.6(c). The shape 
of the deflected beam is to some scale the influence line of M E . Thus 



M P = 



e 



En 



'nE 



#EE #EE 



[6.7] 
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If 0ee = 1 radian, then 

Thus the influence line for M E is obtained by dividing the ordinates of the 
deflected shape by 6 EE or by setting 0 EE equal to unity. 



[6.8] 



6.3.2 Influence Line for Shear Force 

Let it be required to draw the influence line for shear at point E of the beam of 
Fig. 6.7(a). The beam is cut at E such that the shearing resistance is removed 
without impairing the flexural resistance. The beam is cut at E and a linkage 
system or slide device is inserted which cannot carry shear force and thus which 
permits a relative transverse deflection without introducing a change in slope of 
the left and right tangents at E. 




Figure 6.7 



The beam is subjected to a unit load at n resulting with a vertical displacement 
5 En at E. After removing the unit load, a pair of unit loads are applied as shown 
in Fig. 6.7(b). The shear force at E is given by 



See §ee 
If See = unity, then 



[6.9] 



Vv =6 



nE 



[6.10] 
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Therefore, the influence line for V E is obtained by dividing the ordinates of the 
deflected shape by 6 EE . 

EXAMPLE 6.1 Compute ordinates of the influence line for the moment at B 
of the beam shown in Fig. 6. 8. Use intervals of 2.5 m for span AB and 4.0 m for 
span BC. The moment of inertia is constant. 

Using the Muller-Breslau Principle, the capability of the beam to resist moment 
at the section for which the moment influence line is desired is removed by 
introducing a pin. Unit moments, say in kN m, are applied to the beam at each 
side of the pin at B. The modified beam, which will deflect as indicated by the 
dashed line, is shown in Fig. 6.8(b). According to the Muller-Breslau Principle, 
the various influence line ordinates are computed from the relation 

AM - *>»B 

Mb ~ — 
odd 

Thus, the values of 5 BB , as well as the deflections at the necessary sections of 
the modified beam, may be evaluated using the conjugate beam method. 
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Determination of the reactions (slopes): 



_ 2 / lOx 1 

B = 3 V 2 



+ = 8.667 = 5 BB 



^4(^1=2-667 
Determination of moments (deflections): 



. **i /2.5 x0.25\ /2.5\ 

A*2.5 = 1-667 x 2.5 - ^ J ^—1 = 3.906 

m -iaa7 <n , 5.0x0.50 \ /5.0\ 
M s .q = 1.667 x 5.0- I 1 I — 1 =6.250 

. **i nc 1 7.5 x 0.75 \ / 7.5 \ 
^7.5 = 1.667x 7.5 - I 11 —I =5.469 



M IO = 0 



Ml4 = 2.667 x 12 -(^^)(^)= ,4.00 
M 18 = 2,67x8-(^)g)= 16 .00 

Af 22 = 2.667 x 4 10.00 

The value of the influence line ordinate at each point is determined by dividing 
each moment by 5 BB = 8.667. The resulting influence line is shown in Fig. 6.8(d). 

The method of solution is first to assume that the fixed-end moment of 
100 kN m exists at support B of the member BA with no other fixed-end 
moments being considered to exist. 



6.3.3 Influence Lines by Moment Distribution 

The Cross method of moment distribution may be used to obtain more easily 
influence lines for continuous beams and frames. The method is illustrated by 
the following example. 
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EXAMPLE 6.2 Compute the ordmates of the influence line for the moment at 
B of Example 6. 1 using the moment distribution method. 

The method of solution is first to assume that the fixed-end moment of 
100 kN m exists at support B of the member BA with no other fixed-end 
moments being considered to exist. Using the moment distribution method, the 
moments are obtained from the balancing operation. Next, a fixed moment of 
100 kN m is assumed at end B of member BC and again a similar operation is 
performed. 

The moment distribution is shown in Table 6.2 



Table 6.2 



Joint 


A 


B | 


C 


Member 


AB 


BA 


BC 1 


CB 


K 


7/10 


7/10 


7/16 


7/16 


DF 


1.0 


0.6154 


0.3846 


1.0 


m£ b = l oo 


0 


+ 100 
-61.54 


0 

-38.46 


0 




+38.46 


-38.46 




m£ c = i oo 


0 


0 

-61.54 


+ 100 

-38.46 


0 




-61.54 


+61.54 





After finding the final moments due to the 100 kN m moment at every point 
where fixed-end moments can exist, the equation forM BA in terms of the initial 
fixed-end moments is 

M B A = 0.3846Mb\ + 0.6154MFb>; 

The fixed-end moments for a 1.0 kN load placed successively at each of the 
points for which an influence line ordinate is desired, are computed below. The 
fixed-end moment for a propped cantilever loaded with one concentrated load? 
is 

- Pabl b\ 
where a is measured from the pinned-end. 
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For span AB: 



Table 6.3 



x = a 


b = 10 - a 


M£ B =ab(a + b/2)/L 2 


2.5 


7.5 


1.172 


5.0 


5.0 


1.875 


7.5 


2.5 


1.641 


For span BC: 




Table 6.4 






x = a 


b = 16 -a 


Mg A =ab(a + b/2)L 2 


4.0 


12.0 


1.875 


8.0 


8.0 


3.000 


12.0 


4.0 


2.625 



Note: The ordinate* is measured with support C as the origin. The above values 
of the fixed moments are substituted in the equation above for Mba and the 
influence line ordinates are computed in Table 6.5. 



Table 6.5 



Point, x(m) 


A*BA 


A*BC 


M AB = 0.3846Mf A + 0.61 54m£ c 


0 


0 




0 


2.5 


1.172 




0.450 


5.0 


1.875 




0.721 


7.5 


1.641 




0.631 


10.0 




0 


0 


14.0 




2.625 


1.615 


18.0 




3.000 


1.846 


22.0 




1.875 


1.154 


26.0 









The influence line for the moment at support B is shown in Fig. 6.8(d). 
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6.4 PROBLEMS 

6.1 Draw the influence lines for the beam shown in Fig. P6.1 for the support 
reactions R A and R B . 




Figure P6.1 



6.2 Draw the influence lines for the beam shown in Fig. P6.1 for the moments 
M B and A/ 5 . 

6.3 Draw the influence lines for the beam shown in Fig. P6.1 for the shear S 5 . 
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7. Introduction to Matrix 
Analysis 



7.1 INTRODUCTION 

After the introduction of high-speed computers, there has been a revolution in 
structural analysis, not only in the computational methods but also in the 
fundamental theorems. Since digital computers are ideally suitable for 
automatic computations of matrix algebra, it was found desirable to formulate 
the entire structural analysis in matrix notation. Matrix methods of structural 
analysis are based on the concept of replacing the actual structure by an 
equivalent analytical model consisting of discrete structural elements having 
known properties which can be expressed in matrix form. Matrices are useful 
in expressing structural theory and in producing an efficient means for carrying 
out numerical calculations. 

Two methods have been formulated in matrix structural analysis: the 
flexibility and stiffness methods. It will not be possible in this textbook to 
develop the two matrix methods to sufficient depth. The methods are developed 
to the level of manual computation. 

7.2 FORCE AND DISPLACEMENT MEASUREMENTS 

It is evident that the overall description of the behaviour of a structure is 
accomplished through the dual consideration of force and displacement 
components at designated points. There are a number of ways of measuring a 
force applied to a structure or its displacement at designated points in a 
prescribed direction. Such points are commonly known as node points. The first 
step in the analysis of structures is to idealise the actual structure into a 
mathematical model which consists of distinct structural elements interconnected 
through node points. In this text the word force includes moment. 
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J 



2 



Figure 7. 1 



To designate the forces and displacements at the nodes of a given structure, a 
coordinate system is used to identify these measurements. For the frame shown 
in Fig. 7.1, for example, the system consists of four arbitrary coordinates which 
are identified by four numbered arrows shown at the specific nodes or joints. 
The forces are listed in column matrix [P] and is referred to as a force vector 
and represents an ordered array of force measurements. For instance, the force 
vector for the frame of Fig. 7.1 is represented by 



[P} = 



Pi 
P3 
P4 



[7.1] 



Likewise, the coordinate displacement vector, having the same significance as in 
the force vector may be expressed as 



[A] = 



5i 
5 4 



[7.2] 



In a similar manner, the forces and displacements at the nodes of a given 
element may be designated by listing in column matrices [P] and [A] , 
respectively. For the beam element of Fig. 7.2, for example, with direct forces at 




Figure 7.2 
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node 1 and moments at node 2, the force vector is written as 



[p] = 



xl 

*yi 
M x2 



[7.3] 



and the displacement vector as 



[A] = 



"i 



v 



0 



e 



y2 



[7.4] 



A necessary step in the formation of the force and displacement vectors is the 
establishment of the node points and their location with respect to coordinate 
axes. At this stage it is necessary to define two sets of orthogonal coordinate 
systems. The first set is that of the structure, known as the global axes, and 
consists of a single coordinate system. The second set is that of the members or 
elements, known as the local axes, and consists of one coordinate system for 
each member. Since the members are in general differently oriented within a 
structure, these axes originating at member ends will usually be differently 
oriented from one element to the next. Global and local coordinates are 
illustrated in Fig. 7.3(a) for trusses and in Fig. 73(b) for frames. 

When forces are applied to structures, displacements occur. Alternatively, 
when displacements are prescribed, node forces are necessary to produce them. 
The relationships that exist between applied forces and displacements play an 
important role in structural analysis. The force and displacement characteristics 
of a structure are usually described under definitions of flexibility and stiffness 
coefficients. The flexibility and stiffness coefficients depend on the force- 
displacement properties of the structure and the coordinate system used. 

A simple illustration of such relationships is obtained by considering a linear 
elastic spring shown in Fig. 7.4. Single coordinate is indicated for the force and 
displacement measurements. The force P will stretch the spring thereby 
producing a displacement A at the end of the spring. The relationship between P 
and A can be expressed as 



A=fP 



[7.5] 



In [7.5] , / is the flexibility coefficient of the spring and is defined as the 
value of the displacement at node 1. In general, a flexibility coefficient is the 
value of the displacement at a point of the structure, in a given direction, due to 
a unit force applied at a second point in a second direction. 
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Figure 7.3 




An alternative way is to establish a relationship between the force P and 1 
displacement A for the spring of Fig. 7.4. The force P required to produce a 
displacement A units is determined from 



P = kA 



[7.6] 



In [7.6], k is the stiffness coefficient of the spring and is defined as the value 
of the force required at coordinate 1 to produce a unit displacement at 1. In 
general, a stiffness coefficient is the value of the force at a point of the structure, 
in a given direction, due to unit displacement applied at a second point in a 
second direction. 



f , k 

wvwwwv 




Figure 7.4 
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Comparison of [7.5] and [7.6] reveals that the flexibility and the stiffness 
of the spring are inverse to one another. 



[7.7] 



Now consider a more general case consisting of an elastic structure, supported 
against rigid-body motion, and subjected to loads P 1 ,P 2 , . . . ,P n acting at 
nodes 1, 2, . . . , n. The corresponding set of displacements is represented by L x , 
A 2 , . . . , A„. For linearly elastic systems, the principle of superposition is 
applicable. Therefore, the displacement A, at node / is given by 



or more generally, 

j = n 

A, = £ W 



[7.8] 



[7.9] 



By definition, fy is the displacement produced at node i due to a unit load at 
node / (Pj = 1). The coefficients //,-, which are the displacements due to unit 
loads, are known as flexibility coefficients. 

In general, for n nodes, there will be n such displacements which may be 
written in a single matrix equation 



A, 

A 2 



f\\ f\2 

fl 1 /22 



fin 

fin 



fn 1 fnl 



f, 



nn 



p? 



[7.10] 



and which can be written in compact matrix form as 
[A] = [F] [P] 



[7.11] 



where [A] is the column displacement matrix, [F] is a square flexibility matrix 
and [P] is a column load matrix (load vector). This equation is of the same type 
as [2.17]. 

Using matrix operation, one can solve the set of algebraic equations 
represented in [7.10] for forces in terms of displacements. In matrix notation 



[P] =[F]" 1 [A] 



[7.12] 
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where [F]~ l is the inverse of matrix [F]. It is noted that [7.12] has the 
same form as [7.6] since it expresses forces in terms of displacements. 
Consequently, 

[F]' 1 = [K] [7.13] 

where [K] is the stiffness matrix which is the inverse of the flexibility matrix. 
Thus 



[P] = [K] [A] [7.14] 
The expanded form of [7.14] is 




[7.15] 



By definition, k xy is the force required at node / to produce a unit displacement 
at node / only (zero displacements at all other nodes). 

Flexibility coefficients for linear elastic behaviour have the property of 
reciprocity which may be expressed analytically as 



ftr fn 



[7.16] 



This equation defines symmetry of [F] . Since [F] is symmetrical the inverse of a 
symmetric matrix will also become symmetrical. Therefore, [7.13] guarantees 
that the stiffness matrix [K] will likewise be symmetrical. Consequently, 



- kji 



[7.17] 



To illustrate these matrices consider a simple cantilever beam of uniform cross 
section shown in Fig. 7.5(a). To determine the flexibility matrix, the influence 
coefficients must be determined by applying unit loads to the free end. 



Due to axial load N = 1 (Fig. 7.5(b)) 




5vn=0 [7.18] 
6 n =0 
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EA, El 



I 




N 



N = l 





Figure 7.5 



Due to vertical load V = 1 (Fig. 7.5(c)) 
6„=0 

5 -ii 
vv 3EI 



6v= 2EI 



Due to moment M = 1 (Fig. 7.5(d)) 



5„=0 



The above results may be written in matrix form as 



'n 
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EA 
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or, when written in compact matrix form 
[A] = [F] [P] 

In a similar manner the stiffness matrix may be determined by unit 
displacements as shown in Fig. 7.6. 



[7.11] 



1 



EA, El 




N 



(a) 



N 



(b) 




Ic) 




Figure 7.6 



Due to unit axial displacement (Fig. 7.6(b)) 



(d) 



N = 



EA 



[7.22] 



Due to unit vertical displacement (Fig. 7.6(c)) 



V = 



\2EI 



[7.23] 



M=- 



6£7 
L 2 



Due to unit rotation (Fig. 7.6(d)) 



V=- 



6EI 



M = 



4EI 



[7.24] 
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The above results may be written in matrix form as 



TV 



M 



EA 



0 



0 - 
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\2EI 


6EI 


L 3 


L 2 


6EI 


4EI 


~ L 2 


L 



6„ 



V 



e 



[7.25] 



which may be written in compact matrix form as 

[P] = [K] [A] 
The results may be checked by matrix multiplication 

[F] [K] = [K] [F] = [/] 

It is noted that both [F] and [K] are symmetric matrices which is the 
consequence of the reciprocal theorem. 



[7.26] 



[7.27] 



7.3 THE FLEXIBILITY METHOD 

The basic theory of the flexibility method is developed in this section, and the 
concepts are clarified by numerical examples. The development of the method 
rests on the basic principles of equilibrium of forces, compatibility and linear 
force-displacement relationships. 

Consider a structure, which is idealised into a model consisting of distinct 
structural elements interconnected through node points, under the action of 
generalised external forces appUed at the nodes P lt P 2 , . . .P„. These may be 
conveniently represented by a column matrix or force vector [P] 

[P] =\P u P 2 ...P n \ [7.28] 

Let it be assumed that the structure consists of m redundants which are 
forces to be determined, that is 

[X\={X l9 X 2 ...X m \ [7.29] 

which are the redundant forces or reactions. If such redundants are removed, 
the structure becomes determinate and the internal forces are determined from 
conditions of equilibrium alone. In an indeterminate structure, the internal 
forces must also satisfy compatibility in addition to equilibrium. In dealing with 
an indeterminate structure with m redundants, the redundants are treated as 
additional loads on the statically determinate structure. It is assumed that the 
structure is composed of an assemblage of / simple elements. Internal forces 
exist in the structure at the node points. If the internal force members are 
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represented by the vector [S] where 

[S] = \S 1 S 2 ... S f \ 
then, [S] can be related to the applied loads [P] and [X] as 



[7.30] 







r pi i 
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= [Bo] 


• 

• 
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[7.31(a)] 



which is written in compact notation as 

[5] = [B 0 ][P] + [B 1 ]m 
or using partitioned matrices 



[S] = [B 0 \B 1 ] 



[7.31(b)] 




[7.31(c)] 



where, in general, [B 0 ] and [fi, ] are rectangular matrices whose elements are 
obtained from equilibrium conditions of the structure. For example, if P t is 
taken as a unit load with all other loads including [X] held at zero, the internal 
forces in the structure represent the coefficients corresponding to the jth column 
in the [B 0 ] matrix. Likewise, the internal forces which result from a unit load Xj 
with all others held at zero represent the coefficients corresponding to the yth 
column of the [B x ] matrix. 

To formulate the compatibility condition, the principle of least work will be 
utilised which may be stated as: The true values of the redundant forces are 
those which make the strain energy U of the strained structure a minimum. 

The strain energy is given as 



5, 

S 2 



[7.32] 



which is written in compact matrix form as 



1 



U=~ [S] T [F] \S\ 
In order to obtain the strain energy U in terms of the unknown \X\, 



[7.33] 
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substitute [7.31(b)] into [7.33] 
be written from [7.31(b)] as 



In doing so, note that the transpose \S | T may 



[S] 



= [y T [*o i *, 



= [P>,X}[B 0 \B i ) 



[7.34] 



Substituting [7.31(b)] and [7.34] into [7.33] 



1 



U = -[P\X)[H] 




[7.35] 



where 



m^[B 0 \B 1 ) T [F\[B 0 \B l ] 

Since [P] and [X] are the applied and redundant forces, respectively, it is 
convenient to partition [H] to conform to the load vectors, thus, 



[7.36] 



U = \[P\X] 



n pp . 

i 



px 



H 



xp . 



H 



XX 



[7.37] 



After expanding [7.37] 



1 



U = ~ ([P] [H pp ] [P] + [P] [H px ] [X] 
+ [X)[H xp ][P) + [X][H xx )[X)) 



[7.38] 



Utilising the theorem of least work and noting that the [H] matrix is symmetric 
gives 



W 
dX 



= W xp ] [P] + [H xx ] [X] = 0 



[7.39] 



from which the redundants are determined as 



[X]=-[H xx )- l [H xp ][P] 



[7.40] 



Solving for [X] from [7.40] all internal forces can be determined from 
[7.31]. 

Summarising, the essential steps in applying the flexibility method to lead to 
the solution of structural problems may be stated as follows : 

1. Idealise the structural problem to be analysed 

2. Specify the redundant forces and identify the internal member forces 
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3. Find the [B 0 ] matrix for unit values of external forces; only one external 
force must act at a time with all other forces held at zero 

4. Find the [flj matrix for unit values of redundant forces; only one 
redundant force must act at a time with all other forces held at zero 

5 . Find the flexibility matrix [F] for all members following the sequential 
order of the member forces in [B 0 ] and [B , ] 

6. Formulate the [H] matrix of [7.36] 

7. Calculate the redundant forces [X] from [7.40] 

8. Calculate the internal forces [S] from [7.31] 



EXAMPLE 7. 1 Using the flexibility matrix method determine the bar forces 
in the truss with double diagonal system shown in Fig. 7. 7(a). The area of all 
top chord is twice the area of all remaining members. 





(b) 



Figure 7.7 



As can be easily seen, the truss is redundant to the second degree. For the 
selection of the redundant members several choices exist. Here members AE and 
BE and the reaction at A are taken as the redundants, then the truss is reduced 

to a determinate one as shown in Fig. 7.7(b). 

To determine the [B 0 ] matrix, P x and P 2 are set unit values one at a time 
with all other forces including the redundants held at zero, then calculate the 
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internal forces in all members for each case. Thus, 



[B 0 ] = 



P x = \ 
0 

-0.667 
-0.500 
0 

+0.667 
0 

-0.833 

+0.833 

-0.500 

0 
0 



P 2 = l 
0 

+0.500 
-0.375 
0 

0.500 
0 

+0.625 
+0.625 
-0.375 

0 

0 



Member 

AB 

BC 

CF 

FE 

ED 

DA 

DB 

CE 

BE 

AE 

BF 



Similarly to determine the [2? J matrix the redundants X x , X 2 and X 3 are 
set unit values one at a time with all other forces including the applied loads 
held at zero. The internal forces in each case are 



[*i] = 



X,=l 


x 2 = 1 


*3 = 1 
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The flexibility matrix for the members is 



[F] = 



1 

EA 
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From [7.361 

[H] = [B 0 \B 1 ] T [F][B 0 \B 1 ] 

substituting and carrying out the matrix multiplications gives 



[H} = 



EA 



11.111 

1.792 


1.792 : 
6.25 




7.031 
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XX 



The redundants are determined using [7.40] 
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-3.800 -5.275 
Solving for the redundants, 
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[ 4.617] 
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The bar forces are determined using [7.31 ] 
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EXAMPLE 7 .2 Determine the shear force and bending moment values in the 
continuous beam of Fig. 7.8 using the flexibility matrix method. 



I 



V 



5m 



F> = IOkN 
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P 2 =5kN 
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(a) 



(b) 



10 m 



Figure 7.8 

The beam is indeterminate to the first degree and the reaction at mode 5 is 
chosen as the redundant as shown in Fig. 7.8(b). 

To determine the [B 0 ] matrix, P, and P 2 are set unit values one at a time, 
thus 
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Similarly the [/?, ] matrix is determined setting X x = i ? thus 
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The member flexibility matrix from [7.20] neglecting the axial deformation 



is 



[F} = 



0.833 
EI 



50 
15 




50 
15 




50 
15 




Using [7.36] to solve for the [H\ matrix 
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The redundant is determined using [7.40]. Hence, 
* 1= ~ 800 175 " 325] 



10 




50 




= 1.09375 kN 

The shear force and bending moment values at the indicated nodes are 
obtained using [7.31]. 
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7 A THE STIFFNESS METHOD 



As in the flexibility method, the stiffness method considers a structure as an 
assemblage of individual members. The connecting members are called node 
points. The fundamental difference in using the stiffness method is that the 
individual displacements of the nodes are taken as the unknowns. In the stiffness 
method the number of unknowns to be determined is the same as the degree of 
kinematic indeterminateness. 

In this method, the first step is to derive the stiffness matrix for a 
component member by relating the member forces to the member deformations. 
In a similar manner the nodal forces must be related to the nodal displacements 
by the total stiffness matrix obtained from an assemblage of the stiffness matrices 
of the individual members. Finally, from equilibrium conditions the nodal forces 
obtained from the unknown nodal displacements must balance the externally 
applied nodal forces to find the total solution; that is, determining all unknown 
displacements, reactions and member forces. In developing the stiffness method, 
the same coordinate systems are used that were employed in the flexibility 
method. 



Member Stiffnesses 

The relationship between the forces acting at the nodes P t - and their corresponding 
nodal displacements forms the stiffness matrix approach. This relationship is given 
in matrix notation by [7.14] and in its generalised form by [7.15], 

Consider a prismatic axial rod element m the ends of which are denoted as i 
and / as shown in Fig. 7.9. 



p. & ' (EA/ m 



a) 



A' 



b) 



L 



L 



Figure 7.9 



The relationship between the axial forces and the corresponding displacements 
of the rod is 



[7.41] 
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The coefficients of the stiffness matrix are found by considering two distinct 
displacement states. The first state is to let the nodal coordinate displacement 
5/ = 1 as shown in Fig. 7.9(b) while holding the others at zero. Imposing 
equilibrium on the forces gives 



« ■ A") 



[7.42] 



m 



The second state is similar, but distinct from the first. Following the same 
procedure as in the first state provides 



[7.43] 



m 



Combining the results given by [7.42] and [7.43] into a single matrix equation 
yields the force-displacement relationship of an axial rod element 
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[7.44] 



Consider a prismatic beam element shown in Fig. 7.10. Using the same 
procedure used in obtaining [7.23] to [7.25] the force-displacement 
relationship for the given nodal coordinate system may be determined by 
assigning unit values to the displacements as shown in Fig. 7.10(b) and (c). 
The coefficients are shown for unit displacements at the end i; similar 
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Figure 7.10 
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coefficients are obtained for unit displacements at end /. Thus 
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[7.45] 



Transformation Matrices 

If the properties of an element is known in terms of local axes, the transforma- 
tions of these forces and displacements to the global coordinates is a necessary 
step in stiffness matrix formulation. Figure 7.1 1 shows member if described by 
the two coordinate systems. The local coordinates are shown asx' and;/ and 
global coordinates as x andy. In this text, the forces displacements and stiffness 
matrices with respect to local axes are identified by primes. The prime is omitted 
when written with respect to global axes. 




(a) Local 



(b) Global 



Figure 7.11 



Referring to Fig. 7.1 1, the relationship between the quantities in the local 
and global axes for flexural members is established as 
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[7.46] 
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and for axial members, after omitting M and 0, the relationship is 
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[7.47] 



Equations [7.46] and [7.47] may be written in compact matrix form 
[A = [71 [P] 



[7.48] 



where [T] is the rotational transformation matrix, which is a function of the 
direction cosines between the two sets of axes, for the particular system shown. 
Solving [7.48] for [P] 



[P] = [T- l ][P') 
= [T] T [P'] 



[7.49] 



Such a matrix is called an orthogonal matrix , which may be defined as a square 
matrix having an inverse equal to its transpose. 

If the displacements are denoted by [5 ] then it follows that 



[5*] = m [5] 



[7.50] 



The transformation matrix [T] may be applied to obtain the stiffness matrix 
in global coordinates. From the definition of stiffness, that is [P] = [k] [5 J, it 
follows that 

[P') = [k'][d') [7.51] 

Substituting [7.48] and [7.50] into [7.51] and noting that T T = T~ l for 
orthogonal matrices 

[T] [P] = [k'\ [T] [5] 

or 



[P] = [T- i ][k'][T][8] 
= [T] T [k'} [71 [6] 

= [*] [6] 

Hence, the transformed stiffness matrix is given by 
[k] = [T] T [k') [T] 



[7.52] 



[7.53] 



Using the relationship derived above, the stiffness matrix for axial members 
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(Fig. 7.1 1) in global coordinates will be 
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where X = cos a and /i = sin a. Similarly, for flexural members 
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[7.54] 



12m 2 








-12X|i 


12X 2 




Symmetric 




-6LX 


4L 2 j 




-12/i 2 


12X/i 


-6L/i ; 


12/i 2 


12X/i 


-12X 2 


6LX j 


-12X/i 12X 2 


6L/i 


-6LX 


2L 2 j 


-6L/i 6LX 41 2 



[7.55] 



Assembly of Element Matrices 

It is important to form the total assemblage nodal stiffness matrix of a structure 
from the stiffness matrices of the separate structural elements. This involves 
only simple additions when all element stiffness matrices have been expressed in 
the same global coordinate system. 

Consider the axial member system shown in Fig. 7.12 with a total of three 
possible joint displacements one for each node. The members have individual 
stiffness constants (EA/L)i and (EA/L) 2 as shown in the figure. 

0 m 



«-2 



Figure 7.12 

The order of the stiffness matrix for the assemblage will be 3 x 3. The individual 
member stiffness matrices are: 



«-(?), U 






[7.56] 
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The assembled stiffness matrix for the complete system can be formed by 
superposition of the individual element stiffnesses contributing to each nodal 
point. Thus, the assembled stiffness matrix for the system shown in Fig. 7.12 
becomes 



[K] = 




1 



[7.57] 



In a similar manner, it may be concluded that the order of the global stiffness 
matrix of a system is equal to the total number of degrees of freedom of the 
system. The order of the matrix may be expressed as the sum of the unknown 
displacements /and the prescribed (support) displacements s. After reordering 
the rows and columns to separate the elements corresponding to the supports 
from the remainder, the rearranged stiffness matrix may be written as 



UJ lK sf j K„\ |_Aj 



[7.58] 



Method of Solution 

Expanding [7.58] and noting that the support displacements, \A S \ =0 

[P f ] = [K ff ] [A,] [7.59(a)] 

[P s ] = [K sf ] [A f ] [7.59(b)] 

The vectors of all unknown nodal displacements (at unsupported nodes) are 
obtained from [7.59(a)] 

[A f ] = [%]"• [P f ] [7.60] 

When [A^] has been found from [7.60] , the support reactions by substitution 
of the results in [7.59(b)] will be 

[P s ] = [K sf ][K ff ]- l [P f ] [7.61] 

The internal force in any element m may be obtained by substituting the 
calculated degrees of freedom for that element, designated by [A m ], into the 
element stiffness matrix [k m ]. Thus, the joint force component acting on that 
element becomes 

\P m ] = [k m ] [A m ] [7.62] 
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For the case of an axial member shown in Fig. 7,1 1 the member force, 
denoted by P m , is found to be 



P m = P xj cos a + P y j sin oc 



[7.63] 



but 



Pxj = ~[^ b xj ~ 5 x /)cos^q; + {b yj - 8 yi )cos a sin a) 



EA 



P yj = — ((5*,- - 5 x/ )cos a sin 0^(6 w - - 8 y /)sin a a) 
Rearranging and writing in matrix form, [7.64] may be written as 



EA 



[7.64] 



m 



(EA\ ["**/- 
= I — [cos a sin a] 



[7.65] 



Similarly, for the case of a beam element [7.45] and the true nodal 
displacements are used and the internal forces for the beam element taken to be 
the shear and bending moments are 



EI 
3 



m 




-6L -12 -61 



4L 



6L 2L 



Oi 



6y< 



& J 



[7.661 



When these are external loads acting between the joints of a beam element the 
concept of equivalent loads must be adopted. The member action is then 
computed by adding the effects of the member end deformation to the fixed-end 
actions produced by the loads. In a similar manner, the support reactions are 
computed by adding the fixed-end effects of the loads. Thus [7.14] may be 
written as 

[P] = [K] [A] - [1*] 



[7.67] 



where [Z^] is the load vector of the fixed-end actions. 



EXAMPLE 7.3 Determine the bar forces, using the stiffness matrix method, of 
the truss shown in Fig. 7. 13. EA= constant. 



Member data for the truss 



Member 


Member ends 


Member properties 




n cosines 


■ 

1 










sin a 


/ 


A 


L 


cos a 


1 


A 


B 


A 


L 


0 


1 


2 


A 


C 


A 


1.155L 


-0.5 


-0.866 


3 


A 


D 


A 


1 .4 1 46L 


0.707 


-0.707 
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(a) (b) 

Figure 7.13 



The member stiffnesses oriented in global coordinate system are obtained from 
[7.54]. 

Member AB 



[*i] = 



EA 



0 


0 


0 


0 1 


1 


0 


1 


0 


-1 


2 


0 


0 


0 


0 


3 


0 


-1 


0 


0 


4 



Member AC 



[* 2 ] = 



EA 



1.155L 



*3 = 



EA 



1.414Z, 



0.25 


0.433 - 


-0.25 


-0.433 


1 


0.433 


0.75 


-0.433 


-0.75 


2 


-0.25 


-0.433 


0.25 


0.433 


5 


-0.433 


-0.75 


0.433 


0.75 


6 


0.5 


-0.5 -0.5 


0.5 1 


1 




-0.5 


0.5 0.5 


-0.5 


2 




-0.5 


0.5 0.5 


-0.5 


7 




0.5 


-0.5 -0.5 


0.5 . 


8 





After assemblage of the element stiffness matrices, the global stiffness 
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equation becomes 



r p ~\ 




0.570 


0.021 i 


0 


0 


-0.217 


-0.375 - 


0.354 


0.354"] 




A,l 

1 


-Q 




0.021 


2.00 i 


1 o 


-1.0 


-0.375 


-0.650 


0.354 


-0.354 




A 2 


0 




0 


0 ! 


— 

1 0 


0 


0 


0 


0 


0 




A 3 


0 


_EA 


0 


-1.0 


! 0 


1.0 


0 


0 


0 


0 




A 4 


0 


L 


-0.217 


-0.375 


: o 


0 


0.217 


0.375 


0 


0 




A s 


0 




-0.375 


-0.650 


1 o 


0 


0.375 


0.650 


0 


0 




A 6 


0 




-0.354 


0.354 


; 0 


0 


0 


0 


0.354 


-0.354 




A 7 


_ 0_ 




L 0.354 


-0.354 


: o 


0 


0 


0 


-0.354 


0.354 







Note that the displacements A 3 to A 8 are restrained and the elements 
corresponding to A, and A 2 are placed at the top left of the global stiffness 
matrix. Hence the reduced stiffness matrix given by [7.59(a)] is written as 



570 
021 



0-021-1 TAH 
2.00 J LaJ 



The unknown displacements are obtained by applying [7.60], 




L_ 
EA 




755 -0.018] 
018 



0.500J 




Having found the unknown displacements, the reactions are calculated by 
applying [7.61]. Thus 



Pcy 
P Dy - 



r 0 


0 "I 


0 


-1.0 


-0.217 


-0.375 


-0.375 


-0.650 


-0.354 


0.354 


■ 0.354 


-0.354 - 


r 0 


0 -1 


0.018 


-0.500 


-0.373 


-0.184 


-0.646 


0.318 


-0.628 


0.183 


- 0.628 


-0.183 - 




755 - 



0.0181 V PI 

0.500 J L-eJ 



-Q 
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The internal member forces are obtained by using Eq. [7.65]. Thus 



D EA TA xi -A x A 
P m = -T [cos a sin a] 

L LA^y - A yi J 



Since A xj = A yj = 0 because of support conditions, the member forces are: 



Member AB 



'AB 



[-X.1SSP- 0.018(21 

= [0 1] 

L 0.018P+ 0.50(2 J 



= 0.0 18/> + 0.50(2 



Member AC 



1 



1.155 



[_0.5 -0.866] 



r"-1.755/>- 0.01801 
L 0.018P + 0. 



50(2 J 



= 0.746/> + 0.367(2 



Member AD 



AD 



i r-1.755P- 0.018(21 

= 7777 [0.707 -0.707] 
1-414 L 0.018(2 + 0.500 J 



= -0.887P - 0.259(2 



EXAMPLE 7.4 Find the support rotations and the support reactions of the 
continuous beam shown in Fig. 7.14. 

The equivalent fixed-end actions are shown in Fig. 7.14(b). 

The member stiffnesses are (see [7.45]): 

Member AB 



[*i] = 





12 


-48 


-12 


-48 


1 


EI 


-48 


256 


48 


128 


2 


512 


-12 


48 


12 


48 


3 




.^8 


128 


48 


256 . 


4 
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2kN/m 



IkN/m 



B 



8m 



8kN 



C 1 



(a) Actual Structure 

12 kN 

5.33 kN 




l0.67KN-m 



(b) Equivalent Joint Load 




G 



CD 




4 



1 



(c) Coordinate System 

Figure 7.14 



Member BC 





r 12 


-48 


-12 


-481 


3 








-48 


256 


48 


128 


4 








-12 


48 


12 


48 


5 








-48 


128 


48 


256 _ 


6 






The assembled stiffness matrix is 














12 


-48 


-12 


-48 


0 


0 "1 


1 




-48 


256 


48 


128 


0 


0 


2 




















-12 


48 


24 


0 


-12 


-48 


3 




-48 


128 


0 


512 


48 


128 


4 




0 


0 


-12 


48 


12 


48 


5 




_ 0 


0 


-48 


128 


48 


256 _ 


6 



Note that the prescribed displacements are - 5 2 - §5 - 0 and 6 4 and 5 6 
are the unknown displacements. Thus the assembled stiffness matrix must be 
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rearranged by moving the rows and columns associated to 8 4 and 5 6 to the top 
left of the stiffness matrix. Thus 



[K) = 



EI 
512 



512 


128 


-48 


128 


0 


48 


4 


128 


256 


0 


0 


-48 


48 


6 


-48 


0 


12 


-48 


-12 


0 


1 


128 


0 


-48 


256 


48 


0 


2 


0 


-48 


-12 


48 


24 


12 


3 


_ 48 


48 


0 


0 


-12 


12 _ 


5 


1 

f • 

— i— 














r . 


*J 













The submatrix [K/f] is 



[%] ■ 



_£7_T512 128 
512 L 128 256 



and its inverse is 



[%] 



-1 



=-L_[ 

224£Y [_ 



256 -128 



128 512 

The load vector for the beam shown in Fig. 7.14(b) is 



[P] = 



r -8 -| 


1 


10.67 


2 


-12 


3 


-5.33 


4 


-4 


5 


--5.33- 


6 



Rearranging the vector in conformity with the stiffness matrix is 



[p} = 



-5.33 


4 


-5.33 


6 


-8 


1 


10.67 


2 


-12 


3 


-4 


5 
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The unknown rotations at support B and C are computed using [7.60] 



[A/] 



224£Y |_- 



256 -128 
128 512 



I r-5.331 = j_ r-3.051 
J L-5.33J ~ El L9.I4J 



The support reactions due to the displacement contributions are determined 
by using [7.59(b)]. After adding the fixed-end effects of the loads, the support 
reactions will be 

[P 5 )=K 5f [A f }-[P F ] 



\ RA ] 




r-48 


0 "I 


M A 


EI 


128 


0 




512 


0 


-48 


.Rc - 




. 48 


48 . 



1 



EI 



r-3.os-| 

L-9.14J 



r -8 "I 




r 8.29 1 


10.67 




-11.43 


-12 




12.85 


. -4 - 




. 2.86. 



7.5 PROBLEMS 



7.1 Using the flexibility method, find the internal forces on all members for the 
truss shown in Fig. P7.1 . The cross sectional area of the members, in square 
centimetres are shown in parentheses. 

(Ans: F AB = 1 .37 kN F BC = 5.33kN F CE = 5.94kN 
F AC = -1 .6 kN F EF = -1 .98 kN F FE = -7.46 kN 
F BE = -1 1 .47 kN F BF = -4.95 kN F EC = 0.77 kN 
F CF = -1.20 kN) 



3m 



12 kN 



8kN 




4m 



4m 



4m 



Figure P7.1 
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7.2 Find the internal forces on all members of the truss shown in Fig. P7.2 using 
the flexibility method. Assume EA to be constant. 




7.3 Using the flexibility method, determine the support moments of the 
continuous beam shown in Fig. P7.3. Assume EI to be constant. 

(Am: M A = -10.8kNm A/ B = -5.4kNm M c = +2.7 kN m) 



12 kN 



ZZl 



6m _^ 4m 



Figure P7.3 



7.4 Find the internal forces of the frame shown in Fig. P7A Use the flexibility 
method. Assume EI to be constant. 
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7 .5 Using the stiffness method, find the internal forces in all members for the 
truss shown in Fig. P7.5. Assume EA to be constant. 

(Arts: F AC =6.22kN F CD = 4.67 kN F AD = 8.90 kN 
F BC = -7.78 kN F BD = -7.10 kN) 



lOkN 



4m 




Figure P7.5 



7.6 Solve Prob. 7.2 using the stiffness method. 

7.7 Determine the support moments for the continuous beam in Fig. P7.7 
Assume EI to be constant. 

(Am: M A = -1 1 .1 1 kN m M B = -4.77 kN m) 



3m 



12 kN 




B 



t 



Figure P7.7 



I 



4m 



7.8 Find the joint moments and shears of the frame shown in Fig. P7.8. Neglect 
axial deformations. 



IkN/m 





B 










2 I \ 




5kN 


I 




1.5 m 


A 







'////)///, 



6 m 



Figure P7.8 
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7.9 Determine the joint moments of the frame shown in Fig. P7.9. Neglect 
axial deformations. 

(Am: M A = 2.59 kNm M B = -5.19 kN m Mq B = -8.83 kN m 
M CD = 2.76 kN m M CE = -6.07 kN m M DC = 1 .39 kN m) 



3m 



3m 



///////A £ 



B 



IkN/m 



21 



BkN 



1.51 



2>. 



WWW 



I 



6 m 



4m 



-j 



Figure P7.9 
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